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PREFACE 

The object of this book is to provide in a compact form an 
account of the methods of integrating explicitly the 
commoner types of ordinary differential equation, and in 
particular those equations that arise from problems in 
geometry and applied mathematics. It takes the existence 
of solutions for granted ; the reader who desires to look into 
the theoretical background of the methods here outlined 
will find what he seeks in my larger treatise Ordinary 
Differential Equations (Longmans, Green & Co., Ltd., 1927). 
With this qualification, it will be found to contain all the 
material needed by students in our Universities who do not 
specialize in differential equations, as well as by students 
of mathematical physics and technology. 

As one of the first things a beginner has to learn is to 
identify the type to which a given equation belongs, the 
examples for solution have not been printed after the 
sections to which they refer, but have been collected at the 
end of the book. When the contents of the first chapter 
have been mastered, the reader may test his skill by 
attacking examples selected at random from Nos. 1 to 122, 
and similarly for the later chapters. The examples occur 
roughly in the order of the table of contents, so that 
working material is always available as reading progresses. 

In conclusion, I wish to record my thanks to the General 
Editors for their encouragement and help during its growth 
and passage through the press. 

May 1939 E. L. I. 

The Second Edition, through the untimely death of the 
author in March 1941, has not had the advantage of being 
revised by him. The Editors are indebted to Dr I. M. II. 
Etherington and Miss N. Walls for some corrections, and to 
Dr A. Erddlyi for undertaking a scrutiny of the book and 
recasting parts of Chapter VI. 


April 1943 


A. C. A. 
D. E. K. 
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CHAPTER I 


EQUATIONS OF THE FIRST ORDER AND DEGREE 

1. Definitions. Let x be an independent, and y a de¬ 
pendent variable; let y\ y ", . . represent successive 

derivatives of y with respect to x. Then any relation of 
equality which involves at least one of these derivatives is 
said to be an ordinary differential equation. The term 
ordinary distinguishes it from a ‘partial differential equation, 
which would involve two or more independent variables, a 
dependent variable, and the corresponding partial deriva¬ 
tives. The order of any differential equation is the order 
of the highest derivative involved. Thus any relation of 
the form 

F(x,y, y\y\ . . y<">)=0 

is an ordinary differential equation of order n. 

Differential equations, both ordinary and partial, are 
of frequent occurrence in mechanics and mathematical 
physics, but the illustrations that best serve to introduce 
the subject are taken from the geometry of plane curves. 

The equation 

f(x,y,C)= 0, . . . (1-1) 

in which x and y are rectangular co-ordinates and C is a 
parameter or arbitrary constant, represents a family of 
curves, in which one curve corresponds to one value of C, 
another curve to another value. If, regarding C for the 
moment as fixed, we differentiate with respect to x, wo 
obtain 8f 8f 
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. ( 1 . 2 ) 
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Generally Bpeaking, (1.2) will involve C ; if C is elimi¬ 
nated between the two equations, there will result an 
equation involving x, y, and y', say 

F(x,y,y , ) = 0, . . . (1.3) 

that is, an ordinary differential equation of the first order . 
When such an equation is polynomial in y (but not 
necessarily in x and y) the index of the highest power of 
y' involved is said to be the degree of the equation. 

Geometrically, the differential equation (1.3) implies 
that at any chosen point of the (x, y)-plane the derivative 
has a certain value or values, that is to say it symbolizes 
a property of the gradient of any curve of the family (1.1) 
that passes through the point (x, y) considered. 

Example 1. The equation 

y=x* + C 

represents a family of equal parabolas having the y-axis ns 
their common axis. On differentiating with respect to x 
we have 

!/' = 2x. 

The arbitrary constant C has disappeared, so that this is 
actually the differential equation of the family of parabolas. 
It expresses the fact that all the curves of the family have the 
same gradient at the points whore they aro cut by a line 
parallel to the y-axis, namely a gradient equal to twice the 
abscissa of tho lino. 

Example 2. The equation 

y = Cx* 

represents a family of similar parabolas having the y-axis as 
their common axis, and all touching the x-axis at the origin. 
Differentiating, wo obtain 

y' = 2Cx, 

which involves C; if this constant is eliminated we obtain the 
differential equation 

y' = 2 y/x. 
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which expresses the property that any line y =mx through the 
origin intersects all curves of the family m point* where they 
have the same gradient 2m. 


2. Integration. The process of elimination by which 
the differential equation (1.3) was obtained from the 
primitive equation (1.1) is in general not reversible; the 
action of recovering the primitive, or an equivalent 
expression, is known as integration. More precisely, to 
integrate or solve a differential equation of the first order 
is to determine all the relations/(x, y)=0 such that the 
values of y and y' deduced from them in terms of * shall 

satisfy the differential equation identically. 

When an infinite set of such integrals can be grouped 
in one comprehensive formula, involving an arbitrary 

constant, say 

/(*, y . C)=o, 

it is known as a general integral; it is in fact either the 
primitive or an expression equivalent to it.* Any Integra 
that can be obtained by assigning a definite numerical 
value to C is a particular integral. But there may be 
integrals other than those that can be obtained by assign¬ 
ing particular values to C\ these are singular integrals. 


As an example, the equation 

(!/')* -*y' +y =° 

(which is of the first order and second degree) possesses the 
general integral y=Cx-C'. This represents a family of 
straight lines, and any particular integral corresponds to a 
definite line in the family. But the equation is satisfied also 
by y = {x\ which represents not a straight line but a parabola. 

This is a singular integral. 


• General theory proves that a differential °T iat j™ ° ( [ th ® 1 
order has one and only one distinct general integral. If two inUyr *fd 
exist, each of which involves an arbitrary constant, they can bo 
tranflfonnod into one another. Soe § 3, hz. i. 
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The remainder of this chapter will be devoted to 
equations of the first order and first degree; that is, to 
equations that may be written in the alternative forms 

p {x,y)+Q{x,y)y'= 0, . . (2.1) 

P(x,y)dx+Q(x,y)dy= 0, . . (2.2) 

where P and Q do not involve the derivative y'. 

Since from any primitive equation of the form 

u(x, y) =C, . . . (2.3) 

where C is an arbitrary constant, we deduce that 

du cu 

—dx + —dy = 0, 
dx dy J 

it follows that a necessary condition for (2.3) to be an 
integral of (2.2) is 

du du 

P dy- Q di=°- ■ ■ ■ < 2 ' 4 > 


3. Separation of Variables. Among equations of the 
typo Pdx + Qdy -0 the simplest are those in which P is a 
function of x alone and Q of y alone, say P = M{x), Q = N(y). 

1 ho general integral is obtained by direct integration, 
thus: 


!M(x)dx+fN(y)dy=C, 

where C is the arbitrary constant of integration. 

More generally, let P and Q be the products of a term 
involving x but not y and a term involving y but not x, 
so that the equation has the form 


M (x)R(y)dx + N(y)S(x)dy= 0. 

The variables are then said to be separable, for 
h y p (y)S{x) the equation becomes 


on division 
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and has the general integral 



When, as in the above cases, the process leads to an 
expression that involves integral signs the result is said 
to be an integration by quadratures. This implies that the 
problem has been reduced from one in differential equations 
to an equivalent in the integral calculus. If it is found 
impossible to evaluate one or other of the integrals, an 
explicit solution of the equation is impossible, and the 
solution by quadratures must be regarded as the best 
attainable, unless an alternative line of approach can be 
discovered. 


Example 1. 

x(y * - 1 )dx - y(x * - 1 )dy = 0. 
Separating the variables: 

xdx ydy _ 
x» - 1 y 3 - 1 

log | x* - 1 | - log I y 1 - 1 I =0 
x* - 1 


Integrating, 

or 


2 /* ~ 1 


= - log o 


log 

which may be written 

(!/* - 1) =c(x* - 1). 

(Note that replacing the arbitrary constant O by another 
arbitrary form, as in this case G = - log c, may help to simplify 
the general integral.) 


Example 2. 


dx 


dy 


= 0 . 


V(i -*■) Vd-y') 

The variables are separate; direct integration gives 

arc sin x + arc sin y = 0 


u 
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which is the general integral. This may be transformed into 
an equivalent expression by slightly rearranging the terms 
and taking the sine of both members, thus: 

sin {arc sin y } =sin {C - arc sin x } 

or, since 

cos (arc sin x) = ± y/( 1 - x*), 
y = ± \/( 1 - x*) sin C - x cos 0 
and, rationalising, 

(y +xcos C)*=(l -x*) sin* 0, 
i.e. 

x* + y* + 2xy cos C =sin* G 

or, if c =cos C, 

x* + y* + 2 cry = 1 - c*. 

Example 3. A change in variables may sometimes succeed 
in converting an oquation into another with separate variables. 
For instance, in 

(x + y)dx + dy = 0 

the variables are not separable, but if y is replaced by v -x, the 
equation is transformed into 

(v - 1 )dx +dv = 0. 

The variables x and v are separable, thus 

, dv 

dx +-= 0 

v - 1 

which leads to 

x + log | v - 1 | = log c or (t? - 1 )c* =c 

so that the original oquation has a general integral of the form 

(x- l)c* =c, or y =ce~ x -x + 1. 

I ho most important instance of reduction by change of 
\ uriahles occurs in the cuse of the homogeneous equation which 
now follows. 

4. The Homogeneous Type. The equation 

y)dx + Q(z, y)dy= 0 

is said to be of homogeneous type when P and Q are 
homogeneous functions of x and y of the same degree. If 
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the degree is m, the substitution y = vx will reduce P and Q 
to the forms 

P(x, vx) = x rn R{v), Q{x, vx) = x m S(v), 

where R and S are independent of x. Thus the factor x m 
may be cancelled out of the equation, which becomes 

R(v)dx + S{v){vdx + xdv) = 0 


or 


{i?(y) + vS(v)}dx + xS(v)dv = 0. 
Separating the variables and integrating, we have 

S(v)dv 


I 


+ log X = C, 


R(v) + vS(v) 

and when the integral in v has been evaluated, the sub¬ 
stitution v =y/x will give the general integral of the original 
equation. 

Example 1. 

(z* -y')dx +2xydy =0. 

The two terms in this equation aro homogeneous and of the 
second degree in z and y; the above process is therefore 
applicable. Making the substitution mentioned, wo have 

z*(l - v*)dx +2x t v[vdx +xdv) =0; 

z* cancels out, leaving 

(1 +v')dx +2vxdv =0. _ 

Separating variables uiul integrating: 


or 


t.e. 


[“,[-=« 

J 1 + u* J* 

log (1 + u*) + log z = log c. 


(1 + v*)x =c 

which loads to the general integral 



x* +y* =cx 


Example 2. 


(2 ye v,x -x)y' + 2z +y =0. 



8 


DIFFERENTIAL EQUATIONS 


Here each term is of the first degree in x and y, for e v/ * is of 
degree zero. Writing y=vx, y' =xv' +v and cancelling x, we 
obtain 


(2ve v - l)xv' +2(v*e v + 1) 

Separating the variables, this becomes 

2v-e~ v , 2c lx n 

dv + -=0. 


=0. 


v J +e 




Integrating, we have 

log (v* +e _, ’) + 2 log x=0, 
whence the general integral 

t/* +x*e~ v ^ x =c. 


5. The Equation with Linear Coefficients. Although 

the equation 

(ax + by + c)dx + (a’x+ b’tj +c')dy =0 . (5.1) 

is not of homogeneous type, it may be reduced to that 
type by a simple substitution. The equations 

ax + by +c=0, a'x +b'y +c' =0 . (5.2) 

represent a pair of straight lines which will intersect unless 
the condition for parallelism, i.e. a'/a =b'jb or ab' -a'b=0, 
is satisfied. Let (/*,£) be the point of intersection; transfer 
the origin to that point by the substitution 

x = h +X, y = k+Y 

nnd the equation will become 


(aX + b Y)dX + {a'X + b' Y)d Y = 0. 


It is 
followed 
integral 


now homogeneous; the substitution Y= vX 
by separation of variables leads to the general 


log CX + 


fc 


(a' +b'v)dv 


+ (a' +b)v + b'v 2 



whose ultimate form * depends upon whether the roots 


• For which soe § 17 infra. 
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of the denominator of the integrand are real, coincident 
or imaginary, t.e. according as (a'+6) 2 is greater than, 
equal to, or less than 4 ab'. 

In the exceptional case when the lines (5.2) are parallel, 
that is when a'/a=b'lb = k (say), the equation can be 
written 

(ax +by+ c)dx + {Jc(ax + by) + c')dy =0. 

Take z = ax +by as a new variable to replace y, then 

b(z + c)dx + (hz + c')(dz - adx) =0. 

Separating the variables and integrating, we have 


x + 


f_^ 

J(6 - a 


(kz +c')dz 


(b -ak)z +bc -ac' 


= const. 


The above are particular cases of an equation of the 
typo 

l ax + by + _c_ 
y V a'x + b'y + c' 


which may be reduced to a form integrable by quadratures 
by the same routine process. 


Example 1. 



4x - y +7 
2x+y - 1 


The lines 4x -y +7 =0, 2x +y - 1 =0 meet in ( - 1, 3); writing 
x—X — \,y — Y + 3 we have 

(2X + Y)dY =(4X -Y)dX. 

The transformation Y =vX reduces this equation to 

(2 +v)Xdv +(t>* +3y -4 )dX =0 

which becomes, on separating tho variables and tuking partial 
fractions, 


/ 3 2 I, 5dX 

S-- +-- \dv + —=r =0. 

{v - 1 v + 4j X 


Integrating. 

3 log | v - 1 | + 2 log | v + 4 | + 6 log | A' | = C, 
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(v - l) 3 (t> + 4)*.Y 8 =c 
or 

(7-X) J (F + 4.Y) l =c. 

Reverting to the variables x, y we have the general Integral 


(t/-x-4) s (i/ + 4x + l)*=c. 


Example 2. 

(2x - 4y + 5)y' + x - 2y +3 = 0. 

If 2 = x - 2y, 2;/' = 1 - 2 '; with this transformation the equation 
becomes 

(22+5)2'=4z+ 11. 

Separating the variables, we obtain 

(‘-5Tu)* =M j! 

whence 

42 - log | 42 + 11 | = 8x - C, 
giving the general integral 

4x + 8y + log | 4x-8i/ + ll |= C. 


6. Exact Equations. When the primitive of a differential 
equation involves the arbitrary constant C explicitly, as : 

u (z, U) = C, . . . ( 6 . 1 ) 

the operation of taking the differential eliminates C 

automatically, thus: 

* * 

du(x,y)= 0 . . . (6.2) 

or 


du 

Conversely, if a differential equation of the form 



P(x, y)dx + Q(x, y)dy=0 . . (6.4) 

has originated in such a process, and if no variable factor 
h.ii been cancelled out, it must be equivalent to one of 
the form (6.3) and thus to (6.2), and therefore it must 
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possess a general integral of the form (0.1). Such an 
equation is said to be exact. 

Thus, in order that (6.4) may be exact, there must 
exist a function u(x, y) such that 


du _ da 


(6.5) 


and therefore _ 

8P _ d 2 u _8Q' 

dy dxdy dx 

Hence the theorem: for the differential equation (6.4) 
to be exact it in necessary that P(x, y) and Q(x, y) be linked 
by the relation 

dP dQ 

dy dx’ 

This is known as the condition of integrability* for when 
it is satisfied the primitive u(x, y)=C can be recovered by 
the following process. Starting from the relation 

du(x, y) 


( 6 . 6 ) 


dx 


P{z, y), 


integrating and remembering that since the differentiation 
with respect to x was partial, the inverse process of 
integration will introduce, as the arbitrary element, a 
function of y, we have, 

u(x, y) =fP{x, y)dx+<f>(y) 

= S{x, y)+<f>{y) (say). 

But, on the other hand, 

du dS % 

an equation which will give <f>'{y) since Q and S are both 
known; the final integration to obtain <f>{y) will introduce 
the arbitrary constant C of the general integral. 


• Immediate integrability is implied; when the condition is not 
satisfied, the equation is still integrablo, though not without some 
preliminary manipulation. 
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Example 1. 

(1 +y i )ydx + (1 +x')xdy 
(1 +z*+y*) 8/ » 

Since, in the above notation, 

3F_ 9 f y+y * 1 _ 1 +z* +y* +3x*i/* 

% ”0*1(1“ ' (1 + x> +y*) iTi ~ 

_ dj x + z* | 
“M(l +x*+y s ) 3/ *J = Hx 

the equation is exact. Hence we are entitled to write 


■.(*,»)=i +m 

* '(l+x'+y*)*'* 


=1 

f yd* f 

■J(l +x* +t/*) ,/ * J(1 +x i +y')*'' + ^ y) 

■I(T7^)«. + I*4{(I ^ ? + y)4 fc + m 


x'ydx 


V(1 +y ») 

on integrating the second integral by parts. The equation 

&u _ x +x 3 Jf/ 

By-(i+x'+yif* + * (y) 

shows that <}>'(y) is zoro, or <f>(y) is a constant. The general 
integral is therefore 


\/(l + ** +y *) 

Note that although the given equation is exact as it stands, 
it would ccaso to be exact if the denominator of the left-hand 
member were removed. 


Example 2. 

log (y 1 + 1 )dx + —~ = 0. 

!✓* +1 

1 ho condition for integrability is satisfied, therefore we 
write 

«(*. y) = Jlog (y* + l)dx + f(y) 

= xlog (</* + !) + <f>[y). 
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From this we obtain 


du 2 xy 


+ *'(y); 


dy y' +1 

comparison with the coefficient of dy in the given equation 
shows that 

or = - i°g (y* + 1 ) + con3fc - 

The general integral therefore is 

(x - 1) log («/* + 1) = C. 


7. Integrating Factors. It will now be supposed that 
the condition for integrability (0.6) is not satisfied. This 
being so, an explicit general integral u(x,y)=C is not 
immediately obtainable, but, on the other hand, the 
theory of the existence of integrals is able to establish 
the fact that a general integral in which the constant C 
is implicit does exist. This may be interpreted to mean 
that, given any non-exact equation 

Pdx + Qdy =0, . . • (7.1) 

there always exists an integrating factor p(x, y) such that 
the modified equation 

y.Pdx + fxQdy = 0 . . • (7.2) 

is exact. Assuming that this is so, /x(x, y) must satisfy 
identically the relation 

s ±njwi . . . , 7 . 3 , 


or 



du j dP 8Q\ 

Q dx +,X {dy dxj 



(7.4) 


To find a general solution of this partial differential 
equation is a vastly more advanced problem than that 
of solving the ordinary differential equation originally 
proposed. For present purposes, however, any particular 
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solution of (7.4) would suffice, and such a solution may 
frequently be obtained by trial. When any integrating 
factor ix, satisfying (7.4), has been obtained and introduced 
into (7.2), this equation becomes exact and may be 
integrated by the method of § 6. The limitations of this 
method of integration are obvious; nevertheless there are 
particular cases which may be treated systematically; 
these will be indicated in § 9. 


Example 1. The equation 


is not exact. 


or 


ydx - xdy = 0 

Any integrating factor p will satisfy 

_ _ 2(xp) 
by dx 


It can be verified that 


d/i dp 

X 8x +y d y +2 ' l = 0 - 


I 1 1 

.1 

x y *y 

are among the possible values of p. Taking p = 1/a:* we have 



which is exact and has the solution y/x = C. Similarly p = 1/y* 
leads to the equivalent solution x/y=C'. Taking p = 1 jxy we 
sopurate the variables: 

dx dy 

- -=0 

x y 

or 


log | x | - log | y | = const., whence xjy = const. 


Example 2. 

(1 +x* +?/>)-* » i 9 an integrating factor for 

(1 +y%rfx + (l +x*)xdy=0. (§0. Ex. 1.) 

Another integrating factor is tho reciprocal of xy(l +x ! )(l +y*); 
this separates the variables. 
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8. The Quotient of Two Integrating Factors. Let n x 
and fj. 2 be two integrating factors of 

Pdx + Qdy = 0, . . . (8.1) 

whose ratio is not of itself a constant. Then the equation 

Pilpi ~ C or fj. 2 = C[x x • • (8.2) 

is an integral of (8.1). For (8.2) is a primitive of the 
differential equation 

But since [x x and fx 2 are integrating factors of (8.1) 




and hence 


\ dy 

P\h. . & 

Uy 


dx 


+ /x 2 


-g}-- 

-si-- 


(8.4) 


(8.5) 


M p Tu~ Q % 


t.e. 


.{ 

{*%'-5 }'•-{'■ 




dfi x d/M 
dx /Xl dx 


} Q ■ 


which reduces (8.3) to (8.1), thereby proving the theorem. 

This theorem is equivalent to the statement that il 
one integrating factor of (8.1) is known, an unlimited 
number of others may be obtained. Thus if fx x is known 
and we write fx 2 = v Pv we have from (8.5) 




dy " dx 
which, on account of (8.4), reduces to 


[WJQ\ = 

\dy dx] 


0 


Pi 


\ dy H dx\ 


so that v = const, is any general integral of (8.1). But if 
u =<? is one form of this general integral, we have v = F(u). 
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Thus 


th=Pi F (u), 

where F(u) is any arbitrary function of u. 


Example. Referring to § 7, Ex. 1, we see that solutions of 
the equation are 

11 11 11 

* 1 = const., — -r — = const., — -f — = const., 

x y J x 1 xy ,f xy 

all of which are equivalent to y/x=C. Furthermore, an 
infinite number of integrating factors is included in the formula 



9. Special Types oi Integrating Factor. It may 
happen that an integrating factor can be found which 
depends on one variable only. Consider, for example, the 
circumstances in which 

Pdx + Qdy= 0 

admits of an integrating factor p(x) depending upon x 
alone. That being the case, (7.4) will become 



11ms p can be a function of x alone if the right-hand 
member of this equation is independent of y; p is then 
determined by a quadrature. 

Example. 

(1 -xy)dx+(xy -x i )dy=Q. 

In the above notation, 

(8P dQ\f -x-(y-2x) 1 
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Hence 



log | fi |= -log|x| or /i 


H x 

With this factor introduced, the equation becomes 


1 

~ • 
x 


- yjdx + (y - x)dy = 0; 

it is now exact and has the integral 

log | z | - xy + iy*=C. 

In the same way, assuming the existence of an in¬ 
tegrating factor which is a function of x + y, say 

H=f(x+y) =/(z), 

it will he found that 


m 

/(*> 


-ffi-S}/“•-« ■ 


. (9.2) 


Therefore a necessary condition for the existence of such 
an integrating factor is that the right-hand member of 
(9.2) shall be a function of z, or x+y, alone; fi is then 
determined by a quadrature. 


Example. 

(6z* + 2xy + 2y*)dx + 3(x* +xy l + 2y*)dy=0. 

Condition (9.2) hero becomes 

f\z) Gy* - 4x 2 _ 2 

m ~ ~2x*+2 xy - 3zi/* -3 y* x+y z 

and therefore the condition is satisfied. Wo have/(z) =z* and 
therefore the integrating factor is (z + y) x . Introducing it and 
integrating, wo obtain the general integral 

(x* +r/ 5 )(x +i/)* =G. 

Again, if it is assumed that 

H- =/(*»/) =/( 2 ) 
the condition will become 


m 

/(*) 



(Px-Qy),. 


(9.3) 
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where the right-hand member must be a function of t, 
or xy, alone; here again p is determined by a quadrature. 

Example. 

(xy* + 2 x*y* - y*)dx + ( x'y * + 2 x*y - 2 x*)dy =0. 

Condition (9.3) becomes 

JVl xy* - 2x'y -2y+Ax = j _ £ = j _ ? 
f(z) xy'-2x*y xy z 

whence /(z)=c*z _ *. An integrating factor is therefore 
// =e zv x~ t y~ i i it leads to the general integral 



10. The Linear Equation. An equation of the type 

f(x)y' +g(x)y = h{x), 

which involves y and y' linearly, is said to be linear. By 
dividing throughout by f(x) it may be brought into the 
standard form 

y'+py=q, . . . (10-1) 

where p and q are functions of x alone. When condition 
(9.1) is applied to the equivalent dillercntial form 

dy + {py - q)dx =0 

it reveals the existence of an integrating factor which is 
a function of x alone. In fact, if the factor p(x) is intro¬ 
duced the condition of integrability reduces to 


or 


dp V(ppy - pq) 

dx dy ^ 


dp/p = pdx. 

Integrating, 

log p =jpdx, whence 
The exact equation 




. (10.2) 


e\ fdl dy + (py - q)e^ l ' Jl dx = 0 
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(10.3) 


( 10 . 4 ) 


has the general integral 

ye*** 1 -fq^dx-C, . 

which thus involves two quadratures. 

No/e 1. When the integrating factor , haa been obtained 
and introduced, the equation may be written 

fidy + n'ydx =qpdx (for pp = P ) 

° r d(py) =qydx, 

and integrating, 

py=) qpdx+L" 

The general integral is therefore of the form 

y = Cp~ x +p- i !g/ ldx - • 

Note 2. Replacing q by zero in ( 10.1) we th ® 
corresponding reduced equation, homogeneous m y y 

y'+py= 0; • • • < 10 - 6) 

its general integral involves only one quadrature: 

y=Ce~^. 

No/e 3. If any particular solution of (lOdl^ay^n^ 
known, the general integral is obtaina y 4 
For we then have 

yi'+py i =< 7’ 

and eliminating , between this equation and (10.1); also 

writing v for y - y» wo liave 

v' + pv =0, 

t.e. the corresponding reducofl equation. Thus the general 
integral of (10.1) ia 

y =v +y\ 

= Ce-l ,ldz +Vi- • * ' (1 ° -6) 

, * „ 4 Via first is takon from the general 

It consists of two terms, tl known as the com- 

integral of the reduced equat.on and ,» know/^ 

■plernentary function, the second is any I 

w . . The difference between any two particular 
integral is a social case of the complementary function. 
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For since the relation 

y-y,=C<rl’ d ’ 

is general, any particular integral y t is obtained by attributing 
the appropriate value to C, say C': 

Note 5. If y l and y, are any two particular integrals, 

y - 2/1 


= const. 


. (10.7) 


2/a " Vi 

furnishes the general integral. 

Example. 

xy’ -(x + l)y =x* -x s . 

Dividing by the coefficient of y' we reduce the equation to the 
standard form 

, x + ! 

y - y=x~x\ 

X 

If ft is the integrating factor, 

fx + 1 


i •*•+ 1 , 

log /»=- dx = -x -log x 
* 


or 


y = e~ x /x. 

Introducing /x into the equation in its standard (not its original) 
form, wo have 

e z y’ x + 1 

e-*y=e-*(l-x) 


or 


d (e~ x y\ 


i.e. 


e~ x y/x=C +fe- x { 1 -x)dx =C +xe~ z , 
so that the general integral is 

y=Cxe x +x*. 

Application. A condenser of capacity C (farads) is being 
charged from a source of electricity of potential E (volts) 
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through a non-inductive ^tance * jSw 

(coulombs) at time t (seconds) is given by the linear 

equation dQ Q _ „ 

dt C 

If the initial charge is aero, it is »q»tadto And Q at any 
instant (i) when the voltage ® - —t 00 ^ ? ^ 
voltage is alternating, t.c. & o P 

constants. . *_„_otion we find that. 

Following the general method of mtegration 

whatever E may be, 


QjlRC = K + ^ i e , l nc dt, 


where K is the constant of integration. 

(i) When E is constant, 

q<JIM = k + CEc‘l RC . 

The constant K is determined by the initial condition that 
Q =0 when t =0, so that 

a _ if J.C.R. 


Q t t\na = CE(e , ' nc - !)• 

Thus the charge at time t is given by 

Q =CE{ 1 - e~ , l nc )t 

which approximates to the steady value CE when < is large 

(ii) When E= Exempt, 

e*l nc Q =K + ^je‘ //iC sin pi dt 

= K+CE 0 Jl nc (*in pi -RCp cos pO/O +«* C V)- 

Since Q =0 when e =0 we have 

0=K-RC'pEJ{l+I< i C t P t ) 


CE, 


and hence 

Q = .. 

The influence of tho exponential term diminishes r 


-(sin pt-RCp cos pi +RCpe- t,R0 )- 


1 + R x C'p 

The influence of tho exponent.a. ~ a ~ 0 " pC riod E. 
that Q approximates to a form of t c 
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11. The Bernoulli Equation. The equation 

y'+py=qy n . . . (li.i) 

in which p and q are functions of x alone, is associated 
with the name of James Bernoulli. It may be reduced 
to linear form by a change in the dependent variable. 
For the result of division throughout by y n : 

y~ n y' +py 1 ~ n =q 

suggests the substitution 


v=y l ~ n , i/ = (l -n)y~ n y'. 
The equation then becomes 

t/+(l -n)j>y = (l -n)q, 
which is of the standard linear type in v. 
Example. 

dy y r 1 5 

Tx~2x = ° Xy ' 

The above steps give, in turn : 

y‘ _ 1 

y b 2 xy l 


= 5x* 


v=y~\ v' = -ly-*y' 


. (H.2) 


v' + ~v = - 20.c ! . 

X 

This linear equation lias the integrating factor x*i 

x i o' +2rr= -20 r*. 

Integrating: 

x*v=C-4x i or v = Cx -* - 4.c* 
and since y = v~ Vi , 

1 

y ~ - 4 ^)* 


12. The Riccati Equation. The standard form of this 
equation is 


y' =py-+qy + r, . 


• (12.1) 
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„b«....I—"-»; * ■‘35 “"i’ 

by the substitution + \/v 

where . is a new dependel^nl iM»g» 
it will be found that the equation in a is 

v' +{2p>ji +g) y + v =0 • ' 1 

and the integration can bo completed by quadratures. 

Note 1. The general integral may be the 

in terms of any three tunc ‘‘ OD ^" 8 *b.«tution »=». + >/»■ 
equation. For if we ma par ti C ular values of v. say 

y , and y 3 will correspond to two par 

v x and v t$ thus: .. 

y t «* + !/*» ys=yi + l l v " 

, Qro particular solutions of the linear 

and since w = «h and u * P' . iU bo (io.7) 

equation (12.2) its general solution 

!^=C. 

«>, - Vi 


Making the reverse substitutions 

1 1 


v = 


we find that 


y -!/» 


V, = 


yt-yi 


V. = 


Vi -y i 


!/-!/*_ ~y * 


y - 2/» 2/3 - !/i 

which furnishes the general integral. 

No<e2. !£,-**■ a fourth particular solution, wo have 

y«-y« y» -r», 
y* - yi -'/a - y« 

the cross-ratio of any four solutions is a constant. 
Example. *». 
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The solutions y=r, y = -x are evident on inspection. Taking 
the former, we write 

y=x + l/v, y' = 1 - v’/v* 

and derive the linear equation 

2x*(i>' +v) = l -x; 

its general integral is 

v = (Cxe~ x - l)/2x. 

The general integral of the proposed equation is therefore 

2x Cxe~ x +1 Cx+e* 

y X Cxe~ x - 1 ~ X Cxe~ x - 1 - *Cx - e x ‘ 

The cross-ratio of the four particular functions x, -x, 
x+e x x - c* 

x x _ eX > X ^~ +C ~ x satisfy the equation will be found to 

be - 1. 

13. Change of Variable. To sum up, equations which 

may he readily integrated by quadratures are those with 

separable variables, those of the homogeneous type, exact 

and linear equations. In the homogeneous type, the 

variables are made separable by the substitution v = y/x 

by which a change of dependent variable is effected. 

Actually the method of integrating an equation with 

separable variables, or a linear equation, is to introduce 

an integrating factor which renders it exact; this is the 

fundamental principle upon which integration of equations 

of the first order and degree essentially depends. If 

other equations can he integrated it is because a change 

in variable (as in the case of the Bernoulli equation) has 

enabled a reduction to one of the above types to be 
performed. 

When an equation is proposed for solution and does not 
appear to come under one or other of the foregoing headings, 
it is advisable to consider if a change in one or both of 
the variables cannot be discovered which will transform 
t *e equation into one of a recognisable type, just aa a 



EQUATIONS OF THE FIRST ORDER AND DEGREE 


25 


change in the dependent variable of a Bernoulli equation 

rendered it linear. f 

For instance, if an equation is of the torin 

=f(ax +by+c ), 
dx 

the substitution v =ax + by at once suggests itself. It 
transforms the equation into 

< ^ = a + 6/(t>+c), 
dx 

where the variables are separable. 

Again, an equation such as 

{*/(y)+9(y)}-£ -& (y) 

may not be integrate as it stands but if 
independent and * as the dependent vanable, it may 

written 

dx 

k (y)j- y -f(y) x =o(y)> 

when it becomes a linear equation in x. 

It is not possible, even if it were desirable to fotmulate 

special mention, thus n naturally suggests 

ran p g rrpo.r:fo"s ^ *. 

x=r cos 9, y ~ r s ' n 
xdx + ydy - rdr, *dy - ydx -rMB. 

In this case it is usually best to leave the general integral 
expressed in terms of r, 9. 
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Example 1. 

(xy + 1 )dx + 2x 5 (2 xy - 1 )dy = 0. 

The prominence of ry in the equation suggests that v-xy 
might be token as a new variable in place of y. The equation 
then becomes 

(4t’ s - 3c - 1 )dx - 2x(2v -\)dv=0. 

The variables are separable; the usual process of integration 
leads to the general integral 

(4x«, + l) 3 i*«/-l) s =Cx\ 

Example 2. An equation of the type 

xy' - y =f(x)g{ylx) 

is transformed by the substitution y=vx into 

*V =f(x)g(v), 

whose variables ure separable. For instance, 

V 1 - .I s 

xy' - y - 2x-- 

J J x 4 - 1 

becomes, after transformation and separation, 

dv 2xdx 
r* - 1 x 4 - 1 

whence 

U — 1 =C— 

t> + 1 x ! + 1 

If c=(l - C)j( 1 +<?), the general integral is 

x s +e 
y=x -~— 
cx s + 1 

Example 3. 

(x s -f ir)(x<lx + ydy) + (x 1 + y* - 2x + 2 y)(ydx - xdy) = 0. 
Changing to polar co-ordinates, we obtain, cancelling r J , 

dr - (r - 2 cos 0 + 2 sin 0)d0 =0. 

I his is linear in r and has t ho general integral 

r=(V - 2 sin 0. 



CHAPTER II 

INTEGRAL CURVES 


14. Families of Plane Curves When referr 11 to ect 
angular axes of co-ordinates, the equation /(*, j)-0 
represents, for any definite value of the ^ 

plane curve. The aggregate of all values of C for wh cl 
The equation is possible gives a complete on^paramefer 

family of curves. Thus the equation x* + y = C ' 

all positive values are in turn assigned to C ^Presents 

the family of aU circles whose centres are at the origin. 

The corresponding differential equation 

P(x, y)dx + Q(x, y)dy=0 • • ( 14<1 > 

", ™ ^ “5 

:r■= 

of (14 11 The only exception to this statei 
when the ratio P/Q becomes indeterminate, a case 

will be investigated later. nation /14 \) 

Thus the integration of the differential equation (14. ) 
is equivalent to the geometrical problem of discover! g 
the Sy of plane curves, A. inUyral 
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This amounts to the statement that through any point 
(x, y) in general there passes but one integral curve of a 
differential equation of the first order and degree. 

The symbolic expression of any geometrical problem 
that involves, directly or indirectly, the gradient at any 
point is in fact a differential equation of the first order. 
It may, however, be of higher degree than the first, but 
for the moment problems that lead to equations of the 
first order and degree will suffice as examples. 

If P be any point (x, y) on a plane curve, let T, N, 0 
be the points at which the tangent, ordinate and normal 
at P respectively meet the x-axis. Then if ip is the angle 
between TP and the positive direction of the x-axis, 


y' =tan ifi. 

The gradient may be introduced indirectly, e.g. 


Subnormal = NG = NP tan i{j =yy' = 

dx 

Subtangent =TN =NP cot =y/y' =y—. 


The lengths PT, PG are known as the “tangent” and 
“normal ” respectively; they involve y' but not linearly. 


Example 1. Curves whose subnormal has the constant 
value 2a: 


Integrating: 


yy' = 2a or ydy = 2 adx. 
y 2 =4 a(x + C) 


which represents a family of equal parabolas, with axes 
coincident with the x-axis. 


Example 2. Curves whose subtangent has the constant 
value a: 


dx du dx 

y =a or - - - 
dy y a 
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Integrating: 


logt/=logC+- or y-Ce xa , 


a family of exponential curves. 


Example 3. Curves in which the subtangent is twice the 
abscissa : 

y iX = 2X. 

y iy 

The variables are separable: the integral curves are the family 
of parabolas y * =Cx. 

Example 4. The ordinate at any point is equal to the 
perpendicular from the origin to the normal at that point. 

If (X , Y) is any point on the normal at (x, y )• 

X + Yy' =x+yy' 

and therefore the length of the perpendicular from the origin is 

.*+ygL =y . 

vti+y' 1 ) 

This symbolic statement involves y'» as it stands,, but on 
rationalisation and simplification it reduces to an equation 

the first degree in y': 

x*+2 xyy = y l 
2xydy + (x* - y* )<** = °- 

This can be integrated as an equation of homogeneous type 
(S 4 Ex 1) or as one which has for an integrating ac 
function of x alone. The second method is the simpler; the 
integrating factor *-• leads to the general integral 

x* + y* = Cx, 

representing a family of cirelee touching the y-axis at the 
origin. 

15. Trajectories. Let 

f(x,y,a) = 0 and g(x,y,P)=0 
be the equations of two families of curves each dependent 
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upon one parameter. When each member of the family j3 
cuts every member of the family a according to a definite 
law, any curve of either of the families is said to be a 
trajectory of the other family. The most important case 
is that in which curves of the families intersect at a 
constant angle. 

When a curve intersects all the members of a family 
of curves at right angles, it is said to be an orthogonal 
trajectory of that family. Thus any circle of centre 0 
is the orthogonal trajectory of the family of straight lines 
radiating from 0, and conversely any line through 0 is 
the orthogonal trajectory of all circles concentric with 0. 
That is, the families 

x ~ + >r = a 2 , y = fa 

are mutually orthogonal. 

In rectangular co-ordinates let 

P{x, y)dx + Q(x, y)dy =0 . • (15.1) 

be the differential equation of one family of curves; when 
written in the form 

y' = - P[*, y)IQ{x, y), 

it expresses the fact that if an integral curve passes through 
(.r, y) it will have a certain gradient at that point. The 
gradient of the orthogonal curve through that point will 
be the negative reciprocal, i.e. 

y'=Q(x, y)/P(x, y). 

In other words, an orthogonal trajectory will be an integral 
curve of the differential equation 

P{x, y)dy-Q{x, y)dx= 0. . . (15.2) 

Moreover, whenever an integral curve of (15.2) meets an 
integral curve of (15.1) it will cut it at right angles. 

11ms the differential equation of the one family is 
obtained from that of the other by replacing dx by dy 
and dy by -dx, or, what is the same thing, by replacing 
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,J by -1 ly'. If the differential equation of the family 
f(x,y, a) =0 is F(x,y f y')= 0, the differential equation of 
the orthogonal family will be F(x, y, -l/j/')=0, and its 
general integral g(x, y, p) =0 will be the orthogonal family 

itself. . 

When the curves of the families f(x, y, a)=0 and 
g(x, y, (3) = 0 intersect at a constant angle other than a 
right angle, the trajectories are said to be oblique, lo 
make the problem definite, let the tangent at (*, y) to 
the curve of / = 0 make the angle 0 with the a:-axis; let 
the inclination of the tangent to ? = 0 at (x, y) be 0 t and 
let the law connecting them be 0 X = 0 + ^ cn 

tan i/»! — tan cu t an 0t - m 
tan ^ = 1 +tan 0 X tan to ~ 1 + m tan 0/ 

writing m for tan co. . 

Let the differential equation of the original family he 

F{x, y, y') =0, 

or since y' =tan 0 in the original family, 

F(x, y, tan 0) =0, 

or again 

/ tan 0i - m \ 

F \ C} V ' 1 + m tan 0 X ) 

But in the family of oblique trajectories y’ = tan 0, and 
therefore this family consists of the integral curves of the 
diflcrential equation 

If therefore the differential equation of a family of 
plane curves is known, the differential equation ot the 
family that intersects it at a positive angle cu is obtained 
by replacing y' wherever it occurs by (y - tn) /( + m y)* 
where w=tan oj, or, what is the same thing, by rep acing 
dy by dy - mdx and also dx by dx + mdy. 
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Let f(r, 6, a) represent a one-parameter family of 
plane curves in polar co-ordinates, and let the differential 
equation of the family be F[r, 9, dr/d9) = 0. Let (r, 6) 
be any point through which a curve passes, and let <f> be 
the angle which the tangent at that point makes with the 
radius vector, so that 



An oblique trajectory, of angle co, will be such that the 
angle which its tangent makes with the radius vector will 
be </>,=</>+ co, and we have 

tan <f>= tan (<f> l - co) 


tan <f) 1 - tan co rdd/dr - m rdO - mdr 
1 +tan (f> l tan co 1 + mrd9/dr dr + mrdO' 


Thus the differential equation of the oblique trajectories 
will be obtained from that of the original family by 
replacing rdO by rdd - mdr and dr by dr + mrdO. 

in the case of orthogonal trajectories we should have 

<f>i =<f>± 



1 dr 
r dd 


bo that rdO is replaced by -dr and dr by rdd. 


Example 1. The equation 

x 1 + f/ 1 - 2ax = 0 

represents, for varying a, all circles which touch the j/-axis 
ut the origin, the centre of any individual circle being (a, 0). 
Its differential equation (§ 14, Ex. 4) is 

(x* - y*)dx + 2rydy =0. 

lo find the equation of the orthogonal trajectories replace dx 
by dy, dy by -dx: 


(x* -y*)dy - 2xydx =0, 

which is really the same equation with x and y interchanged 
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throughout. Its general integral therefore is 

x* + i/*-2/?y=0, 

so that the orthogonal trajectories are all circles which touch 
the x-axis at the origin. 


Example 2. It may happen that any one member of a 
family intorsects at right angles an infinite set (not necessarily 
all) of members of the same set. Such a family is said to be 
self-orthogonal* For example, consider the family of confocal 

parabolas 

y' =4 c(x +c), 


t.e. parabolas with a common focus at the origin and common 
axis Ox. 

We have _ n. 

y dx 


and eliminating e. 




or 


y(dy' - dx *) + 2xdxdy = 0. 

This, the differential equation of the family of parabolas, 
remains essentially unaltered if dy is replaced by -dx and dx 
by dy, so that the family is, in fact, self-orthogonal, ro see 
how this property originates, consider the two particular 


curves 

y* = Aa(x + a), y % = 4b ( x + b b 

Their point of intersection has the abscissa * = - (a +6) and 
the ordinate is given by y* = - 4a6. Thus the intersect on is 
real only when a and b are of opposite signs. When that is the 
case, the gradients of the two curves at the point of inter- 
section are respectively either \/( -a/b) and v( / 

- V( - a/6) and V( -6/«). which verifio3 th ° ,r <> rtho 8 onill,t y- 

Example 3. To find the oblique trajectories of the family 
of straight linos 

v = ax. 


* This cannot happen 
is of the first degree in x/ \ 
the second degree. 


if the differential equation of the family 
the example given leads to an equation of 
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The differential equation of the family is xy' =y or 

xdy =ydx, 

and therefore that of the trajectories: 

x(dy - mdx) =y(dx +mdy) 
or 

xdy - ydx = m(xdx + ydy). 

Clearly this is a case for a change to polar co-ordinates, which 
gives 

rdO =mdr. 


The general integral is r = fle 9lm , 
family of equiangular spirals. 


i.e. the trajectories consist of a 


Example 4. In polar co-ordinates the family of circles given 
in Ex. 1 is 

r = 2a cos 6, 

and its differential equation is 

dr = - r tan 0 dO. 

Ihe orthogonal trajectories, which are the integral curves of 

rdO =tan 0 dr, 

are the family 

r = 2/? sin 0. 


16. Level Lines and Lines of Slope on a Surface. Let 

the surface be referred to a rectangular system of axes 
with the 2 -axis vertical; the ( x , y)-plane will be the hori¬ 
zontal reference plane. The section made on the surface 
by a horizontal plane z = a is a level line or contour line 
of level a. 1 he horizontal projections of the various level 
lines made by a regular sequence of horizontal planes z = a lt 
a 2 > • • • form a contour map. If a is varied continuously, 
as it does when the secant plane moves uniformly upwards, 
the contour map is replaced by a family of curves in the 
(•*, y)-plane dependent upon the parameter a. If the 
equation of the surface is f(x, y,z)=--0 these curves will 
be contained in the equation/(.r, y, a)=0. 

Ihe lines of steepest slope on the surface are the lines 
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of minimum distance between successive contoui-s. and 

are therefore lines cutting the level lines at right angle 

On account of the fact that when two curves in space, 

one of which lies in a horizontal plane cut orthogonaU^ 

their horizontal projections will cut orthogon y, 

that on the horizontal reference plane the 

the lines of steepest slope are the orthogcinal tra ecto es 

of the contour lines. If the equation of these 

is n(x V (3) =0, the lines of steepest slope on the surface 

is £A X » y> PI . surface bv the cylinders 

will be the intersections of the smtace ny y 

g(x> y> p )= o whose generators are parallel to the ax . 

Example. Consider the paraboloid 

2 z =ax* +by*. 

The level lines project into the family of curves 

ax 2 + by z = 2a, 

whose differential equation is 

axdx + bydy = 0. 

The orthogonal trajectories, whose differential equation 

axdij - bydx =0, 

will be found to be j/" = ^\ and hence the lines of steepest 

slope may be found CQse of a paraboloid of 

In particular, when a-6, t.e. i ^ intersections of 

revolution, the lines of steepest I they are 

the surface by all planes through the 2 -ax,s ; that is, thej 

meridian lines on the surface. 

17. Singular Points. _ Let (x 0 . !/o) b ea point such^ 
its neighbourhood, the functions ( , j), ’ 

continuous and single-valued. Then e 1 ‘ 

P(x, y)dx + Q(x, y)dy = 0 • , ' 

assigns the value y = - ^(* 0 - yo)IQ( x O' !fo) ^ ( 

at that point, and this value is determinate with the 

exception of the case (17 2) 
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Points at which the gradient becomes indeterminate 
through the simultaneous vanishing of P and Q are said 
to be singular. Thus singular points are the points of 
intersection of the curves P(x,y)= 0 and Q(x, y)=0, 
which are the loci of points at which the integral curves 
have tangents parallel to the x-axis and the y-axis respec¬ 
tively.* They are thus, in general, isolated points. 

To illustrate possible modes of behaviour of the integral 
curves in the neighbourhood of a singular point, we shall 
consider the linear fractional equation 

dy ax+by 

• < 17 ' 3 > 


for which the origin is the only finite singular point. 
Making the usual substitution y=vx, we obtain 

a + (b -l)v -mv 2 


dv 
x — 
dx 


l + tnv 


(17.4) 


If v - a, v - ft are factors of the numerator, v = a, v = ft 
will be particular solutions of (17.4), and y=ax , y = fix 
of (17.3) Thus among the integral curves of (17.3) are two 
straight lines through the origin, whoso joint equation is 

ax 2 + (6 - l)xy - my 2 =0. 

These are known as the principal lines, and a, j8 the 
principal directions through the origin. They are real 
and distinct, coincident, or imaginary according as 

(6 - 1) 2 + 4<m >, = or < 0. 

When they are real, (17.4) may be written 



m(v-a)(y- 1 3) 
/ + mv 



It ia nssumod that ^=0 and Q=0 have not a branch in 
common. If /^=0 were a branch common to both, R could be 
cancelled out of the equation, and the above statement would bo 
strictly true. 
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Separating the variables and integrating, we obtain 


log 1 x | + /x log | v - a | + v log | v - (3 | - log C 
or 

x[v - a)^(v - j3) v = C, 

where 

l+ma l+mfi 

v= mp^y M 

and thus the general integral of (17.3), when the principal 
directions are real and distinct, is 


(y - ax'fiy - fa) v = C. 

When [j. and v are both positive, this represents a 
family of curves having asymptotes y = ax, y = fix in 
common. The asymptotes themselves are integral curves, 
but apart from them, no integral curves pass through the 
origin, which is a col of the family.* When /z and v are 
of opposite signs, the equation of the integral curves may 

be written 

y - fix =c(y - ax)P, where p = - /x/v. 

Since p is positive, every integral curve passes through 
the origin, which is a base point or nodal point of the family. 
When p-1, i.e. ,jl= -v, the integral curves pass through 
the origin in all directions; the nodal point is then said 

to be isotropic. When p> 1, t.e. | p- | > I v l» y = aX an< 
y = Px are integral curves; all other integral curves touch 

the principal line y = at the origin. When 0 < p < > 
i.e. | p. | < | v |, y = ax and y = px are integral curves; all 

other integral curves touch y = ax at the origin. 

When the principal directions coincide with that ot 

the line y = ax, (17.4) may be written 


dv Tn(u-a) 2 

x j ~ + — i- 

dx l+mv 



The origin ia a aaddlo-point on the surfaces = (y -uxHy - pr)- 

1 > 
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whence 


kgM+ksl.-al-ij^-O, 


and thus the general integral of (17.3) is 

(y ~ ax)(log \ y-ax\-G) = ( l/m + a)x, 
from which we derive 

, l/m + a 

V a + 1 - C + log \ y - ax \ 

Thus the line y = ax, which is itself an integral curve, is 
tangential at the origin to every other integral curve. 

When the principal directions are imaginary, we may 
write a = K + i\, fi = K-i\, whereupon (17.4) becomes 

(lv 7/i{(v -k) 2 + A 2 } 


f— 4-'- 

dx l + mv 


= 0 . 


Separating the variables and integrating, we have 
log | x | + i log {(y - k ) 2 + A-} + - ^ arc tan —=log G, 
and hence the general integral of (17.3) is 

§ l°g {(!/ ~ xx) 2 + A 2 * 2 } + —arc tan -=l°g 0. 

If l + v ik =0, the integral curves 

(// - kx) 2 + X 2 x 2 = C 2 

are a family of ellipses encircling the singular point at 
the origin, which is their centre and limiting point. In the 
general case we make the substitution * 

y - kx = r sin 0, Xx = r cos 6 

and obtain 

or log r + p6 = log C, where p = (l + m#c)/mA, 

r =Ce~» e . 

* This transformation may ho regarded as a strain t]=y-icx, 
* ~ ' Jr consisting of u shear in the direction of the y-axis and an 
extension in the direction of the x-axis, followed by a change to 
polar co-ordinates. 
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This is the equation of an equiangular spiral which winds 
asymptotically round the origin. The transformation 
back to the x, y -plane involves no deformation other than 
a simple strain, so that the integral curves are in general 
spirals with a common asymptotic point at the origin. 

Note. When (17.1) has a singular point at the origin it 
may be written 

dy _ ax+by+p(x, y) 
dx ~~ lx+my+q(x, y)' 

where vlx v), < 7 (x, y) can be developed as series whose lowest 
l h Z fre of the s^ond degree at least It can be shown that 
unless o=6=0 or l=m = 0 the integral curves behave, m the 
neighbourhood of the origin, ns if p and , were absent, and 
precisely the same cases arise. 

Example 1. 

dy my (n ^ m , „ > 1). 


dx 


x 


, . . _„i „ _ ri r n _ x n ltn —w) shows that y is of the 

The general integral y=Ux x pn , vvi.^n 0 

order of x"* or x" according as m < or > ». When rn < 0, 
the origin is a col of the integral curves ; when m > 0, it is a 
noLl point, with integral curves touching the yax.s when 
m < 1, and tho x-axis when m > 1. 

Example 2. 

dy __ x +2x» 
dx y 

The integral curves „• +x> +*< = C> are closed ovals , the origin 

is their limiting point or centre. 

Example 3. , «. 

dy _ -x+ny -/*»/(**+S/ B ) 

dx~ y +px —/*x(x 2 + y 2 ) 

When transformed into polar co ordinates, this equation 
becomes *. =( . r(f . _ 1 )<10. 

and tho general integral is r*(l +C'”*)-1. The integral 
curves ore spirals with the origin as asymptotic point. 
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18. The General Integral. When the primitive is of the 
form f(x,y,C)= 0, involving the arbitrary constant C 
implicitly , the differential equation will arise through 
elimination of C between the equations 


/ = o, 



Only in exceptional cases will the result of this elimination 
involve y' to the first degree; if the resulting equation 

t/, r/') =0 . . . (18.1) 


is a polynomial in y\ involving y' to the with power, it 
is said to be of degree m. But it may be irrational in y' or 
even transcendental. 

Theoretically, (18.1) may be solved for y', giving a set 
of equations 


!/' = F i( x > y), y' = F t (x t y) .(18.2) 


each of the form hitherto considered. This set may be 
finite or, in the case of a transcendental equation, infinite 
in number. If the general integrals of the equations 
(18.2) are respectively 


/i(*,y.C)-0, Mx,y,C)= 0 , ... . (18.3) 

the general integral of (18.1) is defined to be any equation 
/(*, y, C) =0 which is satisfied when, and only when, at 
least one equation (18.3) is satisfied. In particular, when 
the equation is of dogree m, its general integral will be the 

40 
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product of the m integrals (18.3), that is 

C)Ux, y.G) . ■ O)-0. 

■» 'p* •« ""T' 

the present chapter. 

may be written 

y =mx 

where / is some specified function. Then 

p =y‘ = m - 

Eliminating m, we obtain the diSerontiel equation 

y =px +/(P)- 

In particular, the equation of the family y=m* + a/m, where 
a is a constant, is 

t/ =px +a/p 

or A 

- py + ®- u » 

of the second degree. 

Example 2. p , _ 2px +x , _ =0. 

This equation of the second degree may be decomposed into 

p -y-x= 0, p +y-x= 0. 

These are linear oquatione whose general integrals may bo 
WriU ° n y - Cc x + x + 1 = 0, y - Ce~* - x + 1 = 0. 
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Multiplying, we obtain the general integral of the given 
equation 

<7* - G{e x (y -x + 1) + e~ x [y +x +1)} + {y + 1)* -x* = 0. 

19. The Clairaut Equation. We consider the equation 
obtained in the last section (Ex. 1) 

y=px+f(p), . . . (19.1) 

known as the Clairaut equation, with a view to discovering 
whether or not it admits of any integral other than the 
primitive from which it was derived. Differentiating 
with respect to x, we obtain 

P=P+{x +f'{p)}j^ 

This equation may be satisfied in just two ways. On 
the one hand, we may take dp/dx = 0, whence p=C, and 
thus recover the primitive 

y = Cx +f{C). . . . (19.2) 

On the other hand, the equation is satisfied if 

and if x is now eliminated from (19.1), 

y = -pf'(j>)+f{p)- • • • ( 19 - 4 ) 

Equations (19.3) and (19.4) taken together are the para¬ 
metric equations of an integral curve of (19.1); by 
eliminating p we obtain the equation of this curve in a 
form cf)(x, y)=0 involving no arbitrary constant. This 
integral curve represents a singular solution (§ 2) of the 
equation. Since the envelope of the primitive family 
(19.2) is obtained by eliminating C between that equation 
and x +f'{C) =0, and since the result of this elimination 
is identical with that of eliminating p between (19.3) and 
(19.4), it follows that the singular solution may be inter¬ 
preted geometrically as the envelope of the family of 
integral curves. 
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Example 1. 
Differentiating, 


y = px +a/p 


Dinerenuatmg, p =p+(x _ a/p , )( * p/(te . 

Firstly, dp/* = 0 or p = 0 gives the general integral 

y = Cx +a/C. 

Secondly, the alternative —/P* leads to y=2a/p, whence, 
eliminating p, 

y J = 4 ax. 

Thus the singular 

the lines y=Cx+a/C, anc 101 ab ola together with the 

differential equation consist o th P be confirmed b y 

aggregate of its tangents. This iact y 
considering the equation in the form 

p*x-p»/+a= 0; 

, • « it has real roots only when 

r r>“. szrn 

Sample 2. A geometrical prob^em^smnothnes^leads^ to^ a 

differential equat.on wh^e^JS ^ ^ e^ ^ „, e 

tanyen, a, any ZoftUe ^ at Vl/> on the 

area a-. Since the* respectively, the diSeront.al 

x- and y- axes are x-ylp> J I 

equation is t> , 

lx -ylpM-P*) = *“'- ° r (!, - pX) P 

This is equivalent to the Clairaut equal,on 

y =px ±ay/( -2 p), 

whose singular solution, given parametrically by 

* = *«v/(-l/2/». ±i«v/(- 2 P). 

i9 = 4a v Thus the curve sought is a rectangular hyperbo a. 
20. Generalisation— ^ 1|lore general 

Clairaut equation is a spccia 
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(VAlembert equation (also knowD as the Lagrange equation) 

y=*g(p) +f(p) ; • • • ( 20 . 1 ) 

in fact it is so exceptional a case that it will now be excluded 
by the stipulation that g(p) is different from p. Differ¬ 
entiating with respect to x, we have 

p-g{p)+{ty'(p)+fip)}fc m • • ( 20 - 2 ) 

Since g{p) does not cancel out with p, we cannot have 
dp/dx = 0 identically, and therefore the general integral 
curve is not a straight line. Nevertheless, the equation 
g(p) -p = 0 may have real roots; if p = m is one such root 
dp/dx will be zero, and (20.2) will be satisfied. Thus 
among the integrals of (20.1) there may be some of the 
linear form y=xg(m) where m is any one of the 

real roots of g[p)-p = 0. It will be shown later (§20, 
Note 2) that these are singular integrals, and not merely 
special integrals, of the equation. 

Since dp[dx is not identically zero, (20.2) may be 
divided throughout by it, and thus rewritten 

dx 

{9(p) - p}^ + g'{p)x +/'(P) = 0. 

This is a linear equation for x\ if ' 

log * (p)= fcfc- 

cf)(p) is an integrating factor, and the general integral is 

{'AP) ~P}<f>{p)x=C-ff(p)<f>(p)dp. . (20.3) 

Thus x is expressed in terms of p and the arbitrary constant 
C\ by eliminating x between (20.1) and (20.3), y may be 
similarly expressed, and thus the general integral may be 
obtained in terms of p, regarded as a parameter. It is 
only exceptionally that the parameter can be eliminated 
and the general integral obtained as a single equation 
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general integral. 

Example. 

y =x+p* - 3 P • 

This is of the above form; differentiating with respect to x. 

p = l +(2p - 2 p')dpldx 

° r (p - \){2pdp/dx + 1) =0. 

The alternatives are therefore 

p = 1 and 2pdp + dx = 0. 

thst 

From the second alternative wo obtain 

and, deducing y from the original equation, we may express the 

general integral as . , 

X =0-P\ y=C-\p'- 

In this case P can be eliminated, to give the general mtegral 
in the form 9(0 _ y) . = 4(0-*)■. 

this line with the particular curve (C ) 

Qy* + 4x s =0; 

the absciss ee are given by + ^ + ^ =fl _ 

<e “ 4x» + Ox 1 + Gx + 1 = 0 or <x + l)«(4x + 1) =0. 

Thus the line touches the particular^ curvo^C 0 f° t ho family are 

abscissa - 1; since the other in g . parallel to the 

derived from this Part-u'ar one by translal»» P overy 

line y =x (which leaves the lne y - x+t Uieir 

curve of the family touches the hno y x x t. 

envelope. 
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21. Further Generalisation. The integration of any 
differential equation that can be solved algebraically for 
y in terras of x and p, thus 


y=f[x,p), • • • ( 21 . 1 ) 

may be attempted by the same device of differentiating 
with respect to x. This gives 


df df dp 
^ dx + dp dx’ 


. ( 21 . 2 ) 


which is an equation of the first degree in dp/dx; its general 
integral 

p = <f)[x,C) or x = ip(j),C), . . (21.3) 

as the case may be, associated with (21.1) furnishes a 
parametric general integral of the latter. On the other 
hand, if p - dffdx and df/dp have a common factor involving 
x and p, (21.2) will be satisfied by equating that factor to 
zero, an equation which may furnish a singular integral 
of (21.1). 

An equation soluble for x, i.e. one that can be expressed 
in the form 

* =f(y, P), 

may bo attacked in a similar manner by differentiation 
with respect to y. 


Example. 


x (x + 1 )e x 

y =— p + -. 

x + 1 p 


Differentiating with respect to x: 


p (x + 2)c x j x (x + l)c x Wp 

^ (x + I)» p + lx + l p % Jdx* 
i.e. 

j X (x + ] )f x l jdp x + 2 \ 

tr + -i--p= /te'^v'r 0 - 

The second factor gives a differential equation whose general 
integral is 

p =C(x + l)e*; 
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substituting for p in the original wo have its general integral 

y = Cxe x + 1/C. 

The first factor gives 

p* =(x + 1 )*e x /x. 

On the one hand, integration of this gives y =2 y/(xe*) + const.; 
on the other, substitution in the original equation gives 

y x = 4xe z . 

This is therefore a singular solution; it may bo verified to 
represent tho envelope of the integral curves. 

22. Equations with One Variable Missing. If an 
equation of either of the forms 

F(x,p) = 0 or G(y,p) = 0 

is soluble for p, thus 

p=f(x) or p=g{y), 

integration by quadratures is immediate. 

On the other hand, we may consider cases in which 

solution for p is impracticable, but solution tor x, or y, 

is possible. If, for instance, we have 

x = <j>(p) or 2 /=«A(p). 

the second variable (y or x) is expressible in terms of p by 

quadratures. Thus if * = <f>(p) we have 

y - C = fpdx = fp<f>'(p)dp 

or alternatively ... 

y _C =/y,dx =px - [xxlp =px - $q>(p)dp, 

and if y = i fi(p) we have 

x-C=fdylp=W'{p)dplp 

01 x -c=\iy\p - yip + SydrlP* - ylP + WP^vlv % - 

Another method is to replace the given differential 
equation by its equivalent in terms of a parameter t. 
Thus let G'(y, p) =0 be equivalent to 



48 


DIFFERENTIAL EQUATIONS 


Then 

x - 

which, taken with y=<f>(t), gives a parametric representation 
of the general integral. 

Example 1. 

x * =/>*(<** - x s ). 

Solving for x, 

x =apf y/(l+p')} i.e. dx=adp/(l +p t ) 9,x . 

Hence 

pdp a 

(1 +p*) s/,= "TfT+Tj* 

In this case p may be eliminated, giving the general integral 

x x +(y -O)* =o*. 

Example 2. 

y=p'/(p + \). 

x-C =jdy/p =y/p + !y dp/p' 

= y/P +fdp/{p + 1) =y/p +log | p +1 |. 

Example 3. 

y 1 +p 1 =a*. 

This is equivalent to 

y=aBiat, p=acoat, (dy =a cos t dl). 

Hence 

dx=dy/p=dt, i.e. x - C =t. 

The general integral is therefore 

y =a sin (x -C). 

23. Homogeneous Equations. Just as a homogeneous 
equation of the first degree can be reduced to integrable 
form by the substitution y=vx, so also can an equation 
of the form F(x, y , p)= 0, where F is homogeneous with 
respect to x ami y. For if the degree of homogeneity is m, 

F(x, vx, p) =x m F( 1, v, p) =*x m Q(v, p) say, and the equation 
become* 




G(v, p) =0. 
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If it can be solved for p, say p=4>(v), we have 

dv Xt \ 

V+X H = ^ V) ’ 

and proceed as in § 4. On the other hand, if it i9 soluble 
for v, say v = ifr(p), we write 

y —vx = xifj(p), 

whence, differentiating with respect to x, 

dp 

P=* P<J>) +Xt P'^dx> 

i.e. dx _ if>'(p)dp 

x ~p -ipipY 

and, integrating, we obtain x and thence y in terms of the 
parameter p. 


Example. 

x* - * v + 2 *yp + ( 1 - k * V = °- 

This equation is homogeneous in x and y; writing y =vx, 
cancelling the factor x*. and solving for p, we huve 

dv -(k'-l)v') 

x Tx +v=p= l 

The variables v and x are separable; writing 1 - (A:* - l)o*=z* 
we find 


which leads to 


t.e. 


— + ^=o 

z + k x 


x(z +k)=C, 
x*z*=(C-kx)' 


x* -(Ar* - l)y* =(C -kx)'. 

This general integral may be written 

x*(l - A:*) +2/*(l -k*) +2kCx-C*=0 

representing a family of circles. As C is quite arbitrary 
nothing is lost by regarding k as positive throughout. 
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24. Geometrical Interpretation of a Differential 
Equation. It will be assumed that F(.r, y, y') is a poly¬ 
nomial of degree m in y', and that each coefficient in the 
polynomial is a one-valued function of x and y. 

Let us replace y' by z and interpret x y y,tzs> co-ordinates 
in space referred to rectangular axes, with the 2 -axis 
vertical. The (x, y)-plane will then be spoken of as the 
horizontal plane. Consider any point A of co-ordinates 
(;r 0 , y 0 ) in this plane, such that the vertical line through A 
intersects the surface F(x, y , 2 ) =0 in at least one point P\ 
let the height of P above the horizontal plane be z 0 . Thus 
the point A has a definite direction y' =z 0 associated with 
it; if ^4 begins to move forward in this direction, P will 
begin to move along the surface, but the altering value 
of 2 will involve a change in the direction of motion of A. 

If, then, we suppose that A moves along the horizontal 
plane in such a way that the direction of its motion is 
measured by AP (y'=z ), the path traced out by A will 
be an integral curve of the equation, for it will be a con¬ 
tinuous curve such that at every point the relation 
F(x,y, y')= 0 is satisfied. This integral curve is the 
horizontal projection of a certain curve on the surface 
F(x,y,z)= 0, such that the condition y'=z or dy=zdx 
is satisfied throughout its length. 

In passing, it may be noted that in the case of the 
Clairaut equation y =px +f(p), the non-singular integral 
curves are the projections of the intersections of the 
surface y=zx 4 -/( 2 ) by the family of parallel planes 2 =C. 

r l o return to the point A (r 0 , y 0 ). A vertical line through 
.-1 will in general intersect the surface not in one, but in 
several (though at most /a) distinct points. With each 
ot these points is associated a value of z, its height above 
the horizontal plane, and thus a definite direction y ' at A. 
Consequently, when the points P on the vertical through 
A are simple intersections with the surface, i.e. when the 
roots ot the equation F{x Q , y 0 , z) = 0 are all distinct, 



EQUATIONS OF HIGHER DEGREE 


51 


integral curves, all with different tangents, will pass 
through A, equal in number to the real roots of this 

equation through ^ ^ ^ y#) of the horizontal plane 

there will pass at most m integral curves with distinct 

tan Now consider a point A(x 0 , y 0 ) in the horizontal plane 
such that the equation F(x 0> y 0 , z)=0 has a repeated 
root. Then we should naturally expect two or more o 
the integral curves that pass through A to have the same 
tangent at A. We shall investigate this case more closely, 
with a view to ascertaining whether or not this supposition 

agrees with fact. 

The analytic condition for the equation 

F(x 0 ,y 0 ,z)=0 

to have a double or multiple root is that 

!/o- 2 ) ==0 

simultaneously with it. The geometrical condition is that 
the vertical line through A shall meet the surface in two 
or more coincident points, i.e. shall be a tangent line to 

thC LeW^ move so that the vertical line continues to touch 
the surface; that is to say, in such a way that the equations 


F(x,y,z)~ 0, ~F(x,y,z)=0 


(24.1) 


are simultaneously satisfied. Then the path of A will 
he the trace on the horizontal plane of a vertical cy in 
enveloping the surface. The most obvious cylinder is 
that which touches the surface along what is known in 
descriptive geometry as the visible outline in h °™ ont “‘ 
projection; its trace forms a natural boundary to the 
family of integral curves, but, as will be seen is not 
necessarily an integral curve itself. There may be other 
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vertical cylinders touching the surface along a curve; 
their traces will generally cut across the family of integral 
curves. All possible traces are included in the eliminant 
of (24.1), which represents a curve or curves, known as 
the discriminant locus. 

For any displacement ( dx , dy, dz) on the surface, we 
have * 

F x dx + F v dy + F z dz = 0, 

but since at any point on the curve of contact of an 
enveloping cylinder, F z = 0 (24.1), it follows that for any 
such displacement originating on the curve of contact 

F x dx + F v dy= 0. . . . (24.2) 

For a displacement along an integral curve, we have 

-zdx + dy = 0. . . . (24.3) 

If the horizontal projection of the displacement on the 
surface, at a point (.r, y, z) on a curve of contact, is along 
an integral curve, (24.2) must be consistent with (24.3). 
Consequently, either dx=dy = 0, which implies that the 
displacement on the surface is vertical, or 


that is to say, 



dF dF 

dx + ~dy 



in conjunction with (24.1). 



. (24.4) 


25. Cusp on the Integral Curve. Let us suppose that 
the equation of the integral curve of F(x, y, y')= 0 that 
passes through ( x 0 , y 0 ) can be expressed in terms of a 
parameter t, thus 

x=x(t), y = y{t). 

* A suffix denotes partial differentiation with respect to the 
variable in question. 
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This integral curve will be the projection of a curve 

x=x(t), y = y(0. Z=2 W 

on the surface, with the condition dry= zdx., 
no loss in generaUty if we assume that the point {x 0 y 0 , o) 
on the curve of contact corresponds to the value t 0. 

We shall consider the case of a vertical .displacement 

on the surface from the point (x 0 , y 0 , z 0 ), t.e. d J 

when t = 0; hence 

x-x Q = \x'(0)t 2 + kx m (0)P + • • • 

y - y 0 = y W* + fc/"'(°)* 3 + • • • 

z - z 0 =z'(0)t + • • • 

The approximation to the integral curve in the neighbour¬ 
hood of (*„, </,) is therefore 

x'(0)(y - ».) - V'(0)(x - x 0 ) = «*-(0)y-tO) - y (0)x (°))« + ■ ■ • 

= K (x - Xq)*' 2 + • • • 

where K is a constant. Thus the integral curve has a 

. T^oint (x Va) on the discriminant locus, and 
cusp at the point (x 0 , y 0 ) u "(C\\lx m t 01 

1 .. .. _ tanccnt at the cusp is y 

& rr 4 -- rsr 

g^u^la^so'thd the integral curve does not 
touch the outline. Thus in this most general case, the 

a'rict Vl^cSusscd briefly in a later 
section (§ 27). 

26 Envelope of Integral Curves. Suppose the condition 

dx =dy= 0 is not satisfied, which (24.4) implies tha 

F x +zF v =0 

at all points of the curve of contact on the surface. But 

we also have the condition (24.2) 

F x dx +F v dy =0, 

where d y! ix ia the direction of a displacement along the 
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discriminant locus. Therefore this value of dy/dx may be 
identified with 2 in F(x, y, z) = 0; that is to say, the portion 
of the discriminant locus under consideration is an integral 
curve of the differential equation F(x, y, y') =0. Thus 
when the curves on the surface which project into 
the integral curves do not cross the curve of contact 
vertically, the latter curve projects into a branch of the 
discriminant locus which is itself an integral curve. If 
(, x , y, z) is a point on the curve of contact, (x, y) will be a 
point both on the discriminant curve and on one of the 
general family of integral curves, and these two curves 
will have a common value of y' at the point, i.e. they will 
have a common tangent. Thus at every point on the 
discriminant curve in question it will be in contact with 
one of the general integral curves, and will therefore be 
an envelope of the family of integral curves. I 11 very 
special cases the envelope is itself a member of the general 
family, but in general it is not; it is then know-n as a 
singular integral curve , representing a singular solution of 
the equation. To sum up, reverting to the p-notation : 

The equation obtained by eliminating p between 

F(*,y>P)= 0 and F v {x,y,p) = 0 . (26.1) 

is known as the p-discriminant equation ; it represents a 
curve (the p-discriminant locus) of one or more branches 
in the plane of the general integral curves. A branch 
tvill be a singular integral curve, i.e. an envelope, if and 
only if the further condition 

F x {x, y, p) +pF v (.r, y, p) =0 . . (26.2) 

holds at every point. Otherwise the branch is a locus of 
special points, generally of cusps, on the integral curves. 

Note 1. In the case of the Clairnut equation px-y +f(p) =0 
or the generalised equation f(p) +g(px -y) =0, condition (26.2) 
is identically satisfied, so that the p-discriminant furnishes the 
singular solution and nothing elso. 
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Note 2. In the case of f(p) +xg(p) -y= 0 (20.1) the second 
of equations (26.1) gives f'(p) +xg'(p) =0, which implies (20.2) 
g{p) —p =0. But this is precisely the condition imposed by 
(26.2), showing that the singular integral curves are straight 
lines whose gradients are given by the real roots of 
g(p) -p=o. 


27. Equation of tlie Second Degree. Lot (x 0 , y 0 , z 0 ) be 
a point on the curve of contact of the surface F(x, y, z) =0 
with its vortical enveloping cylinder, so that z=z 0 is a 
double root of the equation F(x 0 , y 0 , z) =0, which implies 
that F(x 0> y 0 , z) will have the factor (z-z 0 ) 2 . For values 
of x, y sufficiently near to x 0 , y 0 , F(x, y, z) will have two 
factors z— z lt z—z 2 , whore z lt z 2 are functions of x,y which 

both become equal to z 0 when x=x 0 , y=y 0 - T,1US 
factor (z -z 0 ) a is the limit, as x -*■ z 0 , y -> y 0 of a quadratic 
expression z 2 -(z t +z 2 )z +z 1 z 2 , and therefore, in the neigh¬ 
bourhood of the ^-discriminant, the differential equation 
F ( Xi 1Jt p ) _o may bo regarded as approximated to by 
a quadratic differential equation p 2 -(z x + z 2 )p +z x z 2 =0, 
and thus consideration of such an equation will confirm, 
and possibly supplement, the discussion of the preceding 

sections. . 

Wo consider, therefore, the following general typo ot 

equation of the second degree 

p*L(x, y) —2pM(x, y)+N(x, y) =0, . (27.1) 

and wo shall assume that L, M, N are functions develop¬ 
able as ascending power series in x and y. Solving tor p, 
we have 

M ± x'{M 2 —LN) (27.2) 

P — i 


Thus the (x, y)-plane 
namely : 

(a) regions for which 
curves exist; 


is divided into distinct regions, 
M 2 < LN in which no integral 
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(6) regions for which M 2 > LN, where two distinct 
real values of p exist for every ( x , y ), i.e. two integral 
curves with distinct tangents pass through every point. 
The two values of p are equal for every point at which 

M 2 -LN =0. 

Now this equation represents a curve which, however, 
may be composed of several distinct branches. Thus it 
includes 

(c) the frontier between the regions (a) and (6) char¬ 
acterised by the fact that passage across the frontier 
involves a change in the sign of M 2 - LN ; 

(d) curves situated within the regions (a) or (6) such 
that M 2 - LN vanishes without a change of sign. 

Cusp Locus. Let the origin 0 be moved to a point on 

a branch T of the discriminant locus which is a frontier 

between regions (a) and (6). Then in the neighbourhood 
of 0, 

L = I 0 + 1 l x + L i y + . . . 

M =7n 0 + ?n l x +jri 2 y + . . . 

N =n 0 + n x x + n.,y + . . . 

and ™ 0 2 = 1 0 n 0 . 

I hen the gradient of the integral curve that goes 
through O is (27.2) 

V = = ”o/ w o- 

But near the origin 

~ LN = m 2 + 2m 0 (m l x + in^y) + . . . 

“ V*o “ U n i* + *&) - njliz + l,y) - . . 
so that the tangent to V at the origin has the equation 

(2m 0 w?, - - n Q l x )x + (2m 0 m 2 - l 0 n 2 - n Q l 2 )y -0 

or 8a y ax+by=Q. 

Since M--LN changes sign we may assume that one 
at least of u and b is not zero. Since - a/6 differs in general 
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from m 0 H 0 , the slope of T at O differs from that of the 
integral curve, so that T is not itself an integral curve. 
To find how the integral curve through O behaves with 
respect to f\ we require a further approximation to (27.2) 


V = 


m 0 + m x x + m^y + . . . ±\/(ax + by + . . .) 


1 0 +l x x+l t y + . . 

The first approximation is p = m 0 /l 0 \ we therefore write 
y=vi 0 x/l 0 + Y, where Y is of higher degree than the first 
in x, and obtain an equation of the form 

dY ax + . . . ±'y/(/fc + • • •) 
dx Jo "*■ • * * 

where a, p, y are constants. Integrating, we have 


= + yx 1/2 + . 


y = ± ijy^ 3 ' 2 + • • • 


or 

y = m 0 x/l 0 ±lyx a ' 2 + . . 

so that the integral curve has a cusp on T. Hence a.ny 
branch of the discriminant locus that separates a region 
of existence from a region of non-existence of integral 
curves is in general a locus of cusps on the integral curves. 

Envelope. Now consider the special case when 
- a/6 - mjl 0 . When the substitution 

y = m^Ho + Y = - ax/b + Y 

is made, the term in x under the radical disappears, leaving 
dY a x -f . . . ± x /(cx 2 +bY +.■•) _ 
dx lo + • • • 

It will be found impossible to solve this equation by 
y = Ax" + . . . with p < 2; hence it is of the form 

y'-(/3±y)*+ • • •> 

which gives 

y = rn 0 a:// 0 + i(^±y)2: 2 + • • •> 

showing that two distinct integral curves touch one another 
at O. When this happens continuously along a branch 
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of the discriminant curve, since the slope of the latter is 
everywhere the slope of an integral curve, it is itself an 
integral curve. So in this case one of the curves having 
contact is a branch T of the discriminant curve, the other 
is one of the family of integral curves, which are therefore 
enveloped by I\ 

An application of (26-2) at 0 gives 

(;p 2 /, - 2pm 1 + »j) +p(p 2 / 2 -2pm 2 + a 2 ) =0. 

If each bracket is divided by p and if p is replaced by 
n 0 /m 0 and 1/p by lJm Q this reduces to a+pb=0, and 
therefore if (26.2) holds at any point of T, T will there 
touch an integral curve; if it holds at all points of T, the 
T will be an envelope of the integral curves. 

Tac-Locus. Now suppose that 0 is moved to a branch 
T of the discriminant curve for which M 2 - LN vanishes 
without changing sign; the sign on either side of T will 
be assumed to be positive. Then in the neighbourhood of 0 

M 2 - LN ={ax + by) 2 + . . . 

The form of the integral curve through 0 is given ap¬ 
proximately by 

_ wt 0 + m t x + m 2 y + ... ± \/{(ax+by) 2 + . . .} 

/ 0 + I x x + hy + . . . 
from which we find that. 

y =m Q x/l Q + (a± fi)x 2 + . . . 

Thus two distinct integral curves touch one another on the 
curve r. Since the same is true at all its points, T is a 
tac-locus. In general -a/6 is not equal to m 0 fl 0 and 
therefore the tac-locus is not an integral curve. 

A'o/e 1. The above assumes M* - LN to be positive, and 
the integral curves real on either side of the tac-locus. In the 
contrary case when A/* - LN is negative, there will be a real 
tac-locus of imaginary integral curves. 

Note 2. Since the approximation to the tac-locus, at any 
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point on it, is a squared term, the tac-locus occurs ns a squared 
factor in the p-discriminant. 

Example. 

xp* + (y - 3 x)p + my = 0. 

This equation may be solved for y and integrated by the method 
of § 20. The p-discriminant equation is (y-3x) 2 - 4mxy =0; 
the corresponding curve consists of two straight linos through 
the origin which are real except when 0 < m < - 3, and are 
coincident when m=0 or -3. The condition for a singular 
integral (26.2) gives 

p* - 3p + p(p + m) =0. 

The first root p =0 would lead to y =0, which is no part of the 
p-discriminant; the second root 2p = 3 - rn gives 2 y =(3 - m)x. 
On substituting this in the p-discriminant equation, wo find 
m = l (twice) to bo the only admissible case, giving p = 1, and 
note that y =x is a solution of the differential equation. 

Integrating the differential equation for m = 1 in terms of 
the parameter p, we have 

**= C(p +1 y/p*. y'=C( P -3yip 

or, eliminating p between this expression for y' and the 
differential equation, wo obtain the general integral 

(xy* + Cy + ZCx)(y* + 16 Cy - 27 Cx) + C’(i/ - 9x)» =0. 

We find that y =x has, in fact, contact with all integral curves 
for which C > 0, touching each at two distinct points x =y 
= ± y/O. This contact at two points accounts for the double 
value of m; it is as if the family of integral curves had two 
coincident envelopes. Now consider the other branch, 
y = 9x, of the p-discriminant locus; squaring and substituting 
the above parametric values of x, y we find that the va ues 
of p at the points at which y =9x meets any integral curve are 
given by 

(p* - 12p-9)(p*+0p+9)=0. 

The two single values from the first bracket are of no interest, 
they merely correspond to simple intersections. The double 
value p = - 3, however, shows that every integral curve 
encounters y - 9x in points x* = - 4C-/27 at which it has a 
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double tangent. To investigate the behaviour of the integral 
curve near such a point we put p = -3 +t, where t is small, in 
the parametric expressions for x and y. If k=±y/( -27/4 C) 
we find 

x=C l/ '( - 2 +/)( - 3 
or 

= 1 — tj-'j /* — J ^ t 1 + . . . 

(3 P -P')/(P + 1) =9(1 - .U+ r V J )(l "JO" 1 

= 9(1 + T V*+ . . .) 

fy=te(3p-p*)/(p + l)=9 + £/*+ 2 y* + . . . 

Hence 

3(kx - 1) +(ky - 9) = -+ . . . 

and 

{3(fcc-l)+(Ary-9)}*=5 3 a(tc-l)« + . . . 

Thus every integral curve for which C is negative has a cusp 
on the lino y =9x, the double tangent at the cusp being parallel 
to 3x+y = 0. This is an instance whore a branch of the 
p*discriminant locus is a cusp-locus. 

To illustrate the circumstances that arise when the p- 
discriminant locus has a double line, take the case m= -3. 
The general integral is now 

(xy -C)(y -3x+C)=0 

and thus consists of a family of rectangular hyperbolas and a 
family of parallel straight lines taken together. The p- 
discriminant locus is the double line (y +3.c)» =0. Now even- 
one of the parallel straight lines touches one hyperbola of the 
family; in fact t/ = 3x-6y touches xy-t-3y*=0 at the point 
x =y» y = - 3y which is a point on the discriminant line. This 
is a tac-point, or point of contact of two distinct integral 
curves, and the double line of the p-discriminant is a tac-locus. 
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EQUATIONS OF THE SECOND AND HIGHER ORDERS 


28. Reduction of the Order of an Equation. When we 
turn to equations of higher order than the first, we discover 
that there exist certain well-defined types which admit 
of transformations whereby the order may bo lowered. 
In particular, there are equations that may he reducec 
by a transformation of the dependent variable, to allied 
equations of the first order; if the latter can be integrated, 
a reversal of the transformation (which usually amounts 
to one or more quadratures) will lead to the general integral 


of the former. . 

There is one small point that may be noted in passing, 

namely that whereas in an equation of the first older 
e.g. Pdx+Qdy =0, there ia no natural discruninati 
between the variables, in equations where second and 
higher derivatives of one variablo (y) with respect t 
other (*) occur, the distinction is evident In cq.urt iciiis 
that arise from problems in physics and mechanics, in 
particular, this distinction arises from a difference m the 
character of the variable; as when one represents a length 

and the other represents time. , 

We shall deal mainly with equations of the second oicle 

and shall consider, in particular, those ola.ses .n which a 

change of dependent variable enables a reduction to the 

first order to be effected. No general discussion of higher 

orders will be attempted, but when t 10 sc<>l ' c ... ^ 
process extends beyond the second order, the tact 

mentioned. 


61 



62 


DIFFERENTIAL EQUATIONS 


The most general differential equation of the second 
order in which x is independent, and y dependent variable, 
may be written 

F(x,y,y\y')= o. . . . (28.1) 

If it may be derived from an equation 

f(x,y,C,C 1 )=0 . . . (28.2) 

in which C, C 1 are arbitrary constants, by differentiating 
twice with respect to x and eliminating these two constants, 
(28.2) is known as its primitive. The elimination of C and 
G \ may be performed in either order, but can only lead to 
the one equation (28.1). For if two such equations were 
found, y could be eliminated between them, leaving an 
equation of the first order having the primitive (28.2), 
which cannot be the case unless C and C x are specially 
related. Integration consists of recovering (28.2) or any 
equivalent expression containing two arbitrary constants, 
which is a general integral. 

29. Equations that do not Involve y. The simplest 
case where reduction of order is possible occurs when the 
dependent variable y itself is absent from the equation, 
which may be written 

F(x,y\y')~0 . . . . (29.1) 

We replace y' by p and regard p as a new dependent 
variable, temporarily replacing y, and thus have 

F(*,P,P')= 0, . . . (29.2) 

an equation of tlie first order in p. 

Let us suppose for the moment that this equation can 
be integrated explicitly, thus 

P =/(•*•, C), (29.3) 

introducing one arbitrary constant C. We then obtain 
y by the quadrature 

y~jpdx=ff(x, C)dx + C lt 


. (29.4) 
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introducing the second arbitrary constant C l and thus 

arriving at the general integral. 

If, on the other hand, the integration of (29.2) gives z 

more naturally in terms of p, thus 

x=g(p,C) . • • ( 29 - 5 ) 

we proceed to take the differential 

dx=g'(p, C)dp 

and thus we have 

y = \pdx =jpg'(p, 0)dp + C X . • (29.6) 

or, alternatively, we may write 

y =J pdx =px -fxdp =pg{p, C) -fg(p, C)dp +C X . 

Thus a: and y are expressed in terms of the parameter p. 
When x is absent as well as y, and the equation can be 

written as 

V =/(*/') 

the process is simplified, for we have 

p'=f(p) or dx —dplf(p). 

Also dy =pdx=pdp/f(p) and the general integral is given 
parametrically by 

x=Sdplf(p)+C, y=fpdplf(p)+C i. • ( 29 - 7 > 

In the case of the equation of order n 

yin) =f(y in ~ l} ) 

we may write z=y {n ~ l) and thus obtain as in (29.7) 
x =fdzlf(z) +C, i/'"- 2 ’ =J*fe//( 2 ) +C (- 

and then in turn 


yl f»—3) 


= Jyln-2)dx =f,/ n -*'dzlf(z) 


’ dz f dz 

fwnw 


+ C (J.' + 0 2 


y 


(n-4) 


f dz r dz f dz L 
}f(z))f(z)) Z f(z) 


\C x x* + C 2 x + 
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and so on until we arrive at y, which will involve an 
arbitrary polynomial of degree n- 2 (involving n — 1 
arbitrary constants). 

An equation of the type F(x, r/<">) =0 may also 

be reduced to the first order by writing z =y(n-i). 


Example 1. 

* 2 y'=y' % -2xy'+ 2x*. 

Writing y' =p, we have 

x s p'=p» -2xp + 2x*, 

an equation homogeneous in x and p ; if p =vx it reduces to 

xv' =t’ s -3v +2 

whoso integral is 

Cx = (v - 2)/(v - 1) whence v =(Cx - 2)/{Cx - 1). 

Thus 


and finally 



= vxdx = 


x(Cx - 2) 
Cx-1 



y = Lc- 



- Io e I c* -11 



Example 2. 


2xyt/' -xy m +y'=0. 
Writing y' =p and dividing by x 1 we have 



which is exuct and has the integral 


p'-?+C = 0. 

We thus obtain 

whence p* + C 

V =fpdx =px - jxdp 


P * f p * 

+c ■ - i l0 « I p‘ + 01 

It ia hero possible to eliminate p . 
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Example 3. 

The substitution y m 
order 


xy" - 2y" =r*. 

gives the linear equation of the first 
xz' - 2z = x\ 


whence 


y m —z =x* + Cx* 


and in turn 

y’ = %x* + JCx* + Ci 

y' = 2 l ox 5 + + 

y= r b*‘ +AC** +C..T + C,. 


30. Equations that do not Involve x. VY lien the in¬ 
dependent variable is lacking, the equation takes the 

form i on l \ 

F(y,y',y')=0, • • • * 301) 

and we again write p for y'. But as y, p are now the 
variables involved, y" requires transformation as follows : 

dp dp dy _ dp 
y = dx~dydx~ P dy 


Thus the equation becomes 

F(y, p, pdp/dy) =*°; . • • ( 30 * 2 > 


it is of the first order with y as independent, and p as 
dependent variable. If an integral of the form 

p =f(y, C) 

is obtainable, we separate the variables, thus : 


dy 


/(y. O 

and obtain the general integral 

dy 


dx 


f_«y_ 

J/(y. C) 


+ C 
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On the other hand, when the first integral appears in 
the form rr 

. y=g(p,c), 

we have 


J^JJ + [ydp_yJg(p, i 

h v ip 2 p J ' ~p 


C)dp 


+ 


If y' is absent from an equation of the above type, 
and it can be expressed in the form 

y" =f(u), 

we write y"=pdp/dy and multiply by dy, obtaining 

pdp =f(y)dy. 

Integrating: 

whence P ‘ =2 ^ +0 > 


-pf-f._ Jl _ +c 

)P Jv'{2J7(y)dy + <7} +c i- 


; qUati ° n that i '"' olves three consecutive derivatives, 

te redo d g h 8 f’r ay ^ ^ "°. "V si 'nilarl y 

be reduced by tak.ng y <->> = „, yU-n =z> yl „, =zdz / du _ 

Example 1 . 

y(y - l)y'+ y' 1 = 0 . 

The transformed equation is 

y(y - 1 )pdp/dy +p*= 0 . 

The variables p and y are separable and the fust integral is 
Hence Cp=y/(y_ 1) . 


~lp- c \ s ~jr d » =c *-°'°B 


ly I +c t . 
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y" + m*y =0. 

This becomes 

pdp/dy + m*y =0. 

Separating the variables and integrating, we have 

p x + m*y* = m’C', 

whence 


or 


p = ±mV(C* -y*) 

-—-— ± mdx 

V(C*-y*) 

and integrating, arc sin (y/C) = ±mx +C U 
i.e. 

y =C sin (C, ±mx) 

which may also bo written as y =A cos tux + B sin mx. 
Example 3. 

yln-Dyin) 

The transformation mentioned above gives 

uzdz/du=z x or dz/z =du/u. 


whence 
Thus we have 


z=Cu or u'=Cu. 


"-*> =u=C l e Cx 


and integrating n - 2 times in succession wo obtain the general 
integral in the form 

y =Ke Cz +an arbitrary polynomial of degree n - 3. 


31. First Homogeneous Type. The expression 

F(x, y, y\ y") 

is said to be homogeneous and of degree in in y and its 
derivatives if. when A is any constant, 

F(x, Ay, Ay', A if) =A m F(x, y, y', y”). 


. ( 31 . 1 ) 
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When this is the case, the equation F(x, y, y', y’)=0 may 
he reduced, by division throughout by y m , to one still 
homogeneous in y, y', y ', but of degree zero. Its form 
would then be 

J ?/' y"\ 

f Y'y'y) =0 ' • • • (312) 

The order may now be lowered by taking y'jy = u as 
a new dependent variable, i.e. by the transformation 

y = ef udx 

which imphes y'= uy, y' = y(u' + m 2 ). For the equation 
then becomes 

f(x, u, u' +u 2 )= 0 . . . . ( 31 . 3 ) 

It is now of the first order, but here as in most cases 
the lowering of order is paid for by an increase in com¬ 
plexity. Thus, when the method is applied to the linear 
equation 

t/'+p(%'+5(%==0, 

which is homogeneous and of the first degree in y, y' } y ', 
it becomes 

u' + u 2 +p(a,’)w+g(x)=0, 

a Riccati equation (§ 12) which is actually less manageable 
than the equivalent linear equation. 

\\ hen an equation of order n is expressible in the form 
/( j: > V i'll, • • y {n) ly) =0, it may be transformed by the 
above process into an equation of order n - 1. 


Example. 

X 'J'J ~ x y'* + yy’ = o (second degree in y , y’, y’). 
This muy be written 


y y'Y.y 


- -n - +"-=0. 

y y'y 

The transformed equation is 

a-(u' +u J ) -xu* +u =0 or xu' + u=0. 



equations of the second and higher orders 6» 

Thus we obtain the first integral 

xu=C or xy'=Cy 

and finally arrive at the general integral 

y=C 1 x c . 

32 . second Ho^eou. ^ ^ 

say.ThTxTreplaned by Ax and dx by \dx we have 

F(Ax,y, A-V,A-V)=A"F(x,y, »,!/)• ■ ( 321 > 

The equation can then always be written in the form 

xV)=0- ' • • (3 - 2) 

We make the transformation z-e‘, so that 

_ ».» 1 >< 1 /lit 


dy dyidx 1 dy ^Jd^,_\_dy 

£ = dt\it~xdt' dx* x*dl* x-dl 

and we also write 

d'J ^ = y (as in § 30) 

Tt~ ' <& d y 

bo that (32.2) becomes 

f( y ,v,vdvldy-v)=0, . . 

which is of the first order. expressible in the form 

to an equation of order n - 1. 

Bxampie ' 

Changing the variables as above, we have 


d 2 y dv 
— =VT" 
dt 2 dy 


(as in § 30) 


. (32.3) 


which reduces to 


dv \* „ dv 

— + 2?/- v — u. 

> du 


70 


DIFFERENTIAL EQUATIONS 


If v* is taken as dependent variable, this equation becomes of 
Clairaut type; its integral is 

v*=4(Cy+C*), 

whence, writing v =xdy/dx and separating variables, we have 

dy _ dx 

± 2y/(Cy + C*)~ 7’ 

whence 

± V(Cy + C l ) = C log C x x. 

Rationalising, the general integral may be written as 

y = C{(log <?»*)•-1>. 


33. Third Homogeneous Type. The characters of the 
first and second types are here united; F(x, y, y\ y ') 
is homogeneous and of degree m in x and y even when 
these symbols are associated with the differential operator d. 
Thus y' =dy]dx is of degree 0 and y" =d i y/dx' 1 is of degree 
-1, so that 


F (^ Ay, y', X~Y) = \™F(x, y, y' } y’). . (33.1) 

The equation F(x, y, y\ y") =0 may thus be written 

f(yl*,y',xy')= o. . . . (33.2) 

The transformation y = xu changes it into an equation 
of the second type, which may be dealt with as in § 32. 
Another procedure is to make the double change of variable 

y=xu, v = xdujdx, 

so that 

y' = U + V, 


» / dvdu 

and thus the equation is reduced to 


dv 

v+v r 

du 


f(u, u + v, v + vdvfdu) = 0 
of the first order in u and v. 


. (33.3) 
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The order of any equation of the type 

f{y/x, y', . . ., x n-1 y ,n) ) =0 
may be reduced by a unit in this way. 

Example. 

x*y" + m(xy" - y)* = 0. 

Writing this equation as xy" +m(y ' -y/x)* =0 and making the 
above transformation, we have 

dv dv , . 

v + v — + mu* = 0 or ——h mv + 1=0. 
dx du 

This linear equation has the integrating factor e mu ; we thus 

obtain 

e mu ( 1 +nw)=C, 

i.e. 

du\ d 

1 +mx-— I =C or —, 
dx / dx 

Integrating aguin, we arrive at the general integral 

xe mu = Cx + C , 

or 

e mv/x = C +CJX. 


mu 


( 


(xe mu ) =C. 


34. A Special Case of Homogeneity. An equation 
homogeneous in y, xy' and x 2 y” may be written 

F(xy'ly, xV/y)=0 or xhy" = yf(xy'/y). . (34.1) 

It is both of the first and the second type, and therefore 
also of the third, and may be reduced by taking a new 
dependent variable u where u=xy'jy. Thus we have 

xy' ~ uy, xy" =u'y + {u - l)y', 
i.e. 

x i y" = u'xy + u(u - 1 )y, 

and so the equation, in its second form, becomes 

xa' +u(u - 1) -/(u). 
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The variables are separable, giving 


dx du 

x ~f( u ) -m(w-1) 


and hence 


dy udu 
y ~f{u)-u{u- 1)‘ 


(34.2) 


Integrating, we obtain a parametric representation of x 
and y in terms of u, involving two constants of integra¬ 
tion. Alternatively we may integrate the first equation 
and obtain x = (f)(u, C) = (f>[xy'ly , C). Solving for xy'fy we 
obtain an equation whose variables are separable. 


Example 1. The Euler linear equation of the second order 

xh/ - (a + P - l)xy' + afiy =0, 

where a, p are constants, is of this type. The above trans¬ 
formation leads to 

xu' +(u - a)(u - P) = 0, 

whence, separating variables and integrating, 


or 


u - a 
u-p 


= Cx a ~? 


Cax a - pxP xy' 

'CV'-x* = “ = y ’ 

i.e. 

d( Cx a -x 1 ) dy 
Cx* - x fl y * 

whence y = C l (Cx a -x B ), that is to say the general integral 
may bo written y =Ax a + Bx B , where A and B are arbitrary 
constants. 


Example 2. 


x, jy' =y'(2ry' + ay). 
Taking this equation in the form 



and applying the same transformation, we obtain 

xu' =u{u + u + 1) 
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whence 

dx _ du dy _ du 

x u(u + a + 1) ’ y u+a+l’ 

x*+\ = C --— y=C x (u +a + l). 

u + a + l 

By eliminating u, the general integral may be expressed in the 
form 

y = l/(A + Bx a+1 ). 

35. First Integral. The term first integral has been 
used to denote (§ 30) the differential equation of the first 
order, involving one arbitrary constant, that results from 
transforming a differential equation of the second order 
into one of the first, and integrating the latter. When 
the first integral is obtainable by immediate integration 
of the differential equation, the latter is said to be exact. 
Thus 

2 y'y" + xy' + y=0 
is exact and its first integral is 

y' 2 +xy =>C. 

In the same way, a first integral of an equation of 
order n 

F(x,y,y',y m .»/«">) =0 . . (35.1) 

will be an equation of order n - 1 

F x {x, y, y', y\ . . ., y<»-'\ C) =0, . (35.2) 

where C is arbitrary. Moreover, when dFJdx = F, then 
F =0 is exact. 

As in the case of equations of the first order, when 
an equation of any order is not exact as it stands, it may 
be rendered exact by the introduction of an integrating 
factor. For example, the equation 

v" + P(x, y)y' +Q(x, y)y' 2 =0 
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admits of the integrating factor 1 /y 1 whenever Pdx + Qdy 
is an exact differential, and the first integral then is * 

log I V' I +S{Pdx+Qdy) =C. 

A differential equation of the second order has two 
distinct first integrals. For let the primitive be 
f(x, y, A, B) =0, where A, B are distinct arbitrary con¬ 
stants (distinct in the sense of not being replaceable by 
a single constant C, as would be the case, for example, 
if A, B occurred only in the combination A +B). By 
differentiating the primitive we obtain an equation 
<f>(x, y, y\ A, B)= 0, and then by eliminating B and A 
in turn between / = 0, cf> = 0 we arrive at 

*i(*» y> y'> A ) =0, F t {x, y, y\ B) =0 

which are formally distinct because A and B are distinct 
in the primitive. By differentiating again and eliminating 
A, B between the two resulting equations and the primitive 
we reach the differential equation 

F(x, y, y\ y") =0. 

Integrating the differential equation once results in 
the recovery of either F x = 0 or F a = 0 or an equivalent. 
Thus F x =0, F 2 =0 are first integrals, and we see that two 
distinct ones exist. If, by adopting two separate methods 
of integration we recover both F x = 0 and F t = 0, we obtain 
the primitive by eliminating y' between them. There 
cannot be more than two distinct first integrals, for if 
F 3 =0 were a third, a distinct primitive would be obtained 
by eliminating y' between F x = 0, F z = 0. But the existence 
theorems prove, what we shall here regard as a postulate, 
that there cannot be more than one distinct primitive 

* I he equation is also integrable (i) when P and Q are functions 
of x alone, for it is then equivalent to a Bernoulli equation in p , (ii) 
when P and Q are functions of y alone, for then it may bo reducod to 
a linear equation in the variables y and p. 
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dependent upon a set of arbitrary constants equal in 

number to the order of the equation, i.e. there cannot be 

more than one distinct general integral. 

Therefore an equation of the second order has just 

two distinct first integrals. 

Example. The equation 

y" + m*y =0 (§ 30, Ex. 2) 

has, among othem, the following integrating factors : 

cos mx 9 sin mx> 2y 

which lead to the corresponding first integrals : 

y' co 9 mx +my sin mx =A, 
y' sin mx — my cos mx = B, 
y' x +m*y* =C *. 

By eliminating y' between the first two we obtain the general 

integral m y = A sin mx + B cos mx. 

The third may be obtained from the tat two 1g^g.“ d 
adding; the relation between the constants \9 A + 

36 Problems Involving Curvature. The geometrical 
problems previously discussed (§§ 14, 15) depended upon 
^relationship between the slope y' and the co-ordinates 
lx y) which was imposed at every point of a plane curve. 
Such problems thus led to differential equations of the 
first order. We now consider problems in which 
curvature at every point obeys a specified law 
The radius of curvature at the point (x, y) 

( 1 + 2 /' 2 ) 3/2 . 
pmm y" * 

the centre of curvature (f r,) relative to that point >• 

giveDby y'(l +y” 1 ) *+£ 

t-*- —- r >-y + y - • 
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The locus of the centre of curvature is the evolute. A whole 
family of parallel curves have the same evolute; they are 
its involutes. 

As the above expressions involve y a problem that 
depends essentially upon curvature will involve the in¬ 
tegration of a differential equation of the second order. 

As an example, we consider Curves whose Radius of 
Curvature is 'proportional to the Normal. 

The length of the normal is understood to be the length of 
that portion intercepted between the curve and the x-axis, 
which is given by t/(l Thus if the normal is n times 

the radius of curvature, we have 

y(l+y'*)«/*= ± w ( — - ffl* 

y 

or 

yy'= ±n(l + 

The ambiguous sign arises from the fact that each of the 
two lengths in question may be measured in either of two 
directions. Lot us agree to measure them both away from 
the curve. If the curve were, at the point considered, above 
the x-uxis, and concave upwards, both y and y" would be 
positive, but the normal and radius of curvature would run 
in opposite directions. Hence our convention demands the 
negative sign. 

The equation to which the problem has been reduced is 
therefore 


yy’= -n(l + »/'■). 


Since x is ubsent, wo write (§ 30) 

y' =P> y 0 =pdp/dy 

and obtain 

ypdfi= -n(l +p')dy. 

Separating the variables and integrating, we obtain 


1 +p' 
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whence 



where a and b are the constants of integration. 

We take two cases whore the radius of curvature has the 
same direction as the normal (n = 1, J) and two cases in which 
the directions are opposed (n = - 1, - $)• 

When n = 1, 

f- ydy - +f)= _ V{a *- y *)+b 


or 


(x -b)* +y*=a*- 


Thus any circle with its centre on the x-axis 
condition; in fact the normal and the radius are 


satisfies the 
coincident. 


When n = £, 



Write y=a sin* 0, dy=2a sin 0 cos 0: then 

x-b=2af sin*0 (10 = aj(l - cos 20)d0 =a(0 - i sin 20) 


or, if 20 = <£, 

x =b + -sin <f>), y = \°( 1 -<*»*)- 

Thus the curve whose radius of curvature is twice its normal, 
and in the same direction as the normal, is a cycloul. 


When n = - 1, 

x -[ _- + 6 = a urg cosh -+ 6 


or 


x-b 

V =a cosh- 

* a 


Thus when the radius of curvature equals the normal, but is 
oppositely directed, the curve is a catenary. 
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When n = - 


or 


x =^yj +b =2Va[y-a)+b 

(x - b) 1 = 4a[y - a). 


So when the radius of curvature is twice the normal, but is 
oppositely directed, the curve is a parabola with its axis 
perpendicular to the z-axis, and its latus rectum equal to four 
times the ordinate of its vertex. 




CHAPTER V 


LINEAR EQUATIONS 


+p n ~i(*)£ + p»wy =fix) - (371) 


37 Form of the General Integral. An equation of 
order n is said to be linear if it is linear in the dependent 

variable y and the derivatives y , y f . • • • J * 

the most general Unear equation of order n is of the 

d ny d n ~ x y ■ Jy 

Po (z fe +Pl(a:) ^ ri ' * 

Tf D denotes the differential operator d()ldx, the 
left-hand member is derived from , by the eompound 

operator 

Po (x)Dn +Pl (x)D--'+ . ■ ■ +P,- 1 (X)Z3+P„(*). 
and if L is written for this operator, the equation assumes 
the convenient abbreviated form 

Ly =/(*)■ • • • (372) 

We take it for granted that since the equation is of 

.... - »*- 

of (37.2) are known, say y = 2/i and V y '** 

Lyi =/(^)» “/( z )» 

Now since p(*)I>'<P,-*>■-P<>W -***»’ “ d 
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Ly 2 - Ly x is merely a sum of terms like this, we have 

i =0 * 

Thus if u represents the difference between any two 
solutions of (37.2), u will satisfy the reduced equation 

Lu= 0, . . . • (37.3) 

which contains no term free from u or a derivative of u. 
Thus we are led to consider, first of all, the general solution 
of this reduced equation. 

In the first place we observe that if C is any constant, 
p(x)D r (Cu) = Cp(x)D r u, and remembering the composition 
of the operator L we conclude that L(Cu ) = CLu. So if 
u = is any solution of Lu =0, we have L[C 1 u l )=C l Lu- l =0, 
and therefore C l u l is an integral, whatever the constant 
C 1 may be. 

Next we note that if we have a number of particular 
integrals w lt u 2 , u 3 , . . . then whatever the constants 
C„ C 2 , C 3 , . . . may be, 

L(C x u x + C 2 u 2 + C 3 U 3 + . . .) e C^Lai\ + C 2 Lu 2 + C z Lu 2 + . • • 

that is, C x u x + C 2 u 2 + C 3 U 3 + . . . is also an integral. 

Functions « lt u 2 , %, . . . are said to be linearly in- 
dependent if it is impossible to find constants G lt C 2 , C 3 , ... 
not all zero such that 

C x u x +C 2 u 2 + C 3 M 3 + ... =0 

identically. For example, the three functions sin 2 x, 
cos 2 x, sin 2 .r are linearly independent, for no relation of 
the type .4 sin 2 x + B cos 2 x + C sin 2x = 0 links them .* 
But sin 2 #, cos 2 #, cos 2# are not linearly independent, 
for sin 2 # - cos 2 # + cos 2 # = 0 identically. 

Now if u v u 2 , . . ., u n form a set of linearly independent 
integrals of (37.3), then u = CyU x + C 2 u 2 + . . . + C n u n will 
be a solution, and since it is not identically zero for any 

• Assuming the relation, we put x=0 and find B= 0; then 
x = ir/2 and find A =0; it then follows that (7 = 0. 
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, . r C C (not all zero), and also since it 

contains the full complement of n arbitrary constants .t 
will be the general integral of the r educed equation (37.3b 
Any such set of n linearly distinct integrals is said to be 

a/ “Su™"ng now to the original unreduced equation 
(37 2) let y = v be any particular integral, involving no 

arbitrary constant, and let y = Y be the general mtegrah 
Thcn u 1 Y - 71 will satisfy the reduced equation and wil 

contain , arbitrary constants; it will in f«t be the general 
integral of the reduced equation, and we have 

T +CoU 2 + • • . +Cf, u n +r l- 

Thus the general integral of the unreduced equation has 
two components, namely (i) the general Bolut on of the 
reduced equation, involving n arbitrary const:ants^a^ 
known as the complementary function, ( ) V 
integral involving no arbitrary constant. 

Example. , 

x{x , + i,y - (3x* - I)(x* + 1 )»' - I)y -*• - + *■ 

The reduced equation 

x( x« + l)»u' -( 3 x* - l)(x* + IK + 4 x(x* - l)u =0 

has two independent particular solutions 

x* + l, (x* + l)logx. 

Thus ita general solution is 

M= (x* + l)(il + D log x). 

A particular integral of the original is the genera, solution 
therefore is 

y =(x* + \)(A +B log *)+*. 
where A and B are arbitrary constants. % 

rjenression of the Order. When any integral of 

0 

t 
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rdern 1 If the known integral is Uj_, the transformation 
is u -UjJvdz, where v is a new dependent variable. This 
ransformation, and the proof, are general, but for brevity 
we take an equation of the third order 

Po um +PiV+p 2 u'+p 3 u = 0, . . (38.1) 

Substitute Pl> P2 ‘ P3 ^ funcfcions of x > or ' constants, 
g 

u = u> = u i , ! vdx + W u' = u x 'fvdz + 2u x 'v + u x v', 

W" = u/’fvdz + 3 u/v + 3 u x 'v' + v^v', 
and rearranging we have 

(PoV +p x «i' +p 2 V +p 3 v 1 )fvdx 

+ (3p„V + 2ft V +p 2 ft)« + (3p„V + Pl K 1 ) l) ' +?w ' = 0. 

Smce ft is an integral of (38.1) the first term disappears 
leaving a reduced linear equation of the second order in v. 

Example. 

x y ~(2x + l)y' + (a; + l)y = o. 

A particular integral is e* ; writing y=e^vdx, 

y' =e*(JWx +v), y' = e*($vdx+2v + v\ 

we obtain an equation which reduces to 


giving 


xv'-v = 0 

v=2Ax, fvdx = Ax'+B, 

y=e*(Ax'+B). 


a!'„o l„r EqUa ion with Constant Coefficients. 

lo ? s . ln t faults from taking the leading 

coefiicient to be unity, the equation nmy be mitten as * 

(D n +a 1 D n ~ 1 + . . . +n„_,ZJ + o„)tt=0, (39.1) 

where a lt . . ., «„ are constants, or in sliort as 

F ( D )u=0. . . . (39.2) 

but inThfi dfff ^ a not in an algebraic symbol, 

but u the differential operator D; the laws of algebra musi 
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therefore not be applied without justification. Our first 
step will be to show that F(D) can be factorised in the same 
way, and with the same result, as if D were an ordinary 
symbol of algebra. Consider the product 

(D-a)(D-p) 

in which a, j9 are constant numbers, real or complex, and 
let X be any function of x that can be differentiated at least 

twice. Then 


(D-a)(D-p)X=( Tx -a 

d 2 X 


IdX 

dx 


-px 

dX 


- dx » -<“ + %7 +a ^ 

s= {D 2 - (a + fi)D + afi}X. 

But we likewise have 

(D - fi){D -a)X={D 2 -(a + fi)D + aP)X. 

As these relations are independent of X, we conclude 

that * 

Z) 2 - (a + ft)D + a.p = {D- a)(D - p)=(D-fi)(D - a), 

where the sign of equality denotes equality of effect not of 
magnitude. Thus any operator of the second order with 
constant coefficients may be factorised in the ordinary 
algebraical sense, into two linear factors, and these factors 
are permutable ; that is, their order is immaterial. By 
repeating the argument with three factors D a, U p, 
D-y we show that an operator of the third order can be 
decomposed into the product of three linear operators 
taken in any order, and, in fact, should find no difficulty 
other than the labour of writing, in proving that the general 
operator of order n is the equivalent of a product of n 

linear operators. 

„; rvK 5,*-!,™ -f-rns” ts "SS-3 

= -(* + 1 )D+*-l. 
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We shall assume that the coefficients a lf . . a n in 
(39.1) are real, so that the numbers a, y, in the linear 
factors of F(D) are either real or paired as conjugate 
complex numbers. Let D -a be any linear factor of F(D), 
so that the equation may be written as 

F 1 (D)(D-a)u = 0, 

where F X (D) is the co-factor of D - a in F(D). Now this 
equation is certainly satisfied if (D-a)u= 0, t.e. when 

u = Ae°*, 

where A is an arbitrary constant (not necessarily real). 
By rearranging the factors so that D-j3 comes last we 
obtain another solution u = Be**. Continuing in this way, 
wo conclude that (39.2) is satisfied by 

u=Ae“*, u= Be**, u = Ce* 

individually, and is therefore also satisfied by 

u = Ac°* + Be** + Cey* + . . . 

If F(D) has n distinct linear factors D - a, . . D - k, 
we should in this manner obtain the solution 

u = A e** + . . . +Ke KX 

containing n arbitrary constants .4, . . K, and this 
expression wnll furnish the general integral of (39.2) pro¬ 
vided that its n terms are linearly independent. 

Assume for the moment that they are not linearly 
independent; that is, assume that *4, B, . . K can be 
determined so that 

Ae°* + Be** + . . . +A’c«=0 

identically. Differentiating n —1 times in succession we 
obtain the n - 1 further identities 

Aae°* + Bp e '*+ . . . +K K e'*= 0, 


Aa'-'e^ + BP"-h** + . . . +K, c«-V« = 0. 
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Eliminating Ae •*, Be", .... #e' z between the set of 


n identities, we have the condit 

1 , 1 , 

a, Q, 


a 


n—1 


P n ~\ 


on 


• } 


1 

K 


n-1 


= 0 


or (v. Aitken, Determinant's and Matrices, p. 41) 

ri(a - S) =0, 


where n(a - ft) symbolises the product of all the differences 
a _ ft a - y, \8-y, etc. As it is zero, one at least of these 
differences must be zero, contrary to the hypothesis that 
a B, .... k are all unequal. Thus e", e**, . . .,e KX 
are linearly independent, and form a fundamental set, 
provided that no two of a, ft . . k are equal. The case 
of repeated factors in F(D) will be deferred to later sections 

The form of general solution that has been obtained 
Buggests that the most direct method of dealing with an 
equation of the form F(D)u = 0 is to adopt as a tentative 

solution 

u =e rz , 


where r is a constant to be determined. 
Du=re TX = ru, D 2 u =r 2 u, . . 


Since 
D"u =r n u, 


we should have 

F{D)u = F(r)u. 

Thus u=e rx will be a particular solution if F{r) =0. Solving 
this algebraic equation (the auxiliary equation) we obtain 
as many of the exponential terms belonging to the general 
integral as it has distinct roots, but when repeated roots are 
present, the solution is incomplete. 

Example. 

y" - 2y' - 5y‘ +6i/=0. 


o 
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The auxiliary equation 

r* - 2r s - 5r + 6 = 0 

has three distinct roots, r= - 2, 1, 3. Hence the general 
integral is 

y —Ae~' z +Be x + Ce ix . 

40. Properties of the Operator F(D). In what follows 
it will be necessary to apply a polynomial operator to a 
function of the form e^X, where I is a function of x, 
at present unspecified except for the stipulation that it 
shall be differentiable at least n times. 

We have 

D . er*X = ae^X + e^DX = e“*(Z) + a)X, 

Z> 2 . e°*X = Z){e«(D + a)X} = e°*(Z) + a) 2 X, 
and in general 

D r . e°*X =e aX (Z) + a) r Z (r = 1, 2, . . n). 

Similarly, it may be proved that 

(D~a) r . e aX X=e aX D r X. 

More generally, if F{D) is any polynomial of degree n 
in D, 

F(D) . e aJ X ={D n + a 1 D n ~ l + . . . + a„_ 1 D +a n }e~X 

t =D n . e“*X +a l D n ~ l . e«*X + . . . 

+ a n _ x D . e^X +a n e as X 

•=e aX (D + a)"A' -\-a l e a3 {D + a) n-1 X + . . . 

+ a n-i enX (D + a)X +a n c ax X 
= e ax F(D + a)X .(40.1) 

At times it is necessary to operate on cos jSx or sin jS.r. 
We then have 

D cos fix= - fi sin fix, D sin fix = p cos px, 

D 2 cos px = -p- cos px, D- sin fix - fF sin fix, 

and in general 

D~ m cos px = (- p 2 )™ cos fix, D Zm sin fix = ( - p 2 )™ sin px. 
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Thus if 0{D 2 ) is any polynomial in D- (not merely in D) 

" 6 have G . D ^ cog p z = G(- fi 2 ) cos px, 

G(D 2 ) sin Px = G[~P 2 ) sin px. 

Now any polynomial containing both even and odd 
powers of D may be written 

F{ D) = G(D 2 ) + DmD*) = G(D 2 ) + B (D')D, 

and so 

/(D) coa J3x = C(D 2 ) cos jix + B(D*)D ' cos fix 

_ G( - B 2 ) cos px + //(£-)( - P sm P x ) 

= G( - P 2 ) cos px - pH ( - £ 2 ) sin px 

with analogous results for /(D) sin 0r and, more generally, 
for i’(D) sin (px +«), where € is constant. 

41 Pairs of Conjugate Factors. Suppose that among 
. r,f F(D) there occurs the conjugate imaginary 

th8 together constitute tire 

% quadratic' factor /b-a) 2 + F- Then the equat.on 
^(.DJu^O may be expressed as 

F 2 (D){(D-a) 2 + p 2 }u=0, 

where /,(D) is a polynomial of degree n-2 in D\ it is 
Batisfied if u is such that 

{(D - a) 2 + £> = 0- 

Let u o=e“w, so that (D-afe-t =e“D 2 u, and • must 
satlsfy (D 2 + j3>=0. 

The previous section shows (see also JM «>a^the 
general solution of this equation rs 

Thus the contribution of the quadraUc factor <D a) p 
to the complementary function is 

A cos Px + B sin Px). 

The same result may be obtained by the use of imaginary 
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exponents. The factor D-a+pi contributes the partial 
solution 

u 1 =e {a - fii)z =e ax e- iPz =e az {cos pz-i sin fix ); 

the factor D-a-pi contributes similarly 

M a =e a *(cos f$x+i sin fix). 

Thus \{u l -ru^=e ai cos fix and \i(u 1 -u sl )=e az am px are 
constituents of the complementary function. 

Example 1. 

y'" + 6y' + 20y = 0. 

This equation may be written as 

(jD* + 6D + 20)y =0, 

i.e. 

(.D + 2)(D* - 2D +10 )y = (D + 2){(D - 1)* + 9}y =0. 

Hence the general integral: 

y = Ae~ lx +e x {B cos 3x +C sin 3x). 

Example 2. 

y* v -y=0 

(D* - \)y = (D 1 + \ )(D - 1)(D + l)y =0 
y =A cos x + B sin x + Ce x + C'e~ x . 

Example 3. 

t/ ,v + 4m*y =0 

(D* +4 m*)y=(D* - 2mD + 2m*)(i>* +2mD + 2m*)y =0 
y =e mx (.4 cos nix + D sin nuc) +e- mx (A, cos nix + 2?, sin mx). 
This may also be written 

y =C i cosh nix cos nix + (7, cosh nix sin nix + C 3 sinh nix cos nix 

+ C 4 sinh nix sin nix. 

42. Repeated Real Factors. A repeated real factor 
in the polynomial F(D) implies that two or more of the 
numbers a, j8, . . . are equal and that the number of 
linearly distinct exponentials e OJ , c 8 *, . . . falls below n. 
We shall now indicate how this deliciency is to be supplied. 
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Let the factor D - a occur f fold, so that the equation 

F(D)n = 0 may be written 

F r {D)(D-a)'u= 0, . • • < 42 - l > 

where F,(D) is a polynomial in D of degree n-r. Thus 
a complete solution of the equation 

(D - a) r u =0, . • • ( 4 —' 

Which is What we now seek, ^satisfy the original- 

£nt 6 lt a ZZZ ^-nt u'on *, thus 




u =ve 

But (D -a) r w =e a *D T v 

A (42 21 will be satisfied if v is such that D'u= 0 ; 
that It fniay he any arbitrary polynomial of degree r - 1 

m *’ 8ay v = a+Bx+Cx*+ . ■ ■ 

and therefore the contribution of the factor (fl - a)' to 

the complementary funcUon;s ^ 

«. *— “rKht - 

ssrsrsr «■- 

plementary function are given by 

F r (D)u=0. 


Example 1 


J/'" -3/ +2t/ =0 

(/J + 2)(D - l)*!/=0 
y=e*(A +Bx)+Ce 


- 1 / 


y = 


y lv -2y' +!/=<> 

(1J- !)*(/> + l) , '/=° 

e*(yl + Bx) + e- z (/l i + # »*)• 


Example 2. 
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43. Repeated Complex Factors. Let the quadratic 
factor (D-a) 2 + fP occur r-fold; then F(D) will have the 
two repeated complex factors {D - a + j8t) r , (D -a- pi) r . 
The contribution of the former will be 

ii 1 =e aX (c os pz + i sin fix)M{x), 

and of the latter 

u 2 = e"*(cos fix -i sin px)M 1 (x), 

where M(x), M^x) are independent arbitrary polynomials 
of degree r - 1. Thus their united contribution, which is 
that of the repeated quadratic factor, will be of the form 

u = (A +Bx+ . . . +Hx r ~ 1 )e? s cos fix 

+ (^ 1 + fi 1 x + . . . +H l x r ~ 1 )e aX sin px, 

containing, as it should, 2r arbitrary constants in all. 

Example. 

y vl + m*y ,v - m 4 y' - m*y =0 
(D* - m*)(D' +m*)*y =0 

y = (A + fix) cos nix +(.4! +fi,x) sin mx + Ce mx + C,c~ mz . 

44. Inverse Operators. In preparation for the dis¬ 
cussion of a particular integral of the non-reduced linear 
equation, we shall establish certain properties of the 
operator inverse to F(D). 

The operators d()/dx and /()dx are inverses of one 
another; the former has been denoted by D, and it is 

natural to denote the latter by D~ l or This is justified 

‘■y * 

D~ l {Dy) = D[D~ l y) -y, 


y 

* Noto that 2 ^ must not be written, as it is doubtful whether it 
1 1 , , 

means ^ y or y —, and the operator is understood to have effect only 
on that which follows it. 
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indicating that the processes of integration and diilerentia- 
tion, successively performed, cancel one another. This is 
analogous to the relation between functions and their 
inverses; e.g. if exp denotes the exponential and log the 
natural logarithm, 

exp log x = e log x =x, 

a relation which holds identically, not for any particular .r, 
and which justifies the symbolic relation 

exp log = 1. 

But a word of caution must be added. Since differentia¬ 
tion ( D) is unambiguous, but integration (-D -1 ) is ambiguous, 
it may happen that if the integration is performed after the 
differentiation, an ambiguous element may remain, tor 
example we have 

D{D~H2x + 1)} = D(x 2 + x) = 2x +1, 

justifying the statement D.D~ l = l. But, on the other 
hand, 

D~'{D(2x +1)} = Z>-»(2) = 2x + C, 

where C is an arbitrary constant, so that the statement 
l)-i D = i i 8 subject to the qualification that the constant 
which D _1 introduces has been properly chosen. 

Just as D r indicates the operation of r-ple differentiation, 
so D~ r will indicate r-fold integration. This is a process 
which has a clear meaning, and the relations D r (D ') = 1 
and D~ r (D r ) = 1 are true, the former invariably, the latter 
subject to the proper choice of the coefficients in the 
polynomial of degree r -1 introduced by the integration. 
But the operator inverse to the polynomial differential 
operator F{D) has not yet been identified with any manip¬ 
ulative process, and can be interpreted only in special 
cases which will now be considered. 

We recall the relation (40.1) 

F(D) . e°*X «= e“F(D + a) X. 
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Now F(D + a)X is a function of r, say Y, dependent upon 
X, t'.c. 


F(D) . e^X = e** Y, where }' = F(D + a)X. 

Thus e°*Y is obtained from e^X by the operator F(D)\ 
conversely e^X may be obtained by applying the inverse 
operator on e^Y, or symbolically 


By the same argument, since F(D + a)X = Y, we have 

x _L_y 

F(D+a) ’ 

and so we have the important relation 


F(D) 


{e* 1 Y) = e aI 


1 


D(D + a) 


. (44.1) 


and in particular 



(44.2) 


In the case X = 1 
= F(a)e ax that 


we deduce from the relation F{D)e ax 



. (44.3) 


provided that F(a) ¥= 0. This interprets the inverse operator 
F~ i (D) with respect to the exponential function. The 
exceptional case F(a)= 0 implies that D- a is a factor of 
F(D)\ let it be a factor repeated r times, and write 

F(D)=G{D){D-aY, where (7(a) *0. 

Let F(D) be compounded of any two polynomial 
operators (7(D), H[D) and let 


Y = F(D)X = H(D)G(D)X. 
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Then 




1 1 V 

G(D)II(D) F(D) 


Thus when F(D) =G{D){D - a) r t 

1 1 


_ 


| 1 e°* 

j = ( 7 T^r W 


^(2))° (D-a) r \<?(/)) . 

Putting 7-1 in (44.1) and noting that by ordinary in 
tegration 


D '-*• 


, *. 1 r> i 1= *: 

• 1 = 2 • Ip'“0. D '' ~ H ’ 


ignoring the arbitrary constant that enters at each step, 
we find that 


1 


1 


e. ax = 


e ax = 


e aX x r 

rlfl(7)‘ 


. (44.4) 


T(Df (D-a)'G(D) 

Thus we have shown that 

(i) a particular integral of F(D)„ is y-**IH*) 

provided F(o) #0. 

(U) when F{D) is of the form G(D)(D-a) r . whore 

0(.) Jo, or when the auxiliary equation has «aa an r-ple 
root, the simplest particular integral is e x /{ )}■ 


Example 1. » 

/+y'-6«/=8e«*. 

Auxiliary equation : 

r j +r _G=0, roots r =2, —3. 

Complementary function : 

u=Ae %m + Br* m . 
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This does not contain a term in e 3x , so that a particular integral 
is 

y=(D'+D-6)-'(Se 3 *) 

= 8e sx /(3* +3 -6)=$e 8X . 

General Integral: 

y=Ae tx + Be~ 3x +$e 3x . 

Note. When a particular integral satisfying certain speci¬ 
fied conditions is required, it may be deduced from the general 
integral by assigning the proper values to the constants 
introduced by the complementary function. The number of 
independent conditions imposed cannot exceed the number of 
disposable constants which is given by the order of the equation. 

For example, to determine the particular solution in Ex. 1 
such that y = 0, y' = 3 when x = 0. We have 

y = Ae 3x +Be~ 3x + £ e 3X , y' = 2 Ae 3x - 3 Be~ tx + 4e 3X ; 

hence 

0 =A + B + $, 3=24-35+4. 

Solving, we have A = - 1, B = - }, and the particular integral 
satisfying the condition imposed is 

y = ^e 3x -e tx -%e- 3X . 

Example 2. Solution of 

y'-3y'+2y=e x +e 3x 
such that y=y' = \ when x = 0. 

Complementary function u = Ae x + Be 3x . This contains both 
terms of the right-hand side of the equation, i.e. F(a) =0 in both 

cases. 

For a particular integral 

y =(D - 2)- , (-D - + e ix ). 

Taking the exponentials separately, substitute the coefficient 
of x for D in the non-critical factor only, thus: 

y= -(D- l) -, e x + (D - 2)~ , e ,z 
= - e z D ~ l . 1 + c 3x D~ 1 . 1 = - xe x + xe 3x . 

Hence the general integral 

y = (4 - x)e x + (B + x)e tx . 



linear equations 

To obtain the particular integral, make the required aub 

atitutions in this and in 

y' = (A - 1 - x)e x + (2 B + 1 + 2x)e tz , 

obtaining 

\=A+B, l=A+2B, whence A-1, B-V. 

The required particular integral therefore is 

y =(1 -x)e x +xe tx . 


45 Inverse Operators Relative to a Periodic Function. 
Equations of the form F(D)y = sin (mx + e) are of sufficient 
importance in applied mathematics to deservei spec 
mention. We separate F(D) into its even and odd parts, 
F(D)=0(D*) + DH{D*). Then since 

F(D) sin (mx + 0 = {<?(»*) + DH(D*)} sin (mx + e) 

= {G( - rn 2 ) + DI1( - m 2 )} sin (mx + 0 

= (g+hD) sin ( mx+e ), 

where g and h are constants, we conclude that, if ^ 0, 

h + 0 , x 

sin ( 1 ^+ 0 =^ sin (mx + e). 

To evaluate this expression, we proceed as follows 
, g-hD . . 

I % v /tnr i. r 1 


^sin^+O-^siu^+e) 


g - hD . , 

l , n 8in (w^+ € ) 

4- h*rnr 


g sin (mx + e) - h m cos (mx+e) 

g 2 + hhri 1 

cos tl sin ( mx+Q -s in t i cos (mx + Q 

y/tf+Wm 2 ) 

Bin(mx + e-eO wher0 tan 

1 'Tf^ 2 ™ 2 ) * 
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Alternatively sines and cosines may be replaced by 
imaginary exponentials. Since e iB = cos d + i sin 9 we write 
cos 9 = Re <fl , sin 9 = Ie ie , where R denotes the real part, and 
I the coefficient of t in the imaginary part of e' 9 . The 
above equation can therefore be written F(D)y =Ie t(rnz+f) . 
We shall, however, drop the symbol I for the time being 
and consider 

F(D)y = e« m,+c >; . . . (45.1) 

the coefficients in F(D) are still supposed to be real, but y 
will now be an imaginary function of the real variable x. 
Following the methods of § 44 we find that 

l e t'(mx+0 

y 


F(D) 


gi(mr+ «) 


F(im) 


But since 


F(im) = G( - m 2 ) + imH( - in 2 ) —g+ imh, 


we have 


gt(mx+«) g-imh 
_ __ _ y_ e i(mx + t) 

g + imh g 2 + m 2 h 2 


or 


(g 2 + m 2 h 2 )y = (g - imh){cos (mx + e) + 1 sin ( mx + «)} 

= g cos [mx + c) + mh sin (mx + e) 

+ i{g sin (mx+e) - mh cos (mx+c)}. 

This furnishes a particular integral of (45.1); taking in 
turn the imaginary and the real components we see that 

psin (mx+e) - mh cos (mx + e) sin (mx + e -€ x ) 

V= g 2 + m 2 h 2 = y/(g 2 + m 2 h 2 ) 

is a particular integral of F(D)y = sin (mx+e), and 



g cos (mx +€) + mh sin (mx + e) 
g 2 + m 2 + h 2 


cos (mx +e — c x ) 
V (9 2 + m*h 2 ) 


is a particular integral of F(D)y= cos (mx+e). 
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Example 1. Differential Ejuatum for a «"&*£** 
Circuit. If i (amperes) is the current intensity * at Ume ! 
(seconds) in a circuit containing inductance L (henr.es), 
resistance R (ohms) and capacity C (farads), under at 
pressed electromotive force e (volts), we have the differential 

equation 

dH di i de 
L di' +R dt + C~dt' 

A particular solution of the reduced equation is c r ‘. where 
and three cases are possible: 

(i) If R'C > 4 L, the radical is real but numerically less than 
unity; hence the complementary function is of the form 
At-*' +Be~ bl , where a and b are real and positive. 

(ii) If R'C=iL. the complementary function has the form 
t~ at (A +Bt), where a = R/2L and is positive. 

(iii) If R'C > 4 L the radical is imaginary! the complementary 
function takes the form .-'M cos 6, + B sin bt). where o » 

positive. 

In all tlireo cases, the exponential rapidly becomes small 
in fact negligible in practice, so that the complementary 

function represents a transient effect. . ., . 

The most important form of . consists of a sum ofj»^ 
terms; os their effect is additive wo confine our attention to 

a single term 

e =E sin {ait + a) 

when a particular integral of the differential equation is given 
by 

;_ CD -— E sin (tut + a)' 

CLD 1 + CRD + 1 

-a* 

mathematical symbols i and t . 
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Evaluating this expression by the above method, we find 
•_ CEo) 

~ \/{(l -CLw*)*+(CR<o)*} Sm {0)t + ° + 1” ~ A’ 

where tan p=CRw/( 1 - CLto *). 

This gives the ateady-state or permanent current, whose 
phase is behind that of the e.m.f. by the angle 90 3 - (}. 

Example 2. 

y ' +m s y -a cos mx +6 sin mx, 

where a and 6 are given constants. This is a case where the 

lght-hand member is contained in the complementary function, 
if we take the equation 

y’ +m t y=e imx , 

wo have the particular integral (c/. 44 . 4 ) 

1 e imx 


y= --- g<mx 

D>+m> e ~ 2im 
x . 

=«-(*in mx - i cos mx). 


2 m 

totegrflTof ira " 8inary PaFta “ tUm W6 308 that P articular 

y' + m'y = cos mx, t/' + m *y =sin mx 

are respectively 

y = ' 27 n mmX> y= cos mx, 
and a particular integral of the proposed equation is 

X 

y ~2m‘ a Sm uuc ~b cos mx )• 

46 Development of an Inverse Operator. Let the 
poljnonnal operator F(D), with constant coefficients 
ntain the factors D - a repeated r times, D -\3 repeated 
* times and so on. Then we may write P 

__ ^1 a >2 a £ 

r(-D)" D~~a + Id - „)«' + • • • + (0^p + flT2 + • • (46.1) 
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for if we operate on both members by F(D) the left-hand 
becomes 1; the right-hand becomes a polynomial in D 
which will reduce to 1 if the constants A v A z , . . 
A r) B x . . . are determined just as they would be if the 
right-hand member were the algebraic development of the 

left-hand in partial fractions. 

We may now interpret (D - a) -1 . x *, where k is a 



for if we operate on the right-hand member by D a we 
recover x\ In similar fashion we can interpret (D - a)~ r 


1 ( D v-» 1 

-J 1 — *' = . 


a*. a / a 2 


27) 3Z) 2 


1 + — + 


a* 


and expanding each term of (46.1) as a series of ascending 
powers of D, we arrive at the following interpretation of 

F~'(D)x *: 

F~ 1 {D)x t = (c 0 + c x D + c 2 D 2 + . . . +C'D')x*, (4G.2) 

where the series on the right contains the first k + 1 terms of 
the algebraic development of 1 IF(D) in ascending powers of 
D. It terminates in D' since D'x* =0 for r > k. It 
follows that if P(X) is a polynomial of degree *, 

F-'{D) .P{x) = {c 0 +c l D+c 2 D 2 + . . . +c k D*)P(x). (40.3) 

Example. Part : cular integral of 

+ 12 y=e ,jr (x*-5x*) 

y =(3 -D)- , (4 -D)- x . e«*(x» - 5x*) 

= e* z (l - D)~'( x * -6z*) 

=xJn\ + D+D'+D a ){l+hD+i D2 I J D, M X ' 5x ‘> 

D + lD t +\b D3 )(x 3 -W) 
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47. General Integral by Quadratures. When the 
factors of F(D) are distinct, the equation F(D)y=f(x) 
may be integrated completely by quadratures. For we 
have 


1 


F(D) 

Now 


• /(*) = 


1 


1 


F'{a)[D -a) F'(P){D - (3) 


1 •/ (x) = n^- € °*{ e ' aI / (a;)} 


D — a 




= £ “{f 


e~ aX f(x)dx + C 



The general integral therefore is 

/”(aj{j e ’ “feW* + C, i| + + • • • ( 47 • 1 ) 

or it may be expressed in the form 


■ ■ " (47 ' 2) 

where the lower terminals x v x 2 , . . . are arbitrary, or 
again as 

n ga(Z-l) € gj*-q 1 

~F 7 {a) + 'F T {p) + ‘ * *|/^» ( 47 * 3 ) 

where u (x) is the complementary function, and the lower 
terminal of integration (which would merely lead to terms 
already existing in w(x)) is ignored. 

If D - a is a double factor, let F(D) = (D - a) 2 G(D). 
Then 


l 


•/(*) 


J. 


1 


G'(a) 


F(D) • J [0(a)(1) - a) 2 {G(a)} 2 (Z> - a 


Now 



1 , 1 

■ f(x) = e “^( e ' “/(*)! = «“{//«- “Mi* + C lX + CJ. 
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Thus a double factor leads to a double integral and in 
general an r-fold factor to an r-fold integral. 

Example. Particular integral of 

y’ +a ty=e“ x CO? ax 

1 1 f 1 _1—Ve 1 ** cos ax 

V ~ C ° 3 01 ~ 2ai\ D - ai D+ai\ 

_ }__[*wid - 1 ) .g-aitx-OjeM 1 cos at dt 

~2ai) 

= Ij*sin o(x-«K* cos atdt 
__L fV'fsin ax +sin (ax - 2 at))dt 

2aJ 

, 1 


= 2aif ‘ X 8m ° X + 2a(/i , + 4a *) 

i f 2o . \ 

_ 1 — cos ax +— s* n ax b 

-»*+4a* l ^ ' 


c' ,x l2a cos ax - /» sin ax) 

t\ 


48 The Euler Linear Equation. 

This equation may be reduced to on* -i* — 
coefficients by the substitution *-«• 

dy dy d< 1* ti.lti.i-A 

H-dt-di-xdf dx* **dt‘ x-dt 

and in general, if D now represents dO/dl, 

J?1 = D(D- 1 ) • • • (D-r + \)y. 

dx T . 

and the transformed equation can * “"^ e ‘ (48 2) 

^ + A>Dn 7- ' U A a can ife expressed 

whose constant coefiicients A v . • •» n 

linearly in terms of a v . . -, a n- n 
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Alternatively we may introduce a new operator {h 
dy d 2 y 

*£-%>> . 

ir 0 =w-1) • • • ( & - r+1 te 

so that the transformed equation is 

(& n +A 1 &"-i+ . . . +A n _ l & + A n )y=f(x) (48.3) 

whose coefficients are identical with those of (48.2). 

Taking as basis the fact that the operators # - a and 
$ - P are permutable when a, )3 are constants, we may 
establish a calculus of the polynomial operator F(&) and of 
its inverse F~ l (&) similar to that of F(D) and its inverse. 
In particular 

F(&)xr = F(r)x\ F(&)x'X = x'F(& + r)X, 
where X is a function of x, and conversely, 

F-'i&jx^x'/Fir) provided F(r) ^ 0, . (48.4) 

F-'{&)xrX=x'F- 1 (&+r)X. . . (48.5) 

If F(&) = (&-r)°G(&), 

then 

F~ 1 (&)x r = (# - r)~ a G~ 1 (&)x r 

= (& -r)~ a x r JG(r) =x r &~ a (\)/G(r). 

But since 

&>J = xdy/dx, ft-iy = jydx/x. 

Hence 

&~ l (l)=fdxfx=\og | x | 

& ~ 2 d) =JIog | x | dx/x = J(log | x |) 2 , 
and in general 

#"*(1) = (log | x |) a /a ! 

The complementary function will contain terms like 
Cx T , where r is a root of the indicial equation 

r»+4 1 r , -» + ■ . . 
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If this equation has n distinct roots ’-a, ft • • - * the 
complete complementary function will be 

C l X a +C»X li + . • • + C n X *- 

equation is that in which /(*) » s cci I obtained 

undetermined coefficients (see Ex. 2 below). 

An equation of the more general ty pe ^ 

lnay be brought into the above 

srr %£• % 4- is 

ax + 6 -e'. The ind.c.al equatmn fou'm oy 
(ax +b)’ in the reduced equation and proceeding 

Example 1 . , , 

x l y" - 2xi/' + 2>j —x +X 1 log I x I +x * 

In terms of the operator * the left-hand mombor a, 

xr of caf tor,m 0,1 th0 

Of the equation mag be foun « • = , og , ,, 

( 0 -\)-i (0 - 2 )- i x=xO-H&-') - 1 

<0 - !,-■«> - 2) :^. log i x, =**u ~ w-' ■ [o-i * 1 

=J( 0 -' - 1 ) log I * I = z,{ i (log I X ° e * 

(&- 1 )-»(d -2)- , x» =(2)- , (l)' ,r * = l x • 
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The general integral therefore is 

y=Ax + Bx* -x log | x | +x l {$(log | x |)* - log | x |} +%x*. 
Example 2. 

(2x +1 )Y + (4* + 2 )y' - 4 y =x *. 

Substituting (2x + l) r in the reduced equation, we get 

4r(r - 1)+4r - 4 =0 or r* -1=0. 

Hence the complementary function 

A[2x +1) + B(2x + l) -1 . 

As x* = \(2x + 1)* - $(2x + 1) + I and 2x + l occurs in the 
complementary function, a particular integral of the form 

y =a(2x + 1)* +6(2x + 1) log | 2x + 1 | +c 

may bo assumed. Substituting we find that 

12a(2x + 1)* + 86(2x + 1) - 4c =*», 

whence 

a= l's» b= ~l\’ c= - tV 

49. The Laplace Linear Equation. In this equation 
the coellicients are linear functions of the independent 
variable, thus: 

K + V)y (n) + ( a i+^)'/ (n - 1) + . • . +{a n +b n x)y = 0. (49.1) 

If the constants b were all zero, the solution would be of 
exponential form. It was probably this consideration 
that suggested the possibility of satisfying the equation 
by an integral of the form 

I =je rl f(t)dt . . . (49.2) 

in which /(/) and a range of integration independent of x 
are at our disposal. We have in general 

d r l/dx r =fe x, t r f{t)dt, 

so that if Ly represents the left-hand member of (49.1), 
we have 

LI = JV'{K + b 0 x)l n + (a l + &JX)*"-! + 

=Ie*‘{P(t)+xQ(t))f(t)dt, 


• • • 


+ «n +Kx}f(t)dt 



LINEAR EQUATIONS 


105 


where 


P(f) =<*</" W" 1 + • • * +r '"’ 

Q(t)=b 0 t n +b l t n ~ 1 + • • • + h n- 


Now 


=\Qm4r ,]it=[eZ ' Qm)] " \ cr 'dt mm]d ‘- 

Thus LI =0 if 

(i) f(t) is such that 

(U) /(() having been determined, the range of integration 
is such that 

[e"<2(0/(0] -0 

independently of x. The former eondition is satisfied by 


1 

/«)=w)“ p \J5w 

and the latter becomes 

j-grf + = 0. 


■} 


. (49.3) 


. (49.4) 


Note. It is often advantageous to replace I in (49.2) by a 
jonstant multiple as - 1 or it. 


Example. 
In this case 


xy'+(2n + l )!/'+*!/=°- 


P(0=(2n + 1)*, Q(0 =l* + *• 

f(PIQ)dt=\og (i* + i) n+1 ; +l) 

Then, replacing t by it, 

J=JV*‘( 1 -«■)"•** 

satisfies the differentia, elation provided that , is convergent 
and that the terminals are so chosen that 

[e ixt { 1 -<*) n+ *]=° 
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identically. The former condition is satisfied if n is positive; 
the latter if the range of integration is ( -1, +1). The 
imaginary part of the integral vanishes; taking the real part 
we conclude that the equation is satisfied by 

y = J cos (xt)( 1 - t*) n -idt. 

Writing t =cos 0 we find that an equivalent solution is 

f ff 

y = 1 cos (x cos 0) sin ,n 0 dO, 

Jo 

an integral well known in the theory of Bessel functions. 


50. Variation of Parameters. When the complementary 
function of any linear equation 

t/<"> +p 1 i/ n ~ 1 > + . . . +p n - 1 y'+p n y =f(x) (50.1) 

is known, the complete solution may be obtained by 
quadratures as follows. If the complementary function is 

u=C 1 u i+ . . . +C n u n 

we replace the constants C by symbols V representing 
functions of x, thus 

y-v l u l + ... + v n u n . 

These n symbols are at our disposal, so we proceed to set up 
n linear relations between the derivatives V x \ . . ., F n '. 
First of all we have 


provided that 


Then 


y = *>i' + . . . + v n u n ’ 

Fj'hj + . . . + V n 'u n =0. . .(60.2) 


y' = I>r+ . . . + V n u n ” 

provided that 

FjV+...+F.V-0, . 

and we continue up to the stage 

y ( «-i) = F 1 ti l ‘«-»»+ . . . +F f ,V n - 1) 


. (50.3) 
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provided that ^ 

V l 'u l in ~ 2) + ■ • • + y n ' u n n ~ 2) =0 - ' (o04) 

Thus n -1 linear homogeneous relations have been estab- 

r between V ' V If we substitute the above 

hshed between V x . • ,, ^ als0 substitute 

expressions for y, y, • • •>!/ lu w v * h 
for t/<"> its value 

!>,<■' + . . • + V „«„<*> + W 1 ’ + • • • + 

we shall find that 

+ • • • + W'” =/{x) - • (0 } 

This is the final relation between 1Y> • • •• 1 « ; w ® 7 So1 ^® 
the system algebraically and obtain V i> • • •> n > 
quadratures, thus obtaining a complete solution 

equation^ ^ ^ the equation of the second order ^ 

y'+py'+qy= r > ■ • • (50- ) 

if Wl , « 2 are the constituents of the complementary function 

we write y = V x u x + V 2 u 2 and determine F,, V t by the 
relations 

+ V 2 u 2 = f ) > Ki'u/ + V 2 u 2 =r. 

If W-u x u 2 -u x u t we have 

\W X ' =-ru 2 , WV 2 '=ru lt 

so that the general integral is 

y _ -«,{>+ • • (50 - 7) 


Example. 

y r +y =sec x. 

With the above notation we have 

y z= Vi cos x + V, sin x, 


W = 1 


y « - COB xj sec X sin x dx + sin xj sec x cos x dx 
= (A + log | cos* |) cos x + (B+*)«•“ *• 


Henoe 
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51. Linear Systems with Constant Coefficients. By 
introducing a new variable t upon which x and y are 
dependent, the linear fractional equation (17.3) may be 
replaced by a system of two linear equations 


dx dy 

T. =lx + my, —=ax+by. 


. (51.1) 


This system will be satisfied by x = Ae rt , y=>A x e rt if 


(l -r)A +mA l =0, aA +(b-r)A x =0. 
Eliminating A, *4, we have the characteristic equation 
(l-r)(b -r) -am =0 or r 2 -(6 +l)r +bl-am =0. 


If this equation has two distinct roots (real or complex) 
a and /?, the general solution of the system is 

x = Ae aI + Be*', my = (a - l)A<t* + (0 - l)Bef*, (51.2) 

dependent therefore upon two arbitrary constants A, B. 

When r = a is a double root, we take, by analogy with 
the linear equation with constant coefficients, a tentative 
solution 


x = e aJ (At + B), y = e n, (A x t + B X ). 

On substitution wo find that 

(l -a)A+7nA x =Q, aA + (6 - ct)^ =0, 

(l -a)B +tnB x =A, aB + (b - a)B 1 =A 1 , 

whence 

x=e°'(At + B), my=e al {(a-l)At + (a-l)B+A}. . (61.3) 

We may regard (51.2) or (51.3), as the case may be, as 
a parametric general integral of the linear fractional 
equation. 

An alternative method is to eliminate all of the de¬ 
pendent variables but one, which will be illustrated by the 
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system of three homogeneous equations 

l 2y + L * z=0 1 ... .. 

MiX + Mrf + M 3 2 = 0 |* * ' ( L4) 

N x x+ N 2 y+ N^= O' 

where the 9 coefficients L l to N , represent polynomiets 
in the operator d()/di with constant eoefficen s. It wiU 
be supposed that the system ,3 determinate * 1 Inch nnpi es 

that the determinant A of the eoeffic.ents is not r^e o 

This determinant itself will, in general, be a polynomial 
in dQldt) its degree is the order of the syst<*'»- 

By operating on the left-hand members of (oM) 

by the minors of/,,, -V„ irrespectively, ^ obtain Ax-0, 

Similarly, Ay-0, A«-0. Thus *, y, » « 
integrals of the reduced linear equation °> 

order is that of the system, provided alway * J^ion 

a mere constant, in which case the only possible solution 

would bo x = y =z = 0. f 

Setting this case aside, *, y, z will be of the lorms 

x = + Q X *F + • • •> 

y = P 2 e* + Q 2 e?‘ + • . 
z = + .... 

where a 8 . . . are the distinct roots of the characteristic 

similarly for the rest. By subsUtution in (.>h 4 ) we M a^^.^ 
relationships between each group o > limited 

whereby their coefficients bcconm^enUip^ ^ 

number from among them. It P . • el|Ua i to 

number of constants finally remaining arb.trary is cqua 

the order of the system. 

By exactly the same process, a system 

L x x + M x y + tf| 2 =/i( 0 l (515) 

L £ x + M# N.,z =/,«) • • ^ l o) 

L 3 x + M jD + A 7 3 2 "/a(0J 


(51.5) 
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is transformed into 

Ay *=g 2 {t), Az=g z (t). . (51.6) 

If x = X, y=Y, z = Z is the complete solution of the 
corresponding reduced system (51.4) and if a particular 
solution of (51.5) is x = <ft l , y = <f> 2 , z = <f> 3 , the general solution 
of the system will be 

x = X +<f> lt y=Y + (f> 2> z=>Z + <f> 3 . 


Example 1. 





The characteristic equation of the system is 


I 1-r -2 =0 
j 5 3-r 

or r J - 4r + 13 =0, whence r = 2 ± 3i. Substituting 

x=e il (A cos 3 1 +B sin 3 1), y=e :t (A l cos 3 t + B x sin 3<) 

in the system, wo obtain the four consistent equations 


A+3B=-2A X , 3A-B-2B X , 

A X -3B X = -5A, ZA X +B X = - 5B, 

and thus the general solution is 

x =e ,, (.4 cos 3t +B sin 3t), 
y = $c : '{ - (A + 3 B) cos 3 1 + ( 3A - B) sin 3 1). 


Example 2. 

I />(D + 1 )x + D(D' + D - 1 )y + (D - 1 )z =e* 

•. Dx + (D J -l)y -z=0 
l x +Dy +Dz=0 

where D—d()/dt. By elimination, we obtain 

(D - 1 ) ! x = -c‘, (D - l)*y = 2c t , (D - 1)«2= -e\ 

whence 

I=e V'+B-n r =e‘(^»< + 5, - 
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On substituting for x. y. z «e obtain the six consistent 
equations 

2 A + Ax= -4, 3^1 +2B+1A 1 +B l +A t = -6, 

a _ a = - 3 A + B + 2/41 — B 2 = - 2, 
i4+i4l+4,= -V B4-4 1 +B I +4.+B 1 =0, 

giving 

y4j = - 2>1 -4, x4, =/l + 3, 

Bj = 4.4 -2B + 7, B,= -3.4+B-6. 

Hence the general integral is 

x=e'(At +B - it*), y=c*{- W A + 4)/ + 4 A - 2B + 7 + 1'\, 

z=e t {(A+S)t-3A + B-6 

The system is thus of the second order. 
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52. Solution Developed as a Taylor Series. We 
consider a linear equation of the second order 

d 2 y dy 

te* = • • • (52.1) 

in which we shall regard z as a complex variable. 

This type includes many equations of very great importance 
which cannot be solved in terms of simple combinations of 
elementary functions. Given such an equation, the usual 
procedure is to express the solution (which may be a 
solution satisfying certain initial conditions) in the form 
of an infinite series from which tables of the value of the 
solution may, if desired, be computed. Thus the con¬ 
vergence of the series is important, not only as a basis of 
the validity of the process, but also as an indication of the 
practica value of the result, for a slowly converging series 
is of little use to a computer. 

It will bo assumed that the coefficients p(x), q(x) are 

one-valued, and have derivatives of all orders except 

possibly for certain isolated values of s. Let x = a be a 

va ue for which p, q and all derivatives are finite: we shall 

obtain a solution such that y and y' have assigned finite 
values y 0 , y Q when x=a. 6 

By substituting in (52.1) we obtain the corresponding 
value of the second derivative, namely 8 

y<)' =p(a)y 0 ' +q[a)y 0 . 
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Differentiating (52.1), we obtain 

y"' = p(x)y' + {p'( x ) + qWW + 

and as y 0 , y 0 ', </„' are known, y 0 "' can be obtained immedi- 
ately. Continuing the process, we obtain the values ot 
successive derivatives for x=a, and thus we have the 
coefficients in the Taylor series 

lx - a ) 2 (x-aV» 

y=y 0 +yo'( x - a ) + yo- 2 !~ + * * * +!/o n! 

Borrowing from the language of the theory of functions 
of a complex variable, a point x=a at which the above 
conditions are satisfied is said to be an ordinary point ot the 

equation. 


Example 1. To find the solution of 

y* +xy=0 

such that y = A, y'=B when x=0. 

Wo have 

y" = -xy. y"=-zy'-y, y'-= - xy-- iy'. 

and in general 

y<")= - X y [n - l) -(n -2)y' n -". 

The origin is an ordinary point; putting x=0 in the nbo\ 
we have 

I/o' =0, y," = -!/o. !/o lr = - 2 Vo'. 'f = ' 3 ^' = °’ 


and in genorul 

y 0 ,n) = 

We thus find that 

„,<•»> =(-l)(-4)<-7) . . • <- 3,1 +2)A \ 
y 0 «»"+») =( -2)( -5)( -8) . . . (-3f» + l)£. !/» 


Thus the required solution is 


x 


V " i “ , -»i + OI 


» 4c* 28x» 


9! 


+ 



2-r 4 10x 7 

- 77 + 7TT 

4 ! / ! 


80.r 10 

— + . • 
10! 



and it may bo verified that (i) the ratio test proves the con- 
vergence of each series for all values of (ii) oac 1 8on 
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separately satisfies the differential equation. Since A and B 
may be regarded as arbitrary constants, the above is an 
expression of the general solution. 

Note 1. When p{x), q(x) are polynomials of low degree, a 
more practical way of obtaining the solution is to assume a 
series of ascending powers of x 

y =c 0 +c x x -t-c*®* -F. . .=Ec r z r , 

to substitute this series in the differential equation and to 
determine the coefficients so that the equation is identically 
satisfied; that is to say, so that the coefficient of each individual 
power of x cancels out. 


Example 2. 

y ' -xy' +ny = 0. 

Substituting the above series in the left-hand member and 
equating to zero the coefficients of successive powers of x in the 
result, we obtain the following set of recurrence relatione :— 

2 c a +?ic 0 = 0 , 603 +(n - l)c a =0, 

12c 4 +(n -2)Cj=0, 20c s +(n -3)c s =0, 

and in general, by equating to zero the coefficient of x r , 

(r + l)(r +2)c r+# +(n -r)c r =0 (r=0, 1, 2, . . .). 

We see that the coefficients of even rank (t.e. suffix) depend on 
c 0 but not on c,; those of odd rank on c x but not on c 0 , and 
that c 0 , c, are arbitrary. Denoting by S x (x) the series obtained 
by taking c„ = 1, c x =0 and by S t (x) that obtained by assuming 
c 0 =0, Cj = 1, we have 


SM = l - 2 V - ^ + ■ 

£ l 4 ! O ! 


_ n-l (n-lKn-3) _ (n - I)(n-3)(„-5)^, 

3 1 6! 7 


• • 


Since 



r - n 

(r + 1 )(r +2) 


-> 0 


as r -*■ co 


for any finite value of n, both series convergo for all values of x. 
The general integral is y =AS x {x) +BS t (x), where A, B are 
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arbitrary constant. It wi.l be noted tba, ^ ™ 

* * • PO»— 

° f % S 2T It U sometimes 

rr rtfprnv” tudied. The transformation most 

frequently employed is of the form 

y =uv 9 

where y is the origins., . the new dependent variab.e, and u 

a chosen function of x. 

For example, if the equation 

eoefflcienfta „ (r+ 2 ) c^ +(n + J)c . 1Cr _, =0 . 

This is a three-term recurrence 

difficult to manipulate than one that involve , 

If, however, we make the substitution 

y=e-\*'v. 

BO that ^ =C -J*V -XV' + (*** - tM* 

we find that the equation in u is 

v * -x» +nw = 0, 

which is precisely the equation considered above. 

i When we consider an equa- 

63. Regular Singularity. 

.. . -. »..•;«“ r E“ : 

of both Q(x) and lt(x), , 0 cedure 0 f § 52 is inapplicable. 

value of y'for x = a, an P iw ^ point or sirujuUirity of 

We then say that x - a is a ,/ 
the equation. 
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The form of a solution valid in the neighbourhood of a 
singularity may therefore differ from that of a solution 
appropriate to an ordinary point. It is important to 
investigate this difference, which is not merely a matter of 
theoretical interest, but of practical importance, as the 
solutions valid in the neighbourhood of singular points are 
as a general rule those required in physical applications. 

We find a clue to a fruitful line of approach in the 

equation 

(x - a) Y + b{x - a)y' + cij =0, 

where b and c are constants. This equation has (§ 48) a 
solution of the form j/ = (x-a)', so that the singularity of 
the equation may imply a singularity (e.g. a pole or a 
branch point) in its solution. When (53.1) can be written 

in the form 

(x - ay-P^x),/ + (x -«)$,(%' + Ri( x )y =o. 

such that Pj(a)^0 and P^x), Q x {x), R x (x) are finite in the 
neighbourhood of x = a, it is thus natural to assume a 
solution of the form 

y = (x -a)'F(x), 

where F{x) is a function which can be developed as a Taylor 
series in x-a. It will be found on substitution, equating 
to zero the lowest power of x-a, that a is determined by 
the indicial equation of the second degree 

ct(ct - l)Pi(u) + oQ x (a) + R x {a) =0, 

giving in general two distinct values of a and leading to two 
distinct solutions. Such a singularity is termed regular. 

54. The Hypergeometric Equation. This is the 
equation 

x(l - *)!l' +{y-{a + P + 1 )*)!/' - a-P'J =0, (54.1) 

in which a, />, y are constants, at present unrestricted. On 
multiplying the equation throughout by x, it will be seen 
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t _ _o is a regular singularity; similarly, the singular,ty 
* = 1 is'regular a^d thereto no other finite -ngulanUes. 
Confining our attention to the smgulanty ,-0, we 

assume the series 

< V f+c,z' +I + ■ ■ ■ +V" + • • • (c ” #0) ’ ( 

in which the index o is for the moment regardedl -» 
arbitrary parameter, and denote^ *1*.° ^ o » 

aswra: »—-—■ 

ing powers of x, that 
Ly = a{a -1 + yYa*? 

+ £ {(a +r + l)(o +r + yKrH - (" + r + a >< a + r + P )C ' ]X ’* r - 

„ now, whatever o may be, we choose the coefticients 
r r .so that 

;:;;;i„o + r + y ,c rtl = (o + t + a„o + r + ^(r>0,( 5 4,, 

leaving c 0 arbitrary, we find that 


, “ +fl , ( ff + a)(o+a 1 lKo 1 Wo ± £ ± l), 

(cr + a)(ar + p)^ __ — —r^/TTZuVr X v + Y) 


, nnd these coefficients are finite provided that 

:no a :u + T“ U a negative integer. Sett.ng such 

cases aside, wc have proved that it 
</>(*> = 

f «,(a+a). . .(a+a_+ 

Y + X'(<t + 1) • • • ( CT + 

and if the series is convergent, then 

r4(X, a)=o(a - 1 +yYo x " 1 - 

The question of convergence cun be settled once and for al, 
by the ratio test, for I 


CqXT 


r)(a +y) 7TT(o+y+r-\) J 
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<Wi (g+r+a)(a+r + ffl 
c r (<7+r + l)(a+r + y) 

a, a, p, y being finite, and so the series converges for 
1*1 <*• 

We can make L<j>(x, ct) =0, and thus satisfy the equation, 
by taking 

a(a-l+y)=0 . . . (54.4) 


provided that the stipulation that neither a nor a + y -1 is 
a negative integer is not violated. Thus the roots of the 
indicial equation are a = 0 provided y -1 is not a negative 
integer, and 1 - y provided 1 - y is not a negative integer. 

Take first the index a =0. A solution is provided by 


<f>(x, 0) = 



a(a +1) . . . (a +r -!)/?(/? + 1) . . .(ft+r- l) 
r!y(y + l) . . . (y+r-1) 




provided that y is not zero or a negative integer. 
Next take the alternative index a = l—y, which gives 
1 -y) =c 0 * 1_7 


i+2 


(a-y-f 1) . . . (a-y-fr)Q—y-f 1) , . . (fl-y+Q 
(2-y)...(r + l-y)r! 



provided y is not a positive integer greater than unity. 

When y = 1 the two particular solutions are identical. 
Wo shall suppose that y is not an integer, positive or 
negative, or zero and use the notation 


j 8 ; y; x) = 1 + 


a • , a ( a + l)^(^ + l) 

1 . y 2 ! y(y + l) 


• • 


for the hypergeometric series. Writing A and B for c 0 and 
c x respectively, we may express the general solution of the 
hypergeometric equation as 
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y = AF(a, P; y\ x) + Bx l -’F(a -y + 1, P~Y + 1 2 V > *)• 

The above procedure is known as the Method of 
Frobenius. 

Note 1. It is possiblo by a modification of the above 
procedure to obtain the complete solution when y has one of 
the excluded values, but the results are not of any particular 
interest. The essential feature is that the solution that 
becomes irrelevant is replaced by one involving logarithmic 
terms. This point wiU be illustrated later in the case of the 

Legendre equation. 

Note 2. The general solution valid in the neighbourhood 
of the singularity x = 1 is 


y=AF(a. pi a+^-y+1 5 1 “*> a 

+ I3( \ -x)'t- a ~ p F(,y-a, y~pi y-a-0+1 


1-x) 


provided that a +P -y is not an integer or zero. 


55. The Legendre Equation and the Function P n {x). 
The equation 

(1 -x 2 )y' +n(n + l)i/=0 . . (oo.l) 

has the finite regular singularities x — ±1 ; in order to study 
the solutions valid near x = l it is advantageous to make 
the change of variable x = l-2f, which transforms the 

equation into 

« (1 - {) S +(1_2f, 5f +n( ’ i+1)y=0 ' 

This is the particular case of the hyporgoomotric 

equation in which a=n +1, Y } ’ 111 ,lt 1 Jf 

case in which the two fundamental solutions corresponding 

to (64.4) become identical, and wo are left with 

« = F(n + l, ~~ n '» 1'» 0* 
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The corresponding solution of the Legendre equation, or 
Legendre Function of the First Kind , is defined as 


P n (x)=F(n +1, -n\ 1; £-$x). 

We shall confine our attention to the important case in 
which n is a positive integer; since the term in £ n+1 and 
every succeeding term in the hypergeometric series 
F{a, jS; y; £) contains the factor j3+w, it is clear that 
F(n + 1, -n\ 1; £) is a polynomial in £ of degree n. 
Hence 


P n (*) = 1 - 


n(n +1) 1 -x {n - l)n(n + l)(n +2)/l -x\ a 


+ 


+ . . . 


1 2 ' (2 !) 2 \ 2 

1.2. . . (n - l)?i(n +1) . . . 2»/l - x' 

■ • • + <->"- (,!)■ -(“I 

is a polynomial in x of degree n such that P n (l) = l. In 
particular, 


P 0 (:r) = l, P l (x)=x, P,{x) = \{ 3x=-l), P 3 (z) = -3*). 


The final term in the above expression for P n (x) is the 
only term that gives rise to x"; we see therefore that the 
coefficient of x n is (2n) !/2 n (n !) 2 . The whole expression 
may be expanded and then rearranged in descending 
powers of x, but it is less laborious to substitute the 
polynomial 


y=c 0 + c 1 x+CnX- + . . . + c n x" = Stye*- 

in the left-hand member of (55.1) and to equate the 
coefficient of each power of x to zero. Taking the coefficient 
«»fx r-2 wo have 


r(r - 1 )c r = (r + n - 1 )(r - n - 2 )c r _ 2 . 

As c n is known, we work backwards and find 

n(n-l) n(n - l)(n-2)(n-3)^ 

C "~ 2 “ " 2{2n - 1) C "’ c "~* = 2.4(2n-l)(2n-3) C "’ 

f _/ r »(»-1) . . . (re-2r + l) 

•»- *«■ t ) 2r r! ( 2w - 1 )(2n - 3) . . . (2 n - 2r +1 ) Cn * 
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Tt is easily verified that the coefficient, of 1 in the ex- 
IZZ H P„(*) is zero, Iron, which it follows ,n success,on 

that c„_„ c„_ 5 , ... are zero. Hence 

( 2 n) ! r ”(n - V _ , , Mn - l)(n - 2 )(n J)^_ 4 _ _ 

2(2^1)^ + 2.4 . (2n - l)(2n -3) 

^ * a 


JV 

= T(-)' 


(2n - 2r) ! 


r«0 


_ r n-2r 

r * 


2»r!(»-r)l(n-2r)! 

w her‘eN = *n for n even, = for » odd - , 

Thus the Legendre polgnomial PJ& «> even ™r odd 
according as its degree n is even or odd, smee P.(D 1. 
wo conclude that P n { 1) ( ' 

sou. An important formula for P„(x) may be deduced 
from this last expansion. We observo that 


and therefore 


( 2n-2r) ? n _ 2r _ £L x t»-*r 
(n - 2r) • dx- 

d" v (-) r J > "- ,r 
P^ x ) cLr n ^ 2"r!(n —r)l 


where r ranges originally from 0 ,o S We^may 

range of summation, however, r n( j their n w doriva- 

^v^sh 8 ™w O 0 f carTomploy the Binomial Theorem and 

obtain the Rodrigues Formula 

• In 


P n W 


1 _ 

2" . « ! d.c n 


(x*-l) n . 


We may prove from U.U formula that IV«1 ha., = y,st,nct 

2 ™ nT Bo'lhat ’ /„ - *■’« I P.W- ThCn • "'fJTn 

by induclion’that/ has at .easy,, "“Tt'li and i 

r/^^rr-^rnSinct 1 zeros them. 
The induction now proceeds. 
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56. Solution for Large Values of | x | . Physical 
problems not infrequently demand the solution of a linear 
differential equation for large positive or large negative 
values of the argument x. In such cases we use the 
reciprocal transformation s = l/£ whereby (53.1) becomes 


Pie- 1 ) 


( &+2e^ 

* dp f dt 


+Q(e~ l ) 


_A ftl 

dt 


+-R(£ _1 )y=o» 


retaining its linear form. The point x=co is said to be 
ordinary or singular according as £=0 is an ordinary or 
singular point of this transformed equation. 

The Legendre equation, so transformed, becomes 


(1 




£?+2t# 

dp^ ? df 


r-2^- 1 * 


-p 


& 

df 


+n(n+l)y=0 


or 


% = f ! ( « 2 -1 +2 ?%+»(" +1 )V =»• 


Assuming 
we find that 


(c 0 #0) 


y=Zc r P +r 

L ( y=c 0 {o+n){o-n-l)p 

provided that c 2 , . . . are so chosen that 
(or +r —n +1)((7 +r +n +2)c r+2 =(o +r)(a +r + l)c r (r>0). 

The equation is thus satisfied if a = —n or n + 1 ; when 
a = — n the solution is of the form 

y=c 0 r n +c 2 r n+ 2 +c 4 r n+4 ... 

with the recurrence-relation 

(r+2)(2a -r-l)c r+2 = -(n -r)(» -r -l)c r . 

Reverting to the variable x we have the first solution 


Vi =c o 


| n(n-l)(n-2)(n-3) 
2(2n —1) ~ 2.4(2» —l)(2n—3) 


* t • 


and taking c 0 =(2a)! /2 n (n !) 2 when n is a positive integer 
we identify this solution with the polynomial P n (x). 
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When o=n + l the solution is 

y =c 0 r +1 +c t r+ 9 +'*r"+ ••• 

with the recurrence-relation 

(r+2)(2n+r+3)c f+2 =(n+r + l)(n+r+2)c f . 

Thus we obtain the second solution 

(n + iKn-f 2 ) 3 

+ 2(2/t +3) 

( n +l)(n +2)(w +3)(n +4) r _ n _ 6 + 

“ - - - — . st\ . e\ 


V 2 “ C o1 X 


..1 

2.4(2n+3)(2n+5) J 

When n is a positive integer this series does not terminate, 
but since 

c r+2 _ (w+r + l)(w+r + «j as r _> a,, 

c r (r+2)(2n+r+3) 

r i „ | - 1 Takine c«=2 n (n !) 2 /(2n + l) ! 
it converges for |x|> »• J aKU, S <0 v ' 

we obtain the second solution 

« 2 "(n+r)! (n+2r)l r _ n _ 2f _ 1 

Qn( x )~ Zj r I (2n+2r + l) I 

which is known as tho Legendre Function of the Second 
Kind. It is of opposite panty to P n (x). 

N oU. Tho method explained in 5 38 enables us W d t 

.nine a second solution by quadrature. Sub.tltut.ng (65. ) 

vdx 


= Pn(*)\ 

we find for v the differential equation 


2 P n 'IP» ~ 


+ 1= 0. 


from which 

Thus wo have a so 
form 


v-C( 1 

con d solution of Legendre's equation m the 

f dx 

/%(X) J (I 
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It is possible to prove that this is equal to Q„(x), provided that 
the lower limit of integration is taken at infinity. 


57. The Bessel Equation and the Function J n (x). 
Unlike the Legendre equation, the Bessel equation 

x*y”+xy'+{x 2 —n 2 )y=0 . . (58.1) 

is not essentially a particular case of the hypergeometric 
equation. It has only the one finite singularity x =0, 
which is regular ; the point at infinity is also a singularity, 
but is irregular. 

To obtain a solution in series by the Frobenius method, 
write Ly for the left-hand membor of the equation and 

<f>{x, a) = ^C r Tf+ r (Co^O). 

Then 

Lf>{x, a) =c 0 (<7 3 -n 2 (af +c 1 {(a +1) 8 -n 8 }af +1 , 
provided the coefficients c 2 , Cj, . . ., c f , . . . are such that 

{(a +r) a —n 2 }c r +c r _ 2 =0 (r > 2). 

If y=<f>[x, <r) is a solution, Lf>{x, a) must vanish identically, 
and therefore 

c 0 (o 3 — n 2 )=0, Cl {(a + l) 2 -n a }=0. 

Since c 0 ^0 we must take o = ±n and consequently e^O. 
The relation between c r and c r _ 2 now shows, taking 
r =3, 5, ... in succession, that all coefficients of odd rank 
vanish. 

Taking first of all a=n and writing r=2s we have 

4s(n +s)Co, = -Co,_ a {8 ^ 1) 


and thus wo obtain the solution 


y=<f>[x, n) =cJ 


•n+2 


X n — 


2 2 (n + 1) 

j- 


+ . . . 

(-Yx n * u 


2 2$ (n + \) . . . (n-f$). s I 


+ 
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The coefficients are finite except when ,, is 
integer. Excluding this case, we sUndard.se th > 
by taking c 0 = l/2"n» + l) in general and C.-1/2 
n is a positive integer. vSince 

. 0 as s - 


2 8 


7) 


r 2j _ 2 45(n+5) 

the series converges for any finite value of r, and thus we 
have the first solution //=•/„(*), " herc 

» (-y 

! r(n +« + 1) - 

in general. In particular, when n is a positive integer, 

(-)• ( X \* U 
./„(x)-2-.rr(^+7ri\2/ • 

This function J„(*) is known as the Bessel Fenchnn oj 

* not°an C intager the second solution ntay he 

obtained by replacing n by n, it is t u re 


HL- , 

^-.(*)-2-rrr(*-« + n\2y 

f / nnd •/ (-0 ftr< ‘ respectively 

The leading terms of J.M »" thc ono is not 

finite (non-zero) m “ ltl ^ C3 ,[ r an( , therefore the general 
a mere multiple of the oth cxpre ssed as 

solution of the Bessel equation maj 

But when n » an intC f is no k^s’of gencrahty 

differential equation only a- (x) is not distinct 

in taking it to be a positive in l. ■ . - fll . st n terms 

from J n (x). For the denommat ^ , he factors 

of the scries for ■/_,.(*) ‘ ‘ * are infinite, and so 

HI -n), T(2-»), • • *. 1 
these terms vanish, ihus 
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* (-)* lx' 

= 2,! P(5 - n +T)\2 


25-n 


(-) r+n 


2r+n 




2( r + n )!r(r + l)V2 

In this case a modification of the method of Frobenius 
is necessary to give an independent second solution. 

58. The Function Y n (x). Consider the functions 

(-)* / x \fc + n + e 


«W*) = 2 

*7-n- e (*) = 2 


s ! r(s + n+€ + 1)\2 

(-)' 


h -n- t 


8 ! r(s -n -c + 1)\2 

where n is a positive integer and e a real number numerically 
less than unity. Both satisfy the Bessel equation 

xhf + xy' + [x- - (n + e) 2 }y = 0 
and so does the function 

»-«(*) 


So if Ly = x*y” + xy' + ( x z - n 2 )y, then 

Ltj ={(n + € ) 2 - n 2 }77 =c( 2 n + € )r;. 

If, therefore, 77 has a limit as e -*■ 0, y = lim 77 is a solution 
of Ly = 0, i.e. is a Bessel function of order n. We prove 
that the limit exists by making the following transformation: 

lim ^ n * e ^ ~ ^ ~ W-n-c( x ) 


«-*-o 


€ 
lim 


and evaluating these partial derivatives. Some authorities 
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“i" d.i, «h. i-- - -i 

thU9+ , ipjW ( _ r a7 -> (x) l • 

Since 

^ (_)5 _/*' 2 ' + ' 

J ^ x) = 2 JTf(7Tv + 1 )\2 


LiZ — ( -Y <+ V°g & x) -'i'(* +v+l) }' 

s ! r{8 + v + 1) \2 


-2 


we have 

dJM) 
dv 

where */-( 2 ) = r'( 2 )/r(z). Hence t 
dJAx) 1 ^ ( - )* 


2 *+n 


r&W| ,; K 

L aw _U„ Z 5 !( 5+ n)!\2 i 

x|iog(ix) + y _1 " • • * "^+»J 

- (-)* 

„ J n (z){l°g (J*)+y} - Z-TToTmO”! 


[ \fs' 25+n 

(i+l + ■ • • + I77.\2 / ) 


m 

is 


i ^ Ttriiiah Association Maths- 

• This definition was adopted intho quiva lont definition 

a Heal Tables, vol. vi (Cambridge, 1937). A 4 


Y n (x) = lim 


J^(x) COS VTT -J -y( X \ 


gin vn 

t Since l’(a + 1) -*!’(*>. w0 h * iv ° »’ (f + + 1/z 

and bo when z is a positive integer, i 

v . (2 + l)= V .(l) + l+i+l + • ' • + z 

. na 779 ir.titi . . (Uulor’fl constant). 

where -y(l)-y *= 0-677216b0 • • 
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It will be remembered that the first n terms in J _ v {x) 
vanish when v =n\ we therefore write 

, , , (-)* ! X \" (")* /*Y*" 

" ,f-„ s! r(s-i'+i)\2/ ^3ir(s-v + i)\2j 

In the first sum we use the relation 


r(z)r(l -z) = 7 r/ sin Z 7 r with z-s-v + 1 ; 


in the second we replace 8 by 8 +n, and thus obtain 


n - J ( — )*sin (s-i' + IJtt . T(v-s)/x N ^" r 
- - W = 2a - 

l-O 


TT . 8 ! 


+2 


(_) S+ n 


U + in - • 


,To ( 5 + n )! T( 5 + n - v + 1)\2 
Differentiating with respect to v , we have 

dJJx) n - 1 ( - )• /x\-* -»' . 

—-= T- : \ sin {s - v + 1 )rr . T(v -a) 


?v f 'toW.a!\2 


{ 7TC0S(.S-V + 1)7T , , tx' 

- sin (,,-„ + !)* U 


00 


+2 


( - )‘ + " 


a; 


2t + 2n - r 


+ I> + »-i' + 1)\2 


{ 


x\ifj{s+n-v + 1 )-log ( o 


8J-y(x) 
<d\ 


(*)1 n - 1 ( » 

r~ ^2-Trio r(»-«) cos («-»)* 

Jv-n ,_ 0 « ■ \-/ 


00 


+ 2 


(_)»+n / X \2.<+nf 


2 p+D->ogU 
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(- r ' 


i -0 


al . W ^ 0 «'(*+n)!\2 


2«4 n 


x |log + y - 1 \ 



• • 


-!} 


2i- n 


(-)"Z-(-) n - 7 " (x H(**>+y) 

+ (-D"i s -iiV«)i( 1+l+ --- + *)® ' 


#-0 


Thus the function 7„(x> or AmI Function oj the Seen,I 
Kind of order n is given by 

dJ_ v {x) 

L * ~ y\"' i / \ n » — 

-a 




-2./„(i){ log (1*) +y) - 2t 71 w 


,-0 


« ( _)* 

V 1 ' 


1 


1 


1 + | + • • • + 


+ 1 + i + • • • + 


n +«/ - 



n 


4 2 # 


'• (n + a)! 

and when • is a positive integer the complete solution of 

the Bessel equation is 

y -AJ n (x)+BY n {x). 




EXAMPLES 

T (!■' - ill/' =ty* + v 2. xy’ +y* = l 

7 9 : 8 - 
II. (**+o»)J/'-<»+»«* +y<*'+«*» 

I4 . t,x' - 2 xy -y')y' =*' + 1 *y - V' ' 5 - <* f *5,' 
16. x(2- - by')dy = 4 y(*> + 3^2 17 - » <f + /_ > 

i8 xhidx - (x* + ay 3 W{/ = 0 * 9 - . 

20. (x»+x»/ +ay')y , =ax t + xy +y* 21. x(x+y)dy=(x +y 

22. *(*‘ +a*y +2/*)!/' =y( x ' + bxy+y') 

23. (x+y)V=x*- 2 xy+ 6 r/‘ 24. xy -y-x cos (Jl ) 

25. ( 3 y-x)y'= 3 x-y +4 

26. (*-6w + 5)d* + (6*-!/ + l)«l/- 0 , \_ 

27! ( 4 y +2)7- =y - 4 X 28. (92 + 2 y + 19 )<ly =(22 + 6y - IBM* 

29. ( 9 x + 21 i/ + 3 )i/'= 7 x- 5 i /+45 

30. ( 8 *+y+ 25 )dx+( 7 x- 16 y + 140 )rfy -0 6 

£ s ^ w ’- s * 

35. (*-y)y=(*-y + 1 )* 3 6 . (x+y)*dy-(*+y 

Show that 37-44 aro exact, and integrate them: 

37. (ax + hy)dx + (hx + by)dy = 0 

38. ( 3 x* + \xy)dx + ( 2 x* + 3y*)dy =0 

39. {3ax* + 2(a + 2 h)xy + (b + 2/<)y V* + {(a + 2 h)x + 2(^21)^^ 

40. ( 4 x*y - 12 x'y 1 + 5 x* + 3 x)«/' + Ox V - 8 xy 1 + 10 xy + 3 </ = 0 

41. xV/y+2x log | !/| =° 

42. (cos x - x cos - (sin 1/ + y sin x)dx - U 


43- 


x+a, (x +a\* . 


(1 -y a )<fajj l x l )dy _ ^ 
^ (1 +xi/)* 


131 



132 


DIFFERENTIAL EQUATIONS 


Show that each of equations 45-56 has an integrating factor 
of one or other of the types y(x), n(y), y(x+y), y(*y)> 811,1 
henco integrate: 

4S . (x +y)dx +dy =0 46. (* s +y*)dx+Sxy^y =0 

47. xy'dx + (x'-y'-l)dy =0 48- (1 -xy)dx + (l -z')dy =0 

49. y 3 dx+(xy + l)dy =0 5®- dx+ {l+(x+y)tany}dy -0 

51. (x*- 2 y*- 3 xy)«ix+ 3 x(y*+x)rfy =0 

52. (x*-y* + l)dx + (x* -y t -\)dy =0 

53. (2x* +xy + a 3 )ydx +(x + 2y)(x* +a*)dy =0 

54. ( 3 xV - y')dx + ( 5 x*y 7 +x 3 )dy =0 

55. [y+ \)dx+(xy+y*+y+ l)dy =0 

56. (lx 3 + 3 x*y + 4 y)dx + ( 4 x* + x + 5 y)dy = 0 

57 .xy'-y=x 3 + \ 58. xy' +y=x log | x | 

59 - (x + l)y' -y = 3 x 4 + 4 x* 60. (x*+a)y' -xy=a 

61. (1 — x*)y' — xy — 1 62. (1 -x*)y'-y = l-x* 

63. (1 + x*)y'+xy = 3 x» 64. ( 2 xy' +y)\/(l +*) = 1 +2x 

65. V(1 +x 3 )y'+y = 2 x 66. 2x(x - l)y'+(2x - l)y-x 

67. x*y'+xy=x* +x + l 

68. 2(x s +x + \)y' +(2x + l)y = 8x* +1 

69. (x* - l)y' - xy =x* 70. x(l-2*)y'+(2x*-l)y=x*-x» 

71. (1 + x*)y' + 2xy=tanx 

72. 2(l-x%'-(l+*)i/ = V(l -**) 

73. 2(1 -x)y'-y= 4 x\/(l -x) 

74. (l+x*)*y'+(l+x)(l+x*)y = 2x 75. xy'+2y=smx 

76. y' cos x + y sin x - cos* x = 0 77. y' +y tan x =sin 2x 

78. y' sin x cos* x + y cos 2x cos* x = 1 

79. y' + y cos x = i sin 2x 80. y' sin x - y cos x =e* sin* x 

81. y' -y =x +sin x 82. Jy'=y tan 2x +1+sec 2x 

83. y' +2y =x l + 3 cosh x 84. y' cos x +y sin x = 1 +tan x 

85. y' + y tanh x = Ge* 1 86. y' sin ! x - y tan x = tan x - tan* x 

87. y' cos x - 3 y sin x =cot x 

88. y' + 2y cosec 2x = 2 cot* x cos 2x 

89. y' tan x +y =sin 3 x + cosec 2 x 

90. y' cot x - y = cosec 2 x + co3 2 x 

91. (1 -x*)y' +xy =x arc sin x + (1 -x)\/(l -x*) 

92. x(x + l)y' - (x* + x - l)y =(x + l)(x* - 1) 

93. xy' = 4 y - 4 \/y 94. 2 xy'-y = ( 2 x* - l)/y 

95- xy'=y + 2xy* 96. (1+x*)y'+xy =xy* 

97. x*y' =2x*y 4-y* 98. xy' +y =y* log | x | 
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ioo. x*y' + xy + Vy - 0 

y' cos x +y sin x + y*=0 


102 . 


99 - xy'+y =x y' 

i°„3. TUylw- 104. 2(1 W + **- 3-+»-=0 

Equation 105 has a particular solution y- 1 . 106 has 
y ={x + 1 )/x», and 107-110 have y =2; complete the integration. 

105. y'={y-i)(*y-y- x ) xo6. y'+y 5 -*" 

107. xy' +(m-l)ax m (y-x)*+y-2x-0 

108. y' +xy* - (2x* + l)y +x s +2 - 1 =0 

109. (x*+a)y'+ 2 y*- 3 xy-a =0 
no. 2x*y' = 2xy + (y* - x*)(x cot x - 1 ) 

Integrate 111-118 by means of a change of variable: 

in. x(y'+ 1)+tan (x+y) =0 112. xy -y=x*+y 

113. (1 - x*)yy' + xy* + 2x* = 0 

114. 2(1 + x)yy' +y* + 2 x - 3 x* =0 
1x5. (x*+xy + l)y'=(y*+xy + 1) 

116. (x>+y»-a)yy' + 2(2*+y*+a) =0 

117. {(x* +y* -a)x - y)y' =(x* +y* -a)y +2 

1x8. (x*+y*-a)(x+yy , )=2xy(y-xy') 

119. (x* — y*)y' =xy (Substitute x =ty ] 

120. (xy' - 2 y)* =x 2 (x i - y ') [Substitute y =*x*] 

121. (x*-l)(xy'-y)= 2 x(y*-x») [Substitute y =<x] 

122. Substitutey=x(l + z)/( 1 -z)inx m (xy' -y) =y -* • Hence 

integrate when (i) m =0, (ii) rn = 1, (iii) — 

Find the equations of families of curves having the following 
properties ( 123 - 132 ): 123. Radius vector equals tangent 

124. Normal bisected by y-axis 125. tangent has con 

stant length a (Tractrix) 126. Area of trapezium bounded 
by axes, tangent and ordinate = a* X27. P rojec tion o 

ordinate on normal =0 128. Subtangent NT equals 

abscissa of projection of origin on tangent 129. Mid-point 
of normal lies on y* =ax 130. Ordi n ate N P meets y =ox 4 - 
in Qi R is taken on NP so that PR =NQ; OR is parallel to 
tangent at P 131. NP +PO = 2 a 132. Arcs—2*/a 
If the lino tlirough 0 perpendicular to OP meets the tangent 
in U and the normal in H, OU =r tan <f>=r i dOjdr is the polar 
subtangent, HO =r cot <£=dr/d 0 the polar subnormal, PU the 
polar tangent, and PH the polar normal. Find curves such 

that ( 133 - 138 ): 

X 33 . HO=k'OU 134. OUlHO=y'lz' j 35 . PH =a 

K 
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136. H. describes the line r cos 0 =p 137. « = V( r * + 2or) 

138. 4 > =m 0 

Find the orthogonal trajectories of the following families 
of plane curves in which the parameter is a ( 139 - 157 ): 

139. x 3 -3x2/*= a 140. y* = ax 141. x* + ay* =0 

142. x* — y 3 + 2axy =a* 143. x(x* + y 3 ) +a(x*-y*) =0 

144. (x* - a)* + (y* - b ) 9 =a* 145. xy=a(x-l)* 

146. y = a log | x | 147. y =a(x + a)e* 

148. tanh* x + tanh* (y + a) = 1 149. x*/a + y* =a s /(a - 1) 

150. Rectangular hyperbolas through (a, 0 ) with y-axis as 

asymptote 

151. Ellipses, centre 0, common major axis 2a along Ox 

152. Conics circumscribing rectangle of angular points 

(± a, ±b) 

153. Parabolas with common focus 0 and latus rectum 2 i 

154. r=tan(0+a) 155. (r*+a*) sin 0+ar = 0 

156. r =a/(l +2 cos 0) 157. r=a(l+2cos0) 

Find the trajectories at angle (0 (tan co=m) of 168 - 160 : 

158. x*-y* = a 159. x* + y* =2ax 160. r p =acosp0 

Integrate 161 and 162 , and consider the run of their integral 
curves near the origin: 

161. y' =(y -x*)/mx 162. y'=(y*-x)/y 

Solve 163-166 for y' and integrate: 

163. Tty' 1 - 2yy’ -x =0 164. cos x (dx 2 -dy 1 ) = 2 sin x dxdy 

165. (y'-x)*=y'+x 166. x*y' a -xyy' +y* =0 

Obtain the singular solutions of y=px+/(p) when /(p) has 
the forms ( 167 - 176 ): 

167. p 3 168. p/(p + 1) 169. 2y/(ap) 170. V(P l “ l ) 

171. 5(p + l) 3/l 172. \/{(a+2ltp+bp 1 )l(ab-h 1 )} 

173 - «W \/(l +P*) 174 - p - p log | p | 

175. \/(l -p*) -p arc cos p 

176. y(p* + 1 ) - p log {p + \Z(P* + 1 )) 

177• y = 3 px+p~ , y~ 4 [Substitute y 3 =v] 178. y=p+e'/p 

179. Transform (1 + p 3 )y* =/(x +py) by the substitution 

2s =x* + y\ and hence integrate 

180. Singular solution of (1 + p 1 )y 1 =(x +py) 3 

181. x 4 p 3 -x 3 yp 3 -x*y*p + xy 3 = 1 [Substitute y =vx, x*=f -1 ] 

182. y=2xp+p* 183. y =x+ap-a log | p | 

184. x - y = 2ap -ap 1 185. 3 p s y = ( 2 p 3 - l)x 
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x8 6. y =P’ +P‘ « 87 ' y 'f g ' 0 . + (p!x)' + P<“ : *> = > P 

• then 

192. Show that, j £ (0 * +^ an ()d( ’ S ec ()} CQ3 1 

where tan «=P. ^"Tfl r,r sin 9 + a) =0 

193. Integral curves of r sin « Q f ( ( 19 -t- 197 ): 

Writing express 2 + ', 

194. (P+P> - y ~* 0 ig7 . 9*+p‘=!/P(P+p) 

‘""^"^'n 1 such that 

W ^ + n = 0 siiaU have ^ a. s.nguiar 

201. P* + P ■ 7 1 20 - ipp- + (p* + o)P = 0 

203 • **£“?, _ 206. PP' + P" + 20 V = 0 

20 ?'. P'(P_-1)=2P'/ 20 x io M 2*pp-=2(2p-)‘ + uxpp' +ap' 

zog. xyy" + ny ( x y z(l -*)«/' + 2(1 - “ 1 

2 X 1 . (*W* “»'•)*•+J “® „ n < J the general integral of 

2x3* Show that the for f ^ on the nature of the 

3 —■ - 

throe cases , . F(u '\ =0 has a first integral of the 

2I4 ' ***** whon 

F(P)=(1 +P)(1+P'» +1 IDay b0 reduced 

2,S ' ^r^r formTyCtbsiftution 

226. p/r^T" 216 - 

218. PP'=P-+PP' famil “ 9 'of curves such that ( 219 - 224 ) = 
Find two-parameter famil f cM , #=a 223. (>= r0 

224. Ordinate through T bisects radius of conlro 0 ) 

,25. Find the inverses (with respect to the ■" 2fI 8ec , „ 

of the integral curves of rr -2r +r 
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226. Form the linear equation whose general integral is 

y=A[\ +x) m +B(l -x) m 

Show that 227-229 have a polynomial particular solution 
and deduce the general integral: 

227. xy'-(x + 3 )y'+y =0 228. (1 + x*)y'- 2 y =0 

229. x(x - l)y' + (2x - ljy 7 - 6y =0 

230. x( 1 + x)y ' = {(n - 2 )x + n}y' + ny [part. sol. y = l/( 1 + 2)] 

231. x(l -x) t y M = 2 y [part. sol. y=x/(l -x)] 

232. y' + 3 y' + 2 y =e x + sin x 233. y' - 3 y' - 4 y = 10 cos 2 x 

234. y'+ 2 y'+ 5 y = 8sinhx 235. y" + 6y' + 9 y =e~ >x cosh x 

236. y" - by' + 6y = 4 x*e x 237. y'+ 2 y' + 3 y=e-'coax 

238. y" - 4 y' + 4 y =e ,x cos* x 239. y'" - 3 y'+ 2 y=x*e a 

240. y* -2y' sin a +y =2 cos a sin x [ct^mrfe) 

241. y' cos* a -y' sin 2a +y =x*e z 1411 a 

242. y ,v + 2 y n +y = 24 xsinx 243. y vl -y lv -y 11 +y = 4 e* 
Particular integral of y“ +y=/(x) when /(x) has the values 

( 244 - 249 ): 

244. e x sin2x 245. e ,x cosx 246. e x (x*-1) 

247. 8 cos x cos 2 x 248. 4 x sin x 249. x cos x - x* sin x 

Ditto for y" - 2 y' +y =/(x) ( 250 - 265 ): 

250. x s -6x* 251. e x (l+ 2 x+ 3 x*) 252. e^x 1 

253. e x sin x 254. 50 cosh x cos x 255. 8x*e 3 * 

256. y*- 2 y'+( 1 +m*)y = ( 1 + 4 m*) cos nix (y = l, y '=0 when 

x =0) 

257. y' + 3 y' + 2 y = 2 xe~* (y =y '=0 when x = 1 ) 

I58. y'- 2ay'+ (o*+m*)y = 2;ne ox cos wu: (y=0, y'=m when 
x = 0) 

259* y'"-y'-y'+y = 8x6-* (y =y' =0, y' = 1 when x =0) 

260. y lv - 2 y* + y = 12 xe r (y = y' =i, y' = 0 , y'" = -3 whenx = 0 ) 

261. 2x l y'-xy'+y =x* 262. x*y'- 2xy'+2y = 2x log | x j 

263. x*y*+ 4 xy'+ 2 y=log | 1 +x | 264. x*y' +xy' +y = log |x| 

265. xhy' - 2xy' + 2y =x s sin x 

266. 2(x + l)*y' - (x + l)y' +y =x 

267. x a+, y* -( 2 a - 1 )x°y' +o*x' ,_1 y = 1 268. 4 x*y' +y = \/x 

269. x*y'- 2nxy'+n(n + l)y =e x x n+l 

Transform 270-273 by the substitutions indicated, and 
integrate: 

270. y' cos x+y' sin x =y cos* x [f =sin x] 

271. (I -x*)y'-xy' +y=0 [x=sint] 
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272. (1+*’)«/'+*«/'= 4 y [x=sinh/] 

271 (1 + x*)V + 2x(l + x*)y' + 4y-° [*- tan vl 
274 Find a substitution t = *(*> such that the coefficients in 
2p(x)y' + p'(*)y' + y = 0 become constants 

275. 4 xy'+ 2 y'+y =0 

\i* sin 2x + 2y 7 cos 2x + 4 y cosec 2x = u 

277. Find a transformation y=v<f>(x) by means o^ " 

77 y' + p(x)y' +g(x)y =0 takes the form a +r(x)»-0 

278. xhf 0 -2nxy'+{n(n + l)+in t z)y -0 

279. y” sin* x - 2y' tan x + (2 + sin* x)y -0 

t, in ,' + W=<> 

282. Relation between p and 9 such that y" + 1>y 

shall have two real solutions (a) whose product is 

unity, (b) the sum of whose squares is unity 

*>■>•> by “ 

integrals 

Particular integrals of 285-287 by variation of parameters: 
285. t/'+4y =2 tan x !/ +y- log | co.-> x 

287. 2x*y' + Ixy' + 3 y = c ° 9 V x 

General integral of systems 288 - 293 . whore x =<*>/*’ 

288. x+J/=sin 2f, y-x=cos 2f 

289. x= 4 x- 2 y +e‘, y= 6 x- 3 y +e _ 

2 2 9 9 x! * + ‘ \*— l J = 0 when < = 0) 

29 , = -0 — 

293. i + 3 y°- 4 x+ 0 y = l 0 cosf, x+y- 2 x + 4 y =0 (* 

x = 4, y = 2 when 1=0) 

Sol u of 

296. xy' +y' - 4 xy =0 (solutio " 9 . ^ “ “'“1 % =0°!uch that 

^ Fi ;i 1 w°r "e^ng + ^y X; , i *>• that 
a second solution is x'-F (a - y + 1; 2 - y i x) 
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298. Transform 297 by the substitution y=e t v and find a 

series for e~ x F(a; y, x) 

299. Second solution of 297 when y = 1 

300. Show that x*y'" + 32 ry’ + (1 - 4 n* + 4 x*)y' + 4 xy = 0 is satis¬ 

fied by the square of any solution of the Bessel Equation 
of order n, express J n (x)J_ n (x) as a power series 


SOLUTIONS 


x. xy=C(x-l)(y + l) 2. 

3. x 3 + 3 y + 6 log | 1 -y | =C 4. 
5. y J + 1 = Cx/(x + 1) 6. 

7. y = tan (C + arc sin x) = {c \/( 1 - 

8. y =(C cos* x - 1 )/(C cos* x + 1 ) 
10. x* + y* =log {Cx*(l +y 1 )} 

xi. y + b =C{x* +a* +x\/(x* +a*)} 
12. y +a log (y* +a‘) =C +x - 1 /x 
14. x* + y* = C(x + y) 

16. x* + 0 x t y i =Cy 
18. y = C exp (x*/ 3 ay s ) 

20. (x*-f xy+y J )‘—= C(x-y)*+ a 


1 +y = ( 7 x*(l -y) 
e** =Cx cos y 
v/(x*-l) + V(!/'-l )=0 
x 1 ) +x}/{V(l “*•) -cx} 

9. y =C(x* -x)/| x* - 1 


13. x* = 2y* log | Cy | 
IS- y(y+ 3 x) 5 =Cx* 
17. x*=C*- 2 Cy 
19. ( y-x) l ~ a =Cxy 
ax. e» /r (x - y)* = Cx 


22. y’ = Cx" exp {x/y - y/x} 

23- (*“!/)' log | C(x -y) | = 4 xy - 2 x* 

24. exp tan (y/x)=C/x 25. (x-y + l)(x+y+ 2 )* =C 

26. (x + y - 1 )* = C(x - y + 1)' 

27. 2 log (x s + y J ) +arc tan (y/x) = (7 

28. (x- 2 y + ll)*=C( 2 x+y + 2 ) 

29. (x- 3 y + ll) 4 (x+y+ 3 )*=C 

30. (x + 2y-10) 6 (x-y + ll)»=C 31. (x - 2y)» + 6x - lOy =0 

32. (« + 1 )(y - x) + 26 log |(a + l)(y +ax) +6(0 - 1 )| =C 

33. y - 4 x - 1 = Cx*(y - x - 1 ) 


34. log | y+a | + 2 arc tan {(y + a)/(x +6)} —C 
35 - (x-y)*+ 3 x+y + i log | 2 x - 2 y + 1 \ =C 
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37 - 

39 - 

40. 

42. 

44. 

46. 

48. 


«.+2ta,+6|r-0 38. x ‘ + y' + l*'y=C 

(*+»)(«■+ 2**»+<>»*) = C *• log I » l=C 

2 x s v* - 4x *y* + 5 x*y +3 xy=C 41. x 1 g 1 y 1 

9 • n 


50 
52 
54 - 

56. 

58 . 
60. 
62. 
63 - 
65 - 
66 . 
67. 

69. 

70. 

72. 

74 - 

75 - 
77 - 

79 - 

81. 

82. 

83- 

84. 

85. 

87. 

88 . 

89. 

90. 

91. 

92. 
94. 
96. 
98. 


y cos 2 - x sin y — C 
1 +xy =C(x +y) 

(x* +’iy*)V x= C 
arc sin x + yV(l 
(x +1/) sec y=C 
(x — 3/)e* < - r+,> * =0 
2xV +x-* -y~ i =C 
(x* +y)(x+y)*=C 
y=C/x+bx log | x | 

y=CV( xt+a ) +X 

“ I . . V 1 


43 . y+ 6 =C(x+a)* 

45. (x+y-\)e x =C 
47. x 1 !/* = log Cy* 

49. ye'" =C ^ 

51. x* + y* + 3xy — Cx 

53 . y(x+y)v/(^+ a *) = C 
55 . e v (x +y) =C(y + 1 ) 

57 . y =Cx + ix s - 1 
59 . y=C(x + l)+x* 

VI.-•;-«*«*■■■ 

; -i%‘ rift'v ’ ’Kiwi 1 .. 

» Q*f| — X“ l cos x J ' - 

y=x *(C +8inx) X _c cosec 2x + sec* x 

y = G C08X-2C08-* 7»- V = 8in x 

y=Cc- B,D *+sinx-l y v 

’ y =Ce* - X - 1 - i( s * n X +C03 x) 

y =(C + log sin x) sec* x + i sec x 

» X X x -C - + K*4 «£ + 1 log | .an lx | 
l Zllc 4 Z I tan (J»+1*1 I + i 31 

" = V<!-;)«;+ J* V( 1 - *■» + « “ * + ,,, 

y*=X*+Cx + l 95 y Jtr/n 

y={C\/(l +X*) + M -1 97- I/=^ {6 J X U 

y ={CX + 1 + log 1 X |) -1 99 - 1/=^ +24 
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xoo. y =x-*{Cxl + 1 }* ioi. y*(z* +C) + 3 z* =0 

102. y = {C sec* z + 2 sec z tan z}103. y ={Ce x + x + 1 }-* 
104. y* = (C+z 3 -z*)/(z + l) 105. y = 1 + l/(Ce I +x) 

106. z*y =x +1 + l/(Ce -2/a - J) 107. y =z + l/(( 7 z + az™) 

108. y=z + l/(Ce- x +z-l) 109. y =x +o/{C'V( a; * +«) + 2:c} 
no. y =z(Cz+sin z)/(Cz-sin z) in. a; sin (x+y)=C 
ii2. y= x tan ( x+C ) 

“3- y 8 =(l +log | 1 -a; | -log | 1 +x |} + 2 x 

114. ( 1 +z)y s =C-z*+z* 1x5. (x +y)* + 2 =C(x -y)‘ 

1x6. (z* + y*)* + 2 a(z*-y*) =C 117. l-o/r*=Ce Jafl 
xi8. log r*+a/r* =C-f cos 20 119. z* =y , (C*-y 1 ) 

120. y=z* cob (z+C) X21. y = (x* +Cx)l(Cx* + 1 ) 

122. Equation becomes z"z' = 2 z (i) y=z coth (C-z), 
(ii) y=z(l +Cz *)/(1 -Cx*), (iii) y=z coth (z" 1 +C) 
X23. xy — G X24. 2 z*+y* = C* 

125. z + G =0 log {a + •/(<** - y*)> - o log y - V(<* 8 -S/*) 

126. y = Cz* + 2 a*/ 3 z 127. y =0 cosh {(z - C 7 )/a} 

128. zl+y*=C 129. y* —Ct x,a + 4 a(z +0) 

130- y=Cz+az log | z | - 6 131. (o — y)(2a + y)* = 9 a(( 7 ±z)* 

132. ay = C + £zv'( 4 z* - a*) - |a* log { 2 z + - a*)} 

133 * r=Ce ±ke 134. r = C sin 0 and r =C cosec 0 
135. r=a sin (0+C) 136. r=0+p log | tan $0 | 

*37* (0 +C)* = 1 +2a/r 138. r“ = C sin m0 

139 - y 3 - 3 z*y = p 140. y 1 + 2 x*=p 

141. z*+y»-2a log | z | =p 142. (z* +y*)* =p + 2a*(z* -y*) 

143. y( 3 z* + y») = / 3 (z* + y*)» 144. x«(y* - &)‘ = /?(** - a) V a 

145. y*+z*- 4 z+log(z + l) 4 = /l 146. 2 y*+z*(logz*-l)=/S 

147. z +a + 1 =P&xp (§y* +z) 148. y=coshz+/l 

149. Self-orthogonal 150. 3 a(y* -z*) + 2 z* =P 

151. z* + y* = 2a*log | /?z | 152. a*logz* +6*logy* =/? + z* + y* 

153 - V( 2 r/J-l)-arctan V( 2 r/i - 1 ) = P ±\6 

154. r - \jr = p ~ 0 i 55 . (r* - a*) cos 0 + /?r = 0 

156. r J (sin 0 - | sin 20) =p IS7 . r = p % /(sm 0) sin $0 

158. z* -1- 2 mzy -y'=P X59. z* + y* = 2 /?(z - ;ny) 

160. r v = p cos (pO + oj) 

x6x. y =G\ x | 1/m -z*/( 2 m - 1 ) (ro*$), y =z*(C -logz*) (w = J); 
m > 0 node, m < 0 col 

162. y* =Ce ,x +x + }, limiting point 163. z*-2Cy=C* 

164. e** - 2 Ce* + C* cos* z = 0 
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x 6c 144 ( 2 y-x*-x +C )*=(1 +8x) * 

•£ «*+“^ + v=; s y=s i”- ^''Crf-oshx 

7 y3 =Cx + 3 /C [Singular Solution y # = 12 x] 

Si. y=Ce* + \IC [S.S. y‘« 4 e-] 

179. (x - C)» + y* =f(C ); also S.S 
tRo 27(x* + y 1 ) 1 - 4x(x* + 9y*) + 4 y -0 
181* 2* - 2 C*xy +C* =0 [S.S. 32 xy* = 27 ] 

Is*. (x, + C)-= 4 (v«W-Cx) 

183 . y +C = e , “ +c,/s [S-S. y =x +a] 

184. (*+C)--4o(y+C)lS.S.y=* -«] 
fir- nr?# + ll* = 4 C’x* [S.S. y+x- 0 ] 

,86. (2x+C)' + ( 2 x +C)' - 18 ( 2 x +C)y = 27 ,' + 16 y 

187. y* =2Cx +C 7 * 

l88 - ;ii pip- + 2P [S.8. y - f +« 

l8 ,.x(l-pi-=0 + Sp--p-. p(.-p)-=cp- + p--jp 

z!;.! 1 ., - 1 *, y =C + !<•* - log | < I - ‘ + J 1 ' + 1 ( * 

,9 °- l =x c± Uc) ... [S.S ,*b=±W« 

i * -** 

**• ■Ttw. ^i’l-Viog +1 )4 '’ 

. 

198. y=C+x + l(C-x) j C . [s .s. x- +»■="). 

Z CWes’enveloped by the con.o 

201. y =log I-^ e + / * 1 ^,,3 +a = Bx A 

z.r + A^y \+*« - 6 - 91 ^ cm2< “* Bsin2ai 

Tos. tut Int - ,yyT= ' ly ' 21 


M2 


DIFFERENTIAL EQUATIONS 


209. y n+1 =Ax n+1 + B [when n = -1, y=Ax B ] 

210. y = 1 /{Ax +Bx') 2zx. y =x(A +B log | x |)* 

212. y-A + B/x+{( 1 - x)/x } log | 1 -2 | 

213. (i) real roots m and n, y = {Ae~ maz +Be- nax )-‘ u *; 

(ii) n= 7 n,y=e m *(A + Bx)~ 1Ja i (iii) complex roots h ± hi, 
y=Ae kx | cos ka(x-B) | -1/o 

214. (sc-B)* + (y -B)*=A* 215. a = 1/(1 -m) 

216. y = l/(x* +Ax + B) 217. y = l/\/(A cos x +B sin x +3) 
218. y =.4 exp (Be*) 219. y =(Ae tx +Bc ~* x ) 3 

220. A*y* -log y* = 44(2+5) 221. Ay* + (Ax -B)* =k 

222. y =a cosh (A +x/a) +B 223. A i y* = l+Be ax 

224. x =A(t + $ sin 2 1 ) +B, y-A sin* t (p= tan t) 

225. r =A cos 0 +B sin 0 +sec 0 

226. (1 -x')y’ +2(m - l)xy' - 1 )t/ =0 

227. y=4(x + 3) + Z?e x (2* -42 + 6) 

228. y =>1(1 +x*) +B{x +(1 +x*) arc tan 2} 

229. y = (6x* - 62 + 1 )(^4 +B log | 1 - l/x |) +B( Ox - 3) 

230. y =(A + Bx n+l )/(l +2) 

231. y = Ax/( 1 - 2) + B{x + 1 + 2x log | x |/( 1 - 2)} 

232. y =Ae~* + Be -* 1 + \c x + T l 5 sin 2 - ^ cos 2 
2 33* y =Ae~ 1 + Be 4 * - £ cos 2x - g sin 2x 

234. y =e~ x (A cos 2 x + B sin 2x - 1) + $e* 

235. y=e~* T (A + Bx + cosh 2) 

236. y =Ae ir + Be 3x +e z ( 2 x* + 62 +7) 

237. y =e~ x (A cos x\/2 + B sin xy /2 + cos 2) 

238. y =e tx (A +Bx + $2* cos 2x) 

239. y=Ae-* r +e x (B +Cx+ J.x* - J.x* + 3 , (J 2 4 ) 

240. y =e x8,n '{4 cos (2 cos a) +B sin (2 cos a)} +cot a cos 2 

241. y = e* ta “ '{.4 cos 2 + B sin 2 + (2* - 2) sec* a} 

242. y =(4 +Bx - 3x*) cos 2 + (4, +B x x -x*) sin 2 

243* y=(A +Bx + Jx*)e* + (>1, + B l x)e~* + C cos 2 +(7 X sin 2 
244. - T ' 0 c r (sin 2x + 2 cos 2 x) 245. $e ,r (cos 2 + sin 2) 

246. e-'( jx* -2) 247. 2x sin 2 - £ cos 32 

248. 2 sin 2 - 2* cos 2 249. jx* cos 2 250. 2 s - 62 - 12 

251. e J (Ax s + + |x 4 ) 252. - e x log | 2 | 

253. - e 1 sin 2 254. - 25e* cos 2 + e _, (3 cos 2 - 4 sin 2) 

255. e ir ( 2 x* - 4x + 3) 256. y =cos mx + 2m(e* - 1) sin mx 

257. y = e~ x (x -l) s 258. y = (2 + 1 )e ax sin mx 

259. !/ = (2 - 1 )e’ + (2 + 1 )*e - ' 260. y=e~ x + Je*(2 - 1)* 
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261. 

262. 

263. 

264. 

265. 

266. 

267. 

268. 

269. 
271. 
273 - 
275 - 
276. 

277 - 

278. 

280. 

281. 

282. 

283. 

284. 


285. 

286. 

287. 

288. 

289. 

290. 

291. 


292. 

293 

294 

295 


y=Ax + BVx + l x ' 

y=Ax*+x{B- log x*-(log | x I) 1 } 

y =A/x + Six' + i( 1 + 1/*)’ lo S I 1 + r 1 " * . 

y =4 cos (log 1 X I) + B sin (log | 2 |) + log | x I 

y = Ax' + Bx - x sin 2 

1/ = {x4 + log | 2 + 1 |)(2 + 1)+BV(* + 1)-1 

y =2 a {A +B log | 2 1} +2 1 - a /(2a 1) 

2/ = B log | x | +i(log I *_l) > n x +fie _ sln , 

5 = ^(1-*»)+Bx>/(l+*‘> 274 .<-JW*» 

y =>1 COS V* +B 6 > n V x 

y =Acost+Bslni (t= log | tan x |) 

4 =exp { - ijpd*); r =9 - ip 5 - m 

y =x n (A cos r/12 + Z? sin mx) 279. 2/ 

y=/rinx + B(cosx + xsinx)-i<x + l)cosx 

x=Ae'+Be~' + le •» 
v + 2po = 0, also (a) <7 < 0, (6) 9 > 0 
(a) 9= -m*/x*; *=*" and x-; (6) 9 -j» /* • 

1/ = cos (m log I X 1) and sin (m log | x 1) m 

y = Je*'*‘ , " 1 (<Terminals (0, 1), also ( - . 

^ if x > 0 or (1, « ) if x < ° 

bin 2x log 1 cos 2 | -x cos 2x 

log 1 cos 2 | +sin x log | tan (J* + **) 1 “ 

-:"■+<*&::I2&: y* 

/M *x = 1 (M^e-‘ /4 +V e '* ,/6 + 21 8in 1 ' 67 003 

_ »* a c -*' /5 - 13 sin t + 18 cos f) 

X = fl -V 7 ) cos I. If = 3(1 - !•) Bin * + 14 t cost 

a . = 2 c’ , - 2 cosf- 3 f sin t, y=<cos*+(l 
y^l-xNy^x-^-*/-^’- • * ‘ 
with (r + 3 )C' r+3 — ( r _ 3 ) r 

y =l+22*+iX*-/ 7 X a +AB® ' * * 

with ( 2 r - \)G r = - ( 2 r - 10 )C,_». 
y 1 = Vx{l+S* , +.1 B 0* 4 -TS6* a+ * * *' 
with 2 rC, = - ( 2 r - 9 ) 0 r-* 
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296. u = 1 +£* +x i l( 2 \) i + x 6 /( 3 !)* + . . . 

v=** + (l+|)^/( 2 !)*+(l+|+J)x*/( 3 !) , + . . . 

a a(a + 1) 

297 ' F <“ ; ” X) = 1+ y X + 2!y(y + lf ’ 

a(a + l)(a+2) s 

+ 3 ! y(y + l)(y + 2 ) x 

298. F(y -a; y; -x) 

299. F( a ; 1 ; x) log | x \ + I>C,x r , where 

_/l _1_ 1 2 2 

\a + a + l + + a+r-l 1 2 

a(a + 1). . .(a +r - 1) 


• • 


(r! ) 


(r>l) 


300. S( - )'(2«)! (**)*•/{(« l)*I> + n + 1)I> -n + 1)} 




INDEX 


Asymptotic point, 39. 

Auxiliary equation, 85. 
complex roote, 87. 
repeated roots, 88. 

Base point, 37. 

Bernoulli equation, 22. 

Bessel equation, 124. 

Bessel function «/«(*)» 125. 

Y n (x), 126. 

Change of variable, 24. 

Clairaut equation, 42. 

Col, 36. 

Complementary function, 19, »i. 
Condenser, charging of, 20. 
Constant coefficients, linear 
equation with, 82-101. 
linear system with, 108-111. 
Curvature, 76. 

Cusp-locus, 63, 66. 

d’Alembert equation, 43. 
Definite integral, solution by, 
104. 

Degroo, 2, 40. 

equation of second, j 5. 
Discriminant locus, 62. 

Envelope, 63, 67. 

Euler linear equation, 7 /, iui. 
Exact equations, 10. 

Families of plane curves, 27. 
Frobonius, method of, 119. 

General integral, 3, 40, 02. 

by quadratures, 100. . 
Geometrical interpretation, 50. 


Homogeneous types, first order, 
6 , 48. 

second order, 67-72. 
Hypergeomotric equation, 116. 

Indicial equation, 116. 

Infinity, singularity at, 124. 
Initial conditions, 21. 
Integrability, condition of, 11. 
Integral curves, 27-39. 

Integral, first, 73. 

Integrating factors, 13, 73. 
quotient of, 15. 
special types, 16. 

Integration, 3. 

Isotropic nodal point, 37. 

Lagrange equation, 44. 

Laplace equation, 104. 
Legendro equation, 119. 
large |x|, 122. 

Legondre function Pn( x )» 120. 

(?„(*). 123. 

Lovol lines, 34. 

Limiting point, 38. 

Linear coefficients, 8. 

Linear equation, first order, 18 

21 . 

in gonoral, 79-107. 

Lincur fractional equation, JO. 
Linear indopondonco, 80. 
Linear systems, 108. 

Lines of slope, 34. 

Nodal point, 37. 

Normal, 28, 76. 

Oporator, linoar, 86 . 
inverse, 90-99. 
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DIFFERENTIAL EQUATIONS 


Order, I. 
of a eystem, 109 . 
reduction of, 01, 81 . 
Ordinary point, 113 . 
Orthogonal trajectories, 30 . 

Parametric integration, 63 . 
Particular integral, 3 , 19 , 81 . 
p-discriminant, 64 . 
Permutable operators, 83 . 
Polar co-ordinates, 25 . 

trajectories, 32 . 

Primitive, 3 , 62 . 

Principal directions, 36 . 

Quadratures, 6, 19 , 100 . 

Recurrence relation, 114 . 
Reduced equation, 19 , 80 . 
constant coefficients, 82 - 90 . 


Regular singularity, 116 . 
Riccati equation, 22 . 

Rodrigues formula, 121 . 

Second degree, equation of, 56 . 
Separation of variables, 4 . 
Series solution, 112 - 129 . 
Simultaneous equations, 108 - 
111 . 

Singular integrals, 3 . 

Singular points, 35 , 116 . 
Subnormal, 28 . 

Subtangent, 28 . 

Tac-locus, 68. 

Tangent, 28 . 

Taylor series, 112 . 

Trajectories, 29 - 35 . 

Variable absent, 47 , 62 , 65 . 
Variation of parameters, 100 . 





